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This  volume  Is  part  of  the  12-voluroe  report  entitled  "Hydrologic  i 

Engineering  Methods  for  Water  Resources  Development,"  prepared  by  The  j 

Hydrologic  Engineering  Center  as  part  of  the  U.S.  Amy  Corps  of  Engineers'  | 

participation  In  the  International  Hydrological  Decade.  Volume  12 
addresses  the  topics  of  river  morphology,  data  collection  and  analysis, 
reservoir  sedimentation,  and  aggradation  and  degradation  In  free 
flowing  streams.  The  emphasis  of  the  volume  Is  on  practical  approaches 

l 

and  techniques  for  analyzing  sediment  problems.  Although  many  of  the  > 

methods  and  procedures  described  herein  have  been  used  successfully  In 
Corps  of  Engineers'  studies,  the  volume  should  not  be  construed  to 
represent  the  official  policy  or  criteria  of  the  Corps. 
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CHAPTE:^  1.  INTRODUCTION 


Section  1.01.  Objective 

The  title  of  this  volume,  "Sediment  Transnort,"  was  chosen  to  focus 
attention  on  the  movement  of  sediment  material  in  flowinq  water.  This 
involves  processes  of  scour,  transport  and  deposition  of  inorganic 
material  both  in  free  flowing  streams  and  in  reservoirs,  '/hile  some 
sediment  transport  formulas  are  included,  they  do  not  form  a major  part 
of  this  volume. 

The  purposes  of  this  volume  are:  (1)  to  identify  potential  problem  , ] 

areas;  (2)  to  identify  which  of  these  can  be  analyzed  with  existing 
mathematical  techniques  and  which  must  be  studied  using  either  movable 
bed  hydraulic  models  or  the  prototype;  (3)  to  identify  the  type  and 
amount  of  data  reouired  for  analyzing  sediment  problems,  and  (4)  to  give 
sufficient  guidance  so  that  comoetent  engineers  may  develop  satisfactory 
solutions  to  the  sediment  problems. 

Section  1.02.  Scope 

Basic  to  all  sediment  studies  is  sediment  yield  from  the  watershed. 

The  paoer  entitled  "Corps  of  Engineers  ’'ethods  for  Predicting  Sediment 
Yield"  by  ’*r.  Robert  H.  Livesey  is  included  as  Appendix  III  to  provide 
i information  in  this  problem  area.  Also,  the  calculation  of  land  surface 

I erosion  is  presented,  but  only  briefly  described.  i 
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Thp  basic  concept  of  calculating  water  surface  profiles  in  natural 
rivers  is  exten^e'1  bcvon'i  ''olune  6,  "Water  Surface  Profiles,"  to  include 
the  river  bed  as  a novable  boundarv.  The  analysis  of  sediment  oroblens 
bv  digital  comuter  simulation  is  useful  for  calculating  the  volume  and 
location  of  sediment  deposits  in  reservoirs  and  for  predicting  aggra- 
dation or  degradation  trends  downstream  from  dams  as  v/ell  as  in  free 
flov/ing  streams.  The  Hydrologic  Engineering  Center  comouter  program 
'/hich  performs  the  simulation  analysis,  "Scour  and  Deposition  in  Rivers 
and  ’’eservoirs ,"  is  included  in  Annendix  VII  of  this  volume.  This  program 
can  be  used  to  evaluate  the  behavior  of  a stream  bed  during  the  passing 
of  a single  flood  event  or  for  a long  oeriod  of  hydrogranhic  record. 

Technigues  are  also  presented  in  Appendix  VII  for  calculating  growth 
of  the  armor  laver,  destruction  of  the  armor  layer  and  hydraulic  sortinn 
of  grain  sizes. 

''  very  useful  technigue  for  calculating  the  volume  of  sediment 
deposits  in  reservoirs  is  based  on  trap  efficiency.  Application  of  the 
technioue  does  not  renuire  the  electronic  computer. 

Since  present  analytical  technigues  do  not  completely  define  sedi- 
ment problems,  a good  understanding  of  the  behavior  of  natural  and 
controlled  streams  is  essential  for  interpreting  any  calculated  results. 

''  brief  discussion  of  river  morphology  is  included  to  aid  in  this  regard. 

In  summary,  attention  is  directed  to'.-.'ard  the  analysis  of  v/ater  born 
sediments  in  defined  channels.  The  regime  concept  is  not  treated  in 
detail.  Hind  blc./n  sediment  is  not  included.  The  impact  of  sediment  on 
v;ater  gualitv  - that  is,  guality  constituents  that  might  adhere  to  and  be 
transported  v/ith  sediment  oarticles  - is  discussed  in  Volume  11.  The 
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grov/th  and  decay  of  bed  foms  and  tv/o  and  three  dimensional  hydrodynamic 
models  are  not  included  in  this  volume.  Coastal  processes  are  not 
addressed. 

In  general,  the  formulas  which  are  included  have  been  incornorated 
into  analytical  methods  and  have  been  anplied  to  a wide  ranee  of  problems 
throughout  the  United  States.  Their  oerformance  has  been  satisfactory 
when  the  methods  were  applied  v/ith  the  appropriate  understanding  of  river 
morphology  and  the  limitations  of  present  theory.  Even  the  most  advanced 
methods  require  a great  deal  of  engineering  judgment  to  insure  a correct 
interpretation  of  results. 

’'lot  all  analytical  techniques  in  use  in  the  Corps  of  Engineers  are 
included  in  this  volume.  This  does  not  imoly  that  techniques  vihich  are 
included  are  the  "best"  ones,  nor  does  it  imply  that  techniques  which 
v/ere  not  included  are  inferior.  It  simply  reflects  the  fact  that  there 
are  a large  number  of  methods  in  use,  there  is  only  a limited  amount  of 
space  available  to  present  them,  and  the  methods  presented  have  been 
found  to  be  generally  satisfactory. 

Section  1.03.  A Summary  Statement 

nature  maintains  a very  delicate  balance  between  the  v/ater  flov/ing 
in  a natural  river,  the  sediment  load  moving  with  the  v;ater  and  the 
stream's  boundary,  ''.ny  activity  of  man  which  changes  any  one  of  the 
follo'/ing  parameters  i 

• water  yield  from  the  v/atershed 

• sediment  vield  from  the  watershed 
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• water  discharge  duration  curve 

• depth,  velocity,  slope  or  width  of  the  flow 

• size  of  sediment  particles, 

or  which  tends  to  fix  the  location  of  a river  channel  on  Its  floodplain 
and  thus  constrain  the  natural  tendency  to  meander,  upsets  the  natural 
balance  and  Initiates  the  formation  of  a new  one.  The  objective  of  most 
sediment  studies  Is  to  evaluate  the  Impact  on  the  flow  system  from 
changing  any  of  these  parameters. 

A classical  example  Is  the  deposition  which  results  when  a reservoir 
Is  Impounded.  In  terms  of  the  aforementioned  parameters,  the  reservoir 
causes  a change  In  the  hydraulics  of  flow  by  forcing  the  energy  gradient 
to  approach  zero.  This  results  In  a loss  of  transport  capacity  with  the 
resulting  deposition.  The  smaller  the  particles,  the  farther  they  will 
move  Into  the  reservoir  before  depositing.  At  times  fine  sediment  In  the 
Inflow  to  a deep  reservoir  does  not  fully  mix.  The  resultant  stratification 
Is  conducive  to  the  formation  of  density  currents,  and  some  material  may 
pass  through  the  reservoir. 

The  obvious  consequence  of  sediment  deposits  Is  a depletion  In  reser- 
voir storage  capacity.  This  Is  represented  schematically  In  fig.  1.01, 


Deposition  in  Deep  Reservoirs.  The  volume  of  sediment  material  in  the 
delta  is  a function  of  the  project  life  among  other  things,  and  that 
delta  will  continue  to  develop,  with  time,  beyond  the  project  life. 
Eventually  a new  channel  and  floodplain  will  exist  in  the  delta  area. 

Identifying  this  deposition  pattern  with  deep  reservoirs  raises  the 
question  "what  is  the  significance  of  deeo  as  ooposed  to  shallow  reser- 
voirs on  the  sediment  deposition  pattern?".  Referring  to  fig.  1.01, 
consider  v/hat  would  happen  if  the  dam  v/ere  moved  in  the  direction  of  the 
delta  denosits  and  the  operating  rule  for  the  reservoir  did  not  change. 
The  configuration  of  the  delta  would  not  change  until  the  dam  actually 
reached  the  toe  of  the  delta  itself.  At  some  point,  the  reservoir  could 
become  classified  as  a shallow  reservoir.  Therefore,  the  delta  formatio!! 
shown  in  fig.  1.01  could  be  as  applicable  to  a shallow  reservoir  as  it  !s 
to  a deep  reservoir. 

The  mean  operating  pool  elevation  is  a major  factor  in  establishing 
the  configuration  of  the  reservoir  delta.  Therefore,  to  shift  the  dam  to 
a location  along  the  present  topset  slope  would  completely  change  the 
mean  operating  pool  elevation  and  thereby  the  shape  of  the  reservoir 
delta  itself. 

A number  of  associated  problems  can  be  listed: 

• Deposits  forming  the  delta  may  raise  the  water  surface 
elevation  during  flood  flows  in  such  a manner  to  require 
special  consideration  for  land  acquisition.  In  deep 
reservoirs,  this  is  usually  not  a problem  within  the 
reservoir  area  because  project  purposes  dictate  land 
acouisition.  Mov/ever,  the  delta  tends  to  develop  in  the 
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upstream  direction  and  cause  problems  upstream  from  the 
reservoir  area  proper.  In  shallow  reservoirs,  on  the 
other  hand,  the  Increase  In  water  surface  elevation  Is 
a problem  even  within  the  reservoir  area.  That  is, 
floods  of  equal  frequency  may  have  higher  water  surface 
elevations  after  a project  begins  to  develop  a delta 
deposit  than  was  experienced  before  the  project  was 
constructed.  Land  acquisition  studies  must  evaluate 
such  a possibility. 

Aggradation  problems  are  often  more  severe  on  tributaries 
than  on  the  main  stem  and  these  locations  are  often  desirable 
for  developing  recreational  facilities.  Analysis  Is  complicated 
by  the  lack  of  basic  data  on  the  tributaries  - usually  less  than 
on  the  main  stem  itself.  However,  the  useful  life  of  recreation 
sites  should  be  evaluated  by  predicting  the  rate  of  delta  growth. 

Because  of  the  high  moisture  level,  reservoir  deltas  often 
attract  phreatophytes  which  In  many  areas  contribute  to 
water  use  problems  due  to  their  high  transpiration  rate. 

Reservoir  delta  deposits  are  often  aesthetically  undesirable. 

Particularly  In  shallow  Impoundments,  the  results  of  reservoir 
sediment  deposits  may  Increase  the  water  surface  elevation 
sufficiently  to  Impact  on  the  ground  water  table. 


• In  existing  reservoirs,  the  United  States  Fish  and  Wildlife 
Service  is  utilizing  delta  and  backswamn  areas  in  the 

nroD  gation  of  v/ildlife.  Since  the  characteristics  of  this 
delta  area  are  so  closely  controlled  by  the  operating  policy 
of  the  reservoir,  any  reallocation  of  storage  would  need  to 
consider  the  impact  on  present  delta  and  backswamn  areas. 

This  represents  a type  of  problem  that  may  be  more  important 
in  the  future  if  changing  priorities  among  oro.iect  purposes 
demand  reallocation  of  storage  in  reservoirs. 

• Looking  do'/nstream  from  the  dam,  degradation  is  usually 
predominant.  It  is  necessary  to  evaluate  the  impact  of 
deoradation  on  a tailwater  rating  curve  at  the  structure. 
Problems  downstream  from  the  dam  are  sufficiently  comolex 
that  thev  are  presented  in  some  detail  in  the  following 
naragranhs. 

Do'.'nstream  from  the  dam  the  hydraulics  of  flow  (velocity,  slooe, 
deoth  and  v.-idth)  remain  unchanged  from  conditions  in  the  natural  state. 
However,  the  reservoir  has  acted  as  a sink  and  trapped  sediment  material, 
esneciallv  the  bed  material  load.  This  reduction  in  sediment  yield  from 
the  watershed  causes  the  energy  In  the  floiv  to  be  out  of  balance  with  the 
boundary  material  for  the  dov/nstream  channel.  Because  of  the  amount  of 
available  energy,  the  water  attempts  to  re-establish  the  former  equili- 
brium sediment  load  from  material  In  the  stream  bed,  and  this  results  In  a 
degradation  trend.  Initially,  the  degradation  trend  may  persist  for  only 
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a short  distance  downstream  from  the  dam  since  the  equilibrium  sediment 
load  Is  soon  re-established  by  removing  material  from  the  stream  bed. 

As  time  passes,  the  degradation  trend  will  migrate  downstream; 
however,  several  factors  are  working  together  to  establish  a new  equili- 
brium condition  In  this  movable-boundary  flow  system.  In  addition  to  this 
re-entrainment  of  sediment  material.  On  one  hand,  the  potential  energy 
gradient  Is  decreasing  because  the  degradation  migrates  from  upstream  to 
downstream  In  direction.  On  the  other  hand,  the  bed  material  Is  becoming 
coarser  and,  consequently,  more  resistant  to  being  moved.  This  tendency 
In  the  main  channel  has  the  opposite  effect  on  tributaries.  Their  poten- 
tial energy  gradient  Is  Increasing  which  results  In  an  Increase  In  trans- 
port capacity.  This  will  usually  Increase  sediment  passing  Into  the  main 
stem  which  tends  to  stabilize  the  main  channel  resulting  In  less  degrada- 
tion than  one  might  anticipate.  Finally,  a new  equilibrium  condition 
will  be  established  between  the  flowing  water-sediment  mixture  and  the 
boundary. 

The  extent  of  degradation  Is  complicated  by  the  fact  that  the  reservoir 
also  changes  the  water  discharge  duration  curve.  This  will  Impact  for 
great  distances  downstream  from  the  project  because  the  existing  river 
channel  reflects  the  historical  phasing  between  flood  flows  on  the  main 
stem  and  those  from  tributaries.  That  phasing  will  be  changed  by  the 
operation  of  the  reservoir.  Also,  the  flow  will  probably  encourage  vegetation 
to  grow  at  lower  elevations  In  the  channel.  The  result  Is  a condition 
conducive  to  deposition  In  the  vegetation.  Actually,  numerous  examples 
of  aggradation  trends  may  be  cited  for  river  channels  downstream  from 
dams.  The  primary  problem  which  results  Is  Inadequate  channel  capacity 

1-08 


resulting  In  Inadequate  levee  height.  Consideration  should  be  given  to 
performing  detailed  simulation  studies  to  determine  future  channel 
capacities  and  to  Identify  problem  areas  of  excessive  aggradation  or 
degradation.  Particular  attention  should  focus  on  all  major  tributaries. 

An  equally  classical  problem  results  when  levees  are  constructed 
for  flood  control  purposes  because  the  position  of  the  river  channel 
must  be  stabilized  to  Insure  that  It  remains  between  the  levees.  River 
flow  entering  the  backwater  curve  to  Its  receiving  water,  whether  that 
receiving  water  Is  the  ocean  or  a lake,  loses  transport  capacity  and 
deposition  of  the  coarser  material  results  In  an  aggradation  trend  In 
the  river  channel.  Consequently,  the  levee  height  becomes  Inadequate 
to  contain  the  design  discharge.  The  rate  at  which  this  process  takes 
place,  although  slow.  Is  not  measured  In  terms  of  geological  time. 
Depending  on  the  sediment  yield  In  the  basin,  significant  changes  can 
occur  during  the  50-100  year  life  of  an  engineering  project. 

A third  classical  sediment  problem  Is  that  of  maintaining  depth  of 
navigation.  This  problem  requires  a detailed  understanding  about  the 
behavior  of  the  movable  boundary  flow  system  since  It  requires  fixing  the 
location  of  a navigation  channel  within  the  main  channel  Itself.  However, 
the  techniques  for  performing  these  studies  can  proceed  under  the  assump- 
tion that  the  location  of  the  river  channel  Is  fixed.  Hydraulic  model 
studies  can  then  be  employed  to  design  the  navigation  channel. 

The  aforementioned  sediment  problems  are  all  associated  with  man's 
activities.  However,  left  to  Its  natural  state,  a river  will  continually 
change  Its  position  on  the  flood  plain.  Its  meander  pattern  and  the  cross 
section  shape  as  It  responds  to  the  flowing  water  and  sediment  mixture. 
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A good  understanding  of  this  process  Is  required  before  one  can  adequately 
Interpret  the  Impact  of  man's  activities,  for  man's  activities  not  only 
change  the  rate  at  which  a river  channel  responds  to  the  flowing  water- 
sediment  mixture,  but  also  change  the  location  of  sediment  sources  and 
sinks  from  those  naturally  existing. 

The  level  of  detail  required  for  the  analysis  of  any  sediment  problem 
depends  on  the  objective  of  the  study.  Considering  a dam  site  as  an  Im- 
portant natural  resource.  It  Is  essential  to  provide  enough  volume  In 
the  reservoir  to  contain  anticipated  deposits  during  the  project  life.  If 
the  objective  of  a sediment  study  is  just  to  know  the  volume  of  deposits 
for  use  In  screening  studies,  then  trap  efficiency  techniques  provide  a 
satisfactory  solution.  The  Important  Information  that  must  be  available, 
then,  is  the  water  and  sediment  yields  from  the  watershed  and  the  capacity 
of  the  reservoir.  However,  If  the  sediment  study  must  also  address  the  land 
acquisition  for  the  reservoir,  then  knowing  only  the  volume  of  deposits  Is 
not  sufficient.  The  location  of  deposits  must  also  be  known,  and  the 
study  must  take  Into  account  sediment  movement.  The  approach  presented 
In  this  volume  Is  to  analyze  such  a problem  with  digital  simulation  of  flow 
In  a mobile  boundary  channel.  Sorting  of  grain  sizes  must  be  considered 
since  the  coarser  material  will  deposit  first,  and  armoring  must  be  accommo- 
dated since  scour  is  involved.  This  simulation  technique  actually  becomes 
a movable-bed  analytical  model.  It  Is  useful  to  predict  erosion  or  scour 
trends  downstream  from  a dam,  general  aggradation  or  degradation  trends  In 
river  channels,  and  the  ability  of  a stream  to  transport  the  bed  material 
load.  Such  a technique  does  not  have  to  define  the  location  of  a navigation 
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channel  In  order  to  be  useful;  in  fact,  there  is  no  analytical  technique 
that  IS  suitable  for  calculating  where,  within  the  river  channel,  a 
navigation  channel  will  be  stable. 

The  simulation  technique  Is  structured  entirely  for  computer  solution. 
The  amount  of  data  that  has  to  be  analyzed  Includes  all  the  basic  geometric 
and  hydraulic  data  required  for  water  surface  profile  calculations  plus 
data  describing  the  size  and  gradation  of  sediment  material  In  the  stream 
bed  and  banks,  the  size  and  amount  of  Inflowing  sediment  material  and  the 
water  discharge  hydrograph.  In  addition,  long  periods  of  hydrograph 
record  are  generally  utilized  since  sediment  studies  attempt  to  predict 
trends  throughout  the  project  life.  The  number  of  calculations  Is 
extremely  large.  For  example,  predicting  deposition  In  a shallow 
reservoir  having  a 50  year  design  life  can  require  a calculation  of 
some  600  to  1000  water  surface  profiles  plus  the  routing  of  sediment 
material  through  the  reservoir  for  the  water  discharge  associated  with 
each  one  of  the  profiles. 

Data  acquisition  programs  have  evolved  to  satisfy  needs.  With  only 
few  exceptions,  data,  when  collected  at  all,  are  collected  for  estimating 
annual  sediment  yield  where  suspended  sediment  Is  the  primary  contributor. 
This  Is  definitely  one  problem  area.  However,  there  are  other  problem 
areas.  For  example,  data  should  also  be  collected  In  a systematic  program 
to  study  channel  morphology.  This  would  require  measuring  the  bed  load 
to  determine  the  amount  moving,  the  gradation  of  the  load  and  the  sorting 
that  goes  on  among  the  various  grain  sizes.  The  link  to  understanding 
the  feedback  mechanism  In  the  water-bed  Interaction  probably  lies  In 


understanding  bed  load  and  the  lateral  movement  of  river  channels.  Attempts 
to  simulate  the  behavior  of  a movable  bed  analytically  will  continue  to 
require  considerable  qualitative  judgment  until  this  link  Is  understood. 

Even  present  techniques  require  data  on  gradation  of  material  In  the  stream 
bed  and  In  the  sediment  load.  Sampling  the  bed  load  will  be  required  to 
obtain  this  type  of  data.  These  data  should  be  collected  by  personnel 
Intimately  familiar  with  both  sediment  movement  and  analytical  techniques 
until  a better  understanding  of  the  nature  of  the  problem  Is  developed. 

All  data  acquisition  programs  should  utilize  "standard"  equipment  and 
techniques. 

Section  1.04.  The  Interpretation  of  Analytical  Results 

The  analysis  of  sediment  problems  Is  not  a simple  extension  of  fixed 
bed  hydraulic  theory  so  that  It  becomes  movable  boundary  hydraulic  theory. 
Sediment  problems  vary  In  degree  of  difficulty  from  the  relatively  simple 
determination  of  the  volume  of  deposits  In  a deep  reservoir  to  aggradation 
and  degradation  studies  In  free  flowing  streams  and  rivers.  Whereas  fixed 
bed  hydraulic  theory  Is  well  developed,  the  analysis  of  movable  bed  problems 
Is  complicated  by  the  fact  that  the  body  of  theory  available  for  performing 
analyses  Is  Incomplete.  The  Interactions  between  flowing  water  and  a 
movable  boundary  are  not  well  understood,  although  the  water  and  boundary 
are  components  of  a closed  loop  system.  At  best,  the  available  theory 
addresses  only  bits  and  pieces  of  that  system. 

It  Is  often  uncertain  as  to  how  to  apply  the  available  theory  to 
obtain  satisfactory  results.  Guidelines  are  conspicuously  absent.  Case 
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after  case  1s  presented  to  demonstrate  the  Inconsistency  of  results  that 
one  can  obtain  for  the  same  problem  analysis  by  using  different  methods 
— all  of  which  are  recognized  in  the  literature  as  acceptable  methods. 

Nevertheless,  the  engineer  is  constantly  faced  with  the  analysis  of 
sediment  problems.  Perhaps  the  following  suggestions  will  be  heloful. 

It  is  good  oractice  to  follow  a three  step  procedure:  (1)  calibrate  the 
analytical  technique,  (2)  verify  the  procedure  by  performing  a "base  test" 
for  comparison  with  observed  conditions,  and  (3)  interpret  the  imnact  of 
any  changes  by  reference  to  that  base  test  rather  than  by  absolute  magni- 
tudes. Due  to  the  difficulty  of  obtaining  consistent  prototype  data 
and  the  uncertainties  involved,  good  engineering  judgment  is  required  in 
determining  the  source  of  discrepancies  between  calculated  results  and 
measured  values.  This  is  especially  important  in  simulation  of  the 
movable  boundary,  but  the  same  approach  is  valuable  even  in  trap  efficiency 
studies. 

Another  useful  technique  is  to  select  several  formulas  for  use  in 
each  phase  of  the  study  and  compare  results.  Sensitivity  studies  will 
often  help  to  identify  the  most  appropriate  formula  for  the  problem  at 
hand. 
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CHAPTER  2.  TERMIUOLOGY  AND  PROPERTIES  OF  SEDIMENT 


Section  2.01.  Introduction 


The  terminology  In  sedimentation  work  Is  somewhat  unique.  Some 
terms  appear  to  be  general  In  nature,  such  as  bed  load,  or  bed  material 
load,  and  yet  they  have  been  defined  rather  precisely.  Other  terms 
appear  to  be  very  well  defined  and  yet  are  somewhat  general  In  nature; 
particle  size  classification  Is  an  example.  The  fact  that  particle 
size  Is  an  Important  variable  Is  obvious,  and  yet  there  are  several 
different  standards  for  classifying  sediment  material  according  to 
particle  size.  Many  of  these  use  the  words  “sand,  silt,  and  clay," 
but  particle  size  diameter  for  material  called  sand  Is  different  from 
one  classification  standard  to  another.  This  Is  also  true  for  silt 
and  clay. 

In  general,  the  physical  properties  of  sediment  are  Important  to 
both  quantity  and  quality  studies.  Consistency  within  the  United  States 
In  determining  these  properties  required  the  establishment  of  a special 
Interagency  sedimentation  project.  The  project  Is  under  the  auspices  of 
the  Committee  on  Sedimentation  of  the  Water  Resources  Council,  and  Is 
engaged  In  the  development  and  standardization  of  sampling  equipment, 
sampling  techniques,  and  methods  for  analysis  of  samples. 

It  Is  useful  to  organize  sediment  properties  In  three  groups: 

Properties  of  sediment  particles 

• Size  and  shape 

* Classification  scale 
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• Shape  factor 

• Specific  gravity 

• Fall  velocity 

• Gradation  of  sample 
Properties  of  sediment  deposits 

• Initial  unit  dry  weight 

• Consolidation  with  time 

• Unit  dry  weight  of  mixtures 
Properties  of  the  water-sediment  mixture 

• Concentration 

• Sediment  load 

• Sediment  yield 


Section  2.02.  Properties  of  Sediment  Particles 

a.  Particle  Size  and  Shape 

Particle  size  refers  to  the  “diameter"  of  a particle.  There  are 
several  ways  to  measure  diameter  but  the  most  comnon  are  by  sieving  or  by 
measuring  the  velocity  with  which  the  particle  falls  through  quiescent 
water  at  a standard  temperature.  The  latter  procedure  produces  an  equivalent 
spherical  diameter  — that  is,  the  diameter  of  a sphere  of  the  same  specific 
weight  that  would  have  the  fall  velocity  that  was  observed.  Sieving  is  the 
most  comnon  technique  for  measuring  the  sizes  of  sand  and  gravel  particles. 

Particle  size  is  probably  the  most  significant  physical  property.  It 
affects  the  resistance  to  erosion,  the  transportability  of  sediment,  and 
subsequent  behavior  in  the  consolidation  of  sediment  deposits.  For  example. 
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very  small  particle  sizes,  up  to  about  0.004  nm,  exhibit  a stronq  influence 
from  electrical  charges  on  their  surface.  They  are  said  to  be  cohesive, 
flaterial  having  this  range  of  grain  sizes  is  called  clay,  and  when  sediment 
problems  involve  clay  a special  body  of  analytical  theory  is  employed  to 
account  for  the  impact  of  electro-chemical  characteristics  of  the  water  on 
those  of  the  sediment  material.  Particle  sizes  between  0.004  mm  and  about 
0.0625  mm  are  in  a transition  range.  They  are  too  large  to  feel  much  influ- 
ence from  the  electromotive  forces  and  too  small  to  mobilize  much  inertia 
against  being  moved  by  flowing  water.  This  range  is  classified  as  silt. 

When  particle  size  exceeds  0.0625  mm  electromotive  forces  are  insignificant. 
These  particles  are  noncohesive  and  are  classified  as  sand,  gravel,  cobbles, 
etc.  Knowledge  of  mechanical  forces  is  sufficient  to  analyze  the  behavior 
of  noncohesive  sediment.  Consequently,  transport  theory  for  noncohesive 
material  is  more  advanced  than  that  for  cohesive  material. 

b.  Classification  by  Grain  Size 

The  grain  size  classification  scale  established  by  the  American 
Geophysical  Union  is  the  standard  used  to  relate  grain  size  to  size  class 
throughout  this  volume.  The  median  diameter  of  a size  class  is  that  for 
which  50*  of  the  material  in  that  class  is  finer  and  50*  coarser. 

Table  2.01.  Grain  Size  Classification 


Size  Class 

Clay 

Very  fine  silt 
Fine  silt 


Grain  Diameter 
(mm) 

less  than  .004 
.004  to  .008 
.008  to  .016 
.016  to  .032 


Median  Di ..meter 

fiSn) 


Medium  silt 
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.0057 

.0113 

.0226 


Table  2.01.  Orain  Size  Classification  (cont.) 


Size  Class 

Grain  Oiameter 
(rril 

Median  Diameter 
(mm) 

Coarse  silt 

.032 

to 

.0625 

.0447 

Very  fine  sand 

.Of  25 

to 

.125 

.0884 

Fine  sand 

.125 

to 

.250 

.1768 

f'edium  sand 

.250 

to 

.50 

.3536 

Coarse  sand 

.50 

to 

1.00 

.7071 

Very  coarse  sand 

1.00 

to 

2.00 

1.4142 

Very  fine  gravel 

2.00 

to 

4.00 

2.8284 

Fine  gravel 

4.00 

to 

8.00 

5.6569 

'ledium  gravel 

8.00 

to 

16.00 

11.3137 

Coarse  gravel 

16.00 

to 

32.00 

22.6274 

Very  coarse  gravel 

32.00 

to 

64.00 

45.2543 

c.  Particle  Shaoe  Factor 

Sediment  particles  are  seldom  spherical.  Clays  are  very  elongated 
and,  to  a lesser  extent,  so  are  the  larger  particles.  The  following  expres- 
<5  sion  has  been  proposed  as  a measure  of  particle  shape  by  using  dimensions 
normal  to  each  other: 

grain  shaoe  factor  = c/ Va*b  (2-01) 

where  a and  b refer  to  the  tm  smallest  dimensions  and  c to  the  largest. 

The  nrimary  influence  of  grain  shape  in  the  noncohesive  particles  is  on 
fall  velocity.  Whereas  this  is  recognized  in  conceot,  it  has  not  been 
formally  incorporated  into  a satisfactory  analytical  expression. 
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d.  Specific  Gravity 

The  specific  gravity  of  sediment  particles  is  another  property  that 
imnacts  on  the  fall  velocity,  and  consequently,  the  hydrodynamic  properties 
of  the  sediment  particle.  The  value  commonly  used  for  quartz  sand  is  2.65. 

e.  Fall  Velocity 

^ property  of  sediment  which  is  very  important  in  transport  calcula- 
tions is  the  velocity  at  which  a single  oarticle  would  fall  through  quies- 
cent water.  This  velocity  is  related  to  the  grain  size,  specific 
gravity,  particle  shape  and  temperature  of  the  water.  Usually,  the  fall 
velocity  is  calculated  from  the  sediment  particle  diameter  as  though  that 
particle  were  a snhere.  Although  it  seems,  intuitively,  that  the  fall  velo- 
city should  be  corrected  for  the  shaoe  of  the  particle,  presently  there  is 
not  enough  information  to  make  such  a correction.  Therefore  the  fall 
velocity  is  calculated  assuming  the  particle  is  a sphere.  The  following 
table  2.02  shows  fall  velocities  for  quartz  sand  (which  has  a specific  grav- 
ity of  2.65)  as  a function  of  temperature  and  particle  size. 

Fall  velocities  for  other  sizes  or  for  oth^r  specific  gravities  may  be 
calculated  with  the  following  equation: 


Cjj  » drag  coefficient 
d ■ particle  diameter 
g * acceleration  of  gravity 
w * fall  velocity 
y * specific  weight  of  fluid 
■ specific  weight  of  the  sphere 
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Sand  Grain  Settling  Velocity  Versus  Tanperature,  SP.G,  2.65,  Shape  Factor  0.9  1 
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For  Reynolds  numbers,  K,  less  than  0.1,  Stokes  law  gives 

Cq  * 24/11  (2-03) 

where 

II  = w • d/v  (2-04) 

V » kinematic  viscosity  of  fluid 
IR  * Reynolds  number 

For  Reynolds  numbers  greater  than  0.1  there  is  no  simple  expression 
for  the  drag  coefficient  relationship,  and  the  following  curve  is  utilized 
with  equations  2.02  and  f'.'OA  to  calculate  fall  velocity  by  successive  approxi- 
mations. 


REYNOLDS  NUMBER,  R 

Fig.  2.01.  Drag  Coefficient  of  Spheres  Is  « 
Function  of  Reynolds  Number 
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f.  Gradation  Curve 


An  important  property  of  the  sediment  sample  is  its  gradation.  Each 
sample  will  usually  contain  a range  of  grain  sizes,  and  it  is  customary  to 
break  this  range  down  into  classes  to  determine  the  percentage  by  weight  of 
the  total  sample  contained  in  each  class  interval.  The  individual  percentages 
are  accumulated  and  the  graph  showing  grain  size  vs.  the  accumulated  percent 
of  material  that  is  finer  than  that  grain  size  results  in  a gradation  curve. 
This  curve  presents  one  set  of  statistics  for  that  sample  as  shown  in  the 
following  example  (fig.  2.02).  Gradation  can  change  drastically  from  sample 
to  sample.  Furthermore,  four  different  gradations  are  significant:  the 
gradation  of  the  suspended  load,  the  gradation  of  the  bed  load,  the  gradation 
of  the  material  comprising  the  bed  surface  and  the  gradation  of  material 
beneath  the  bed  surface.  These  will  all  be  significantly  different,  with 
respect  to  time,  for  the  same  flow  event.  Gradation  curves  are  often  plotted 
on  log-normal  probability  paper. 

Section  2.03.  Properties  of  Sediment  Deposits 
a . Unit  Weight  of  Sediment  Deposits  and  Consolidation  with  Time 
The  weight  per  unit  volume  of  sediment  deposits  is  denoted  "unit  weight." 
The  fact  that  it  could  be  calculated  from  the  unit  weight  of  water,  the 
specific  gravity  of  the  sediment  particles  and  the  void  ratio  of  the  sediment 
deposit  is  academic  because  the  void  ratio  varies  substantially.  Field  tests 
are  required  to  determine  unit  weights.  The  unit  weight  of  a sand  deposit 
can  be  considered  constant  with  time,  but  clay  and  silt  deposits  consolidate 
by  over  200  percent.  Deposits  which  are  always  submerged  take  longer  to 
consolidate  than  those  which  are  occasionally  exposed  to  the  air. 

The  following  unit  weights  are  often  used  when  field  data  cannot  be 
obtained. 
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. STANDARD  SIEVE  OPENING  IN  INCHES  U S.  STANDARD  SIEVE  NUMBERS  HYDROMETER 

i 3 2 1'/^  I ^ VS  H 3 4 6 8 to  14  16  20  30  40  50  70  100  140  200 


Fig.  2,02.  Gradation  Curve  for  a Bed  Sample 
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Table  2.03.  Unit  Dry  Weight  for  Sediment  Deposits  (1) 

Fully 

Initial  Compacted 

Material  Deposits  Deposits  Consolidation 

Classification  ib/cu  ft  Ib/cu  ft  Coefficient Remarks 


Sand 

93 

93 

0 

Always  submerged 

Silt 

65 

82 

5.7 

Always  submerged 

Clay 

30 

78 

16.0 

Always  submerged 

Silt 

74 

82 

2.7 

Moderate  reservoir 
drawdown  (2) 

Clay 

46 

78 

10.7 

Moderate  reservoir 
, drawdown 

Silt 

79 

82 

1.0 

Considerable  reservoir 
drawdown  (2) 

Clay 

60 

78 

6,0 

Considerable  reservoir 

drawdown 


( 1 ) Reference  1 , page  829 

(2)  Significant  because  drawdown  permits  aeration  and  aeration  results  In 
density  changes. 


b.  Unit  Dry  Weight  of  Mixtures 

The  following  equation  Is  utilized  to  calculate  the  unit  weight  of  a 
mixture  of  sand,  silt  and  cley. 


. ^.nd  ■ «s.nd  *[  ’ ’<>9n)(T-1 -O)]- 

100 


where 

^sand  ■ unit  weight  of  sand 
% 

sand  ■ amount  of  sand  In  total  deposit  In  percent 
y'  ■ unit  weight  of  Initial  silt  deposits 
K*  ■ consolidation  coefficient  for  silt 
y"  • unit  weight  of  Initial  cl«y  deposits 


(2-05) 
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K"  ■ consolidation  coefficient  for  clay 

T ■ life  of  deposit  in  years 
% 

clay  ■ amount  of  clay  in  total  deposit  in  percent 

% 

silt  ■ amount  of  silt  in  total  deposit  in  percent 

■ composite  unit  weight  T years  after  deposition  occurred 

c.  Impact  on  Water  Quality 

As  one  might  suspect,  those  sediments  having  electrical  charges,  the 
clay  material,  play  the  most  active  role  in  water  quality  analysis.  There 
is  evidence  that  heavy  metals,  pesticides,  hydrocarbons,  radioactive  wastes 
and  other  pollutants  are  attracted  to  these  sediments  and  consequently 
zones  of  deposition  may  exhibit  rather  high  concentrations  of  these 
pollutants. 

Section  2.04.  Properties  of  the  Water- Sediment  Mixture 

a.  Concentration 

The  dry  weight  of  a sediment  sample  in  milligrams  divided  by  the 
weight  of  the  water- sediment  sample  in  liters  is  the  common  unit  for 
sediment  concentration.  An  alternate  unit  is  parts  of  sediment  per 
million  parts  of  the  water-sediment  mixture.  The  first  definition  is 
preferred. 

b.  Sediment  Load  Terms 

Sediment  load  refers  to  rates  of  sediment  movement  in  the  stream. 
Quantities  should  be  expressed  as  tons  per  day  and  may  be  calculated 
from  sediment  concentration  and  water  discharge  in  tons.  Six  different 
sediment  load  terms  are  encountered  in  sediment  literature:  suspended 
load , bed  load,  wash  load , bed  material  load,  measured  load  and  un- 
measured load.  These  terms  go  together  in  pairs  to  produce  the  total 
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sediment  load  as  follows: 


Based  on  mode  Based  on  availability  Based  on  method 

of  transport  In  stream  bed  of  quantifying 


'suspended' 

load 

wash  load 

• 

measured  load 

Total  load  » 

+ 

bed  load , 

* 

+ 

tied  material  loaij 

s 

+ 

unmeasured  load 

Suspended  load  refers  to  those  sediment  particles  which  are  trans- 
ported entirely  within  the  body  of  fluid  with  very  little  contact  with  the 
bed. 

bed  load  Is  that  portion  of  the  sediment  load  that  moves  essentially 
in  contact  with  the  bed.  There  is  not  a standard  criteria  for  defining 
bed  load.  However,  It  Is  common  to  consider  that  load  moving  within  two 
grain  diameters  of  the  bed  Is  bed  load.  This  point  Is  not  trivial  when 
comparing  one  analytical  technique  with  another.  The  quantity  of  bed 
load  Is  very  sensitive  to  thickness  of  the  flow  depth  used  In  calculating 
It.  However,  It  usually  amounts  to  a small  fraction  of  the  total  load. 
This  definition  recognizes  that  grain  size  alone  Is  not  sufficient  to 
classify  material  as  suspended  load  or  bed  load.  Hydraulics  of  flow  Is 
also  Involved.  Coarse  material  may  move  as  suspended  load  In  high  energy 
flow  and  bed  load  In  low  energy  flow.  Suspended  plus  bed  load  will  equal 
the  total  sediment  load  moving  at  that  point  In  the  stream. 

Wash  load,  on  the  other  hand,  refers  to  that  portion  of  the  suspended 
load  which  Is  not  found  In  the  bed  of  the  stream.  That  Is,  the  gradation 
of  the  material  In  the  stream  bed  Is  coarser  than  the  gradation  of 
material  In  the  wash  load.  This  Is  an  Important  distinction  for  the 
later  discussion  on  analytical  techniques  since  all  transport  formulas 
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are  based  on  the  presence  of  material  In  the  stream  bed.  The  amount  of 
wash  load,  on  the  other  hand,  depends  solely  on  the  supply  of  sediment 
entering  the  stream  and  cannot  be  calculated  with  transport  formulas. 

(In  many  alluvial  streams  the  wash  load  Is  finer  than  0.0625  mm.) 

Additional  Insight  Into  the  meaning  of  "wash  load"  was  given  by  Dr. 
Hans  Albert  Einstein  In  reference  6 (pages  17-36)  as  follows:  "either 
the  availability  of  (sediment)  material  In  the  watershed  or  the  trans- 
porting ability  of  the  stream  may  limit  the  sediment  load  at  a cross 
section.  In  most  streams  the  finer  part  of  the  load,  I.e.,  the  part 
which  the  flow  can  easily  carry  In  large  quantities.  Is  limited  by  Its 
availability  In  the  watershed.  This  part  of  the  load  Is  designated  as 
wash  load.  The  coarser  part  of  the  load,  I.e.,  the  part  which  Is  more 
difficult  to  move  by  flowing  water.  Is  limited  In  Its  rate  by  the  trans- 
porting ability  of  the  flow  between  the  source  and  the  section.  This 
part  of  the  load  Is  designated  as  bed-material  load." 

Wash  load  Is  often  thought  of  as  the  clay  and  perhaps  silt  material 
moving  In  the  stream.  However,  In  coarse  bed  streams,  such  as  those  with 
gravel  cobble  beds,  wash  load  might  also  Include  sand  sizes  since  these  are 
limited  by  supply.  One  should  be  careful,  however,  to  realize  that  the 
amount  of  sand  on  the  stream  bed  at  low  water  may  be  quite  different  from 
the  amount  at  that  same  location  during  high  water.  Therefore,  It  Is  not 
possible  to  Infer  that  the  gradation  of  bed  material  Is  the  same  during  a 
flood  as  It  Is  during  low  water.  In  all  likelihood,  observations  during 
high  water  would  reveal  the  presence  of  much  more  sand  material  In  the 
stream  bed  than  Is  present  In  samples  of  bed  material  taken  during  low 
water.  The  companion  tensi  to  wash  load  Is  bed  material  load.  As  the 
name  Implies,  grain  sizes  In  the  bed  material  load  are  those  grain  sizes 
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which  are  also  found  In  substantial  quantities  k)  the  stream  bed. 

Most  analytical  techniques  concentrate  on  calculating  the  bed  material 
load.  This  will  oftentimes  be  only  a small  fraction  of  the  total  load 
moving,  and  yet  It  Is  the  fraction  that  contributes  the  most  Information 
In  channel  morphology  studies.  The  Importance  of  bed  material  load  Is  so 
great  that  most  of  the  effort  heretofore  In  the  development  of  analytical 
techniques  has  centered  on  being  able  to  calculate  the  bed  material  load 
moving  In  a stream. 

The  term,  measured  load,  refers  to  all  of  the  sediment  load  that  can 
be  measured  with  sampling  equipment,  which  today  Is  the  suspended  load  to 
within  3 Inches  of  the  bed.  Sampling  equipment  that  measures  the  concen- 
tration of  suspended  material  Is  highly  developed  relative  to  other  sediment 
sampling  equipment.  It  can  operate  over  a range  of  depth  or  at  a point; 
be  hand  held  or  operate  from  a cableway;  and  can  measure  a range  of  particle 
sizes  up  to  2 mm.  This  range  of  operation  provides  very  acceptable  results 
when  one  Is  Interested  In  predicting  total  sediment  yield  from  the  watershed 
because  that  type  of  sampling  would  usually  capture  all  but  about  10  percent 
of  moving  sediment.  However,  In  channel  morphology  studies  the  remaining  10 
percent,  the  unmeasured  load.  Is  the  most  significant  fraction. 

The  final  term,  unmeasured  load,  refers  to  that  portion  of  the  sediment 
load  which  Is  not  measured  and  Is  usually  the  material  either  greater  than 
2 mm  or  moving  within  3 Inches  of  the  bed.  It  Is  usually  determined  by 
calculations  using  the  measured  load  as  a guide.  However,  In  coarse  bed 
streams,  the  size  of  material  moving  In  suspension  will  often  exceed  2 mm, 
and  the  unmeasured  load  takes  on  more  significance.  In  terms  of  total 
sediment  yield,  the  unmeasured  load  may  still  be  a very  small  amount. 
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However,  In  terms  of  channel  aggradation  or  degradation  trends.  It  will 
be  a very  significant  load. 

c.  Sediment  Yield 

Sediment  yield  Is  the  annual  amount  of  sediment  delivered  from  a 
drainage  area,  I.e.,  the  total  annual  sediment  load  passing  a gage  on 
a stream.  It  Is  expressed  In  tons  per  year.  A companion  term  Is  unit 
sediment  yield,  and  this  Is  the  annual  sediment  yield  divided  by  the 
drainage  area  - for  example,  tons/year/square  mile.  It  Is  customary  to 
correlate  the  annual  sediment  yield  In  tons  with  the  annual  water  yield 
In  acre  feet  as  shown  In  the  following  figure.  The  relationship  Is 
usually  satisfactory  and  the  scatter  that  does  exist  will  often  exhibit 
a trend  which  can  be  used  to  further  understand  sediment  production  In 
the  basin.  (For  example,  the  five  high  points  In  fig.  2.03  were  ex- 
cluded from  the  regression  analysis  because  they  occurred  during  five 
years  starting  In  the  dry  1930's  whereas  the  remaining  points  are  data 
from  the  late  forties  through  the  fifties  and  reflect  current  land  use 
In  the  basin.) 

Occasionally  annual  sediment  yield  Is  expressed  In  volume  units, 
but  this  Is  not  reconmended . Transport  equations  require  weight  and 
sediment  measurements  and  are  presented  as  weight  per  unit  time.  To 
transfer  from  volume  to  weight  requires  data  on  unit  weight  of  deposits 
which  not  only  Is  difficult  to  determine  but  also  changes  with  time. 
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Correlations  between  annual  water  yield  and  annual  sediment  yield 
indicate  that  a better  approach  for  transferring  sediment  yield  data  from 
one  basin  to  another  is  to  relate  the  intensity  of  sediment  production, 
tons/acre  to  the  intensity  of  water  production,  cfs/acre,  in  the  model 
basin,  determine  the  water  yield  from  the  basin  of  interest  and  then  trans- 
fer the  sediment  data  through  the  observed  correlation.  Even  when  the 
basin  of  interest  is  ungaged,  the  water  production  can  be  approximated  by 
techniques  for  transferring  water  data  from  gaged  areas.  Usually,  water 
data  is  more  readily  available  than  sediment  data. 

Neither  of  these  techniques  establish  the  coarser  load  fraction  of 
the  sediment  yield  (i.e.,  medium  sand  and  up)  because  these  sizes  move  as 
a function  of  their  availability  in  the  stream  bed  and  the  hydraulics  of 
flow.  Fortunately  the  coarse  material  does  not  contribute  substantially 
to  the  total  sediment  yield  — usually  less  than  5 or  10  percent.  The  im- 
portance of  coarse  material  is  discussed  later  in  terms  of  channel  regime. 

Occasionally,  units  of  volume  are  used  to  express  sediment  yield. 

This  is  particularly  true  when  yield  is  calculated  from  measuring  the 
accumulation  of  deposits  in  a reservoir.  Two  difficulties  arise  from  this 
method:  (1)  the  trap  efficiency  of  the  reservoir  must  be  estimated  and 
(2)  the  unit  weight  of  deposits  must  be  available  so  sediment  volumes  may 
be  converted  to  weights  for  subsequent  use  in  analytical  studies.  Trap 
efficiency  is  often  very  difficult  to  establish  since  one  representative 
value  must  be  determined  for  all  magnitudes  of  flood  events.  This  is 
discussed  in  more  detail  in  Chapter  5.  Unit  weight  of  the  actual  deposits 
should  be  measured.  The  "typical"  values  of  trap  efficiency  presented  in 
this  volume  are  not  sufficiently  accurate  for  use  in  such  calculations. 
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CHAPTER  3.  RIVER  MORPHOLOGY 
Section  3.01.  Introduction 

Rivers  are  dynamic  not  only  In  the  sense  that  water  In  motion  Is  dynamic 
but  also  In  the  sense  that  the  size  and  alignment  of  their  channels  are 
continually  changing.  Work  to  set  the  fluid  body  In  motion  Is  provided  by 
the  potential  energy  gradient  from  the  relief  of  the  watershed.  The  river 
channel  serves  as  a focal  point  for  energy  dissipation  as  flow  from  the 
entire  watershed  Is  collected  there.  In  fixed  bed  hydraulics,  energy  Is 
dissipated  by  friction,  expansion  and  contraction  losses  as  discussed  In 
Volume  6, “Water  Surface  Profiles"  In  natural  streams,  however,  part  of  the 
energy  Is  used  to  transport  the  water,  and  part  Is  used  to  transport  the 
sediment  material  which  moves  In  the  water.  The  energy  equation  advanced 
In  Volume  6 to  describe  the  conservation  of  energy  for  flow  within  a fixed 
boundary  does  not  adequately  describe  flow  In  a movable  boundary  because  It 
does  not  provide  a term  for  the  energy  required  to  transport  the  sediment 
material.  At  this  point  In  time,  no  generally  acceptable  equation  has  been 
advanced  to  describe  the  movement  of  water-sediment  mixture  In  a movable 
boundary. 

Section  3.02.  Lane's  Equation  of  Dynamic  Equilibrium 

The  fact  Is  generally  recognized,  however,  that  nature  maintains  a 
delicate  balance  among  the  water-sediment  mixture  flowing  In  a natural 
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stream,  sediment  material  forming  the  boundary  of  the  stream  channel,  and 
the  hydraulics  of  flow.  In  195b  Lane  summarized  this  balance  with  a 
qualitative  statement  whicn  Included  the  bed  material  load,  Q^,  sediment 
size,  D50,  water  discharge,  Q,  and  energy  gradient,  S,  as  follows: 

• DbO  = 0 • S (3-01) 

He  concluded  that  a channel  is  maintained  In  dynamic  equilibrium  by 
balancing  changes  in  the  sediment  load  and  sediment  size,  with  compensating 
changes  in  tne  water  discharge  and  the  energy  gradient.  Although  many 
empirical  equations  have  been  developed  around  the  variables  In  this 
expression,  its  main  value  Is  to  show  qualitatively  the  Impact  of  changes 
in  tlie  water-sediment  mixture  on  the  behavior  of  the  river  channel  conveyance 
system. 

Section  3.03.  Dynamic  Equilibrium  of  Stream  Bed  Profiles 

If  the  yield  of  bed  material  load  should  Increase  while  the  water  dis- 
charge remains  constant,  according  to  expression  3-01,  either  the  effective 
size  of  particles  In  the  sediment  load  must  decrease  (the  sediment  load 
must  become  finer)  or  the  system  will  be  out  of  balance.  A system  which  Is 
out  of  balance  In  this  manner  will  adjust  Itself  by  increasing  the  energy 
slope  until  the  inflowing  bed  material  load  can  be  transported. 

If,  at  some  future  time,  the  Inflowing  bed  material  load  returns  to 
Its  original  value,  the  aggradation  trend  will  cease  and  degradation  will 
begin.  Starting  at  the  upstream  end,  sediment  material  will  be  removed 
from  tne  stream  bed  and  the  slope  will  return  to  Its  original  value. 

Although  this  is  a very  simple  Illustration,  It,  and  numerous  variations. 
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can  be  observed  In  nature.  Whereas  the  sediment  load  was  assumed  to  increase 
in  the  above  illustration,  a more  natural  situation  is  for  the  water  discharge 
to  remain  a constant  value,  the  effective  sediment  size  to  be  fairly  constant 
and  the  energy  slope  to  vary  from  point  to  point  as  in  a natural  river. 
Expression  3-01  indicates  the  sediment  load  must  change  when  the  energy 
slope  does,  and  this  requires  sediment  reservoirs  from  which  material  can 
be  withdrawn  as  the  energy  gradient  increases  and  into  which  material  can 
be  deposited  when  the  energy  gradient  decreases.  The  stream  bed  provides 
those  reservoirs,  and  the  exchange  of  material  involves,  primarily,  the 
coarser  particles  on  and  near  the  bed  surface. 

It  will  be  convenient  for  subsequent  discussions  to  identify  a location 
where  sediment  material  is  being  deposited  as  a sink  and  a location  from 
which  sediment  material  is  being  removed  as  a source.  Part  of  the  great 
complexity  surrounding  river  behavior  is  the  shifting  of  sources  and  sinks. 
Shifting  can  be  related  to  two  factors: 

(1)  The  locations  of  sinks  and  sources  are  related  to  the  location  of 
hydraulic  controls. 

(2)  If  the  volume  of  bed  material  at  a source  is  exhausted,  entrainment 
of  material  will  shift  to  the  next  available  location  downstream 
and  this  may  be  a former  sink. 

As  discussed  in  Volume  6,  "Water  Surface  Profiles,"  the  energy  gradient 
along  a stream  is  controlled  by  cross  section  shape  and  size  in  key  locations 
called  control  sections.  The  energy  gradient  becomes  steeper  as  flow 
approaches  the  control  location.  Reaching  a maximum  at  the  control,  the 
energy  gradient  decreases,  often  abruptly,  only  to  repeat  the  process  at 
the  next  control  downstream.  Consequently,  expression  3-01  would  identify 

the  control  location  as  the  downstream  limit  of  a source. 
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During  the  passing  of  a large  flood,  more  and  more  of  the  hydraulic 
controls  become  Ineffective.  Thus,  new  source/sink  locations  are  estab- 


lished and  old  ones  eliminated.  Center  bars  may  form  where  none  had  existed 
before  simply  because  an  unusually  large  percentage  of  the  total  flow  was 
forced  Into  the  channel  at  one  of  these  new  control  locations.  Just  upstream 
from  a center  bar,  formed  In  this  fashion,  one  will  find  the  deeply  scoured 
area  that  will  become  a new  sink  when  normal  flow  conditions  return. 

Section  3.04.  The  River  One  Observes 

Flood  events  pass  quickly  and  a river  channel  seldom  reaches  an 
equilibrium  condition  during  a single  event.  This  complicates  all  attempts 
to  correlate  observed  channel  conditions  with  flow  and  sediment  load.  The 
river  one  observes  reflects  three  magnitudes  of  flood  events:  (1)  the  low 
flow  events  which  actually  develop^  a low  water  channel  which  meanders  within 
the  river  channel;  (2)  the  normal  flood  events  which  mold  the  river  channel; 
and  (3)  the  most  recent  "superflood"  event  which  molded  the  floor  of  the 
valley.  It  Is  Important  to  recognize  that  conditions  lirmedlately  following 
a superflood  will  reflect  one  extreme  stage  of  equilibrium  in  the  river,  but 
quite  a different  stage  of  equilibrium  will  evolve  after  a sufficiently  long 
period  of  normal  years.  This  Is  particularly  significant  In  levee  design 
where  the  water  surface  profile  for  the  highest  equilibrium  condition  must 
be  determined. 

It  Is  not  safe  to  associate  the  superflood  event  with  geological  time. 
Certainly  a Standard  Project  Flood  (defined  In  Volume  5,  "Hypothetical  Floods") 
would  be  a superflood  and  the  100  year  flood  would  also  most  likely  qualify. 
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As  used  here»  "superflood"  Identifies  a flood  event  which  Is  sufficiently 
large  so  Its  friction  slope  approximates  that  of  the  valley  rather  than  the 
channel  slope. 

Expanding  on  the  second  factor  presented  above,  the  water  discharge 
does  not  have  to  change  for  the  locations  of  sources  and  sinks  to  shift. 

When  all  available  sediment  material  Is  removed  from  a source  area  the  water 
Is  left  with  an  excess  of  transport  capacity.  Material  previously  deposited 
at  a downstream  location,  perhaps,  will  become  the  first  source  for  satisfying 
that  excess.  This  constant  shifting  of  sources  and  sinks  would  eventually 
smooth  out  the  profile  and  produce  a rather  uniform  movement  of  material  If 
the  water  discharge  became  fairly  constant.  However,  fluctuations  In  the 
water  discharge  hydrograph  prevent  this  from  happening. 


Section  3.05.  Rivers  In  Regime 


a.  Width,  Depth  and  Slope  Equations 

Rivers,  It  Is  said,  have  four  degrees  of  freedom:  width,  depth,  slope, 
and  meander  pattern.  The  following  equations  are  advanced  to  describe  width, 
depth,  and  Indirectly,  the  slope. 


(3-02) 

(3-03) 


(3-04) 
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where 


width 


Cg,  Cjj.  Cy 


D 

Q 

V 


Og,  Og,  Oy 


empirical  coefficients  in  the  width,  depth  and  velocity 
equations,  respectively 

depth 

representative  water  discharge 
average  velocity  of  flow 

empirical  exponents  in  the  width,  depth  and  velocity 
equations,  respectively 


To  satisfy  the  principle  of  continuity  it  follows  that: 

ag  + ag  + ay  - 1 

and  S * ^ 


(3-05) 
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Investigations,  reported  by  Leopold,  et  al.(1953),  at  20  river  cross 
sections  in  the  Great  Plains  and  the  Southwestern  United  States  resulted  in 
the  following  average  values: 

Og  ■ 0.26,  Og  = 0.40,  Oy  * 0.34 

It  should  be  pointed  out  that  these  exponents  are  not  necessarily  trans- 
ferable from  one  stream  to  another,  or,  for  that  matter,  from  one  location 
to  another  on  the  same  stream.  Average  values  for  exponents  are  given  for 
158  gaging  stations  in  the  United  States  in  reference  12. 

Undoubtedly,  the  experienced  student  of  river  geomorphology  can  gain 
considerable  insight  about  the  behavior  of  rivers  by  applying  these  regime 
equations.  However,  they  confront  the  novice  with  some  rather  difficult 
basic  questions,  such  as:  “Where  should  one  select  a cross  section  for 
study?  How  does  one  determine  a representative  discharge?  Why  is  sediment 
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load  omitted  from  the  equations?  Why  conduct  a study  at  stream  gage  1oca> 
tions  since  they  are  not  typical  of  the  reach?”  Nevertheless,  this  regime 
approach  provides  a point  of  view  about  the  behavior  of  alluvial  streams 
which  should  not  be  overlooked. 

b.  Meander  Patterns 

The  fourth  degree  of  freedom,  the  tendency  to  meander.  Is  a natural 
fact;  and,  therefore,  river  engineers  usually  avoid  designs  which  produce 
straight  rivers.  Canals  are  often  straight  or  only  slightly  curved,  but 
they  remain  so  because  of  a high  level  of  maintenance  activity  and  because 
the  water  discharge  Is  relatively  constant. 

The  term  "meandering"  Is  used  to  refer  to  the  S-curve  pattern  so 
typical  of  alluvial  streams  as  well  as  to  the  movement  of  the  river  channel 
by  the  formation  and  destruction  of  bends.  Accepting  meander  as  an  independ 
ent  degree  of  freedom,  the  dynamics  of  the  process  become  an  outstanding 
consideration.  Flood  control  with  levees  requires  that  the  meandering  ten- 
dency  be  controlled.  At  the  present  time,  the  meandering  process  Is  not 
well  understood.  "Channel  loops  form  and  migrate,  bends  lengthen,  points 
extend,  cutoffs  occur  and  the  process  Is  repeated  In  such  a manner  that  the 
width  of  the  belt  over  which  this  takes  place  Is  rather  constant  In  a river 
reach.  . . .According  to  Krumbein  and  Sloss,  It  has  been  estimated  that  the 
average  width  of  the  meander  belt  Is  15  to  20  times  the  width  of  the  stream. 
Levees  are  usually  located  within  this  meander  belt. 

Constraining  the  natural  tendency  of  the  channel  to  meander  impacts  on 
the  behavior  and  sediment  transport  capacity  of  the  river.  This  Is  not 
understood  well  enough  to  present  even  a qualitative  analysis  of  the  problem 

^Reference  2,  p.  1-17. 
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However,  It  seems  that  If  the  channel  Is  free  to  shift  laterally,  aggradation 
trends  due  to  deposition  of  the  coarser  sands  are  offset  by  the  river's 
freedom  to  select  a location  through  smaller  particle  sizes.  The  point  bar 
that  results  becomes  a sink  for  the  coarser  sands.  When  the  freedom  to  make 
such  shifts  Is  taken  away  because  of  bank  stabilization,  the  coarser  sands 
deposit  on  the  stream  bed,  rather  than  point  bars,  with  a resulting  aggrada- 
tion trend.  When  designing  levees,  one  should  consider  the  consequence  of 
having  to  raise  the  levee  at  some  future  time  to  just  maintain  the  capacity 
to  pass  the  design  flood.  The  Implication  Is  that  a well -designed  data 
collection  program  Is  essential. 

Modification  of  the  meander  pattern  by  constructing  cutoffs  Is  a design 
consideration.  However,  equation  3-01  Indicates  the  slope  of  a stream  will 
not  change  unless  either  the  water  discharge,  sediment  load  or  effective 
grain  size  of  material  being  transported  Is  changed.  Therefore,  one  might 
be  skeptical  about  affecting  a permanent  lowering  of  the  water  surface  profile 
by  the  construction  of  cutoffs. 

The  study  of  meander  patterns  Involves  the  analysis  of  long  reaches  of 
rivers.  Permanent  controls  need  to  be  located.  If  available,  to  Identify 
these  reaches.  An  artificial  control  established  at  one  location  can  Impact 
on  the  meander  pattern  of  the  river  channel  at  downstream  locations.  It  Is 
Important,  therefore,  to  design  over  long  reaches  and  Implement  these  designs 
to  Insure  that  construction  at  one  point  does  not  cause  undesirable  conditions 
to  develop  at  some  downstream  point. 
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Section  3.06.  Dowlnant  Discharge 


The  possibility  of  selecting  one  representative  Mater  discharge  that 
can  be  considered  to  have  a dominant  effect  on  river  regime  Is  very  attractive, 
not  from  the  standpoint  of  basing  all  studies  of  river  behavior  on  that 
floM,  but  from  the  point  of  view  that  such  a representative  discharge  offers 
a beginning  point  for  studies. 

Whereas  the  concept  appears  to  be  sound: 

"There  Is  not  general  agreement  on  the  proper  value  or  method 
of  determining  this  discharge.  Carlston  concluded  that  the 
dominant  discharge  which  partially  controls  meander  wavelength 
lies  between  the  mean  annual  discharge  and  mean  of  the  maxi- 
mum monthly  discharges.  He  also  states  there  Is  some  evidence 
that  slope  controls  the  meander  wavelength.  On  the  basis  of 
logical  reasoning  that  the  dominant  channel -forming  discharge 
should  be  based  on  the  time-duration-weighted  bed-materlal- 
sedlment  transport  capacity  of  the  streamflow,  an  analysis  of 
this  type  for  the  Arkansas  River  at  Little  Rock  was  conducted 
by  the  U.S.  Am\y  Engineer  Division,  Southwestern,  Dallas,  Texas. 

This  study  Indicated  that  the  dominant  discharge  covers  a range 

from  about  125,000  (half  bank-full)  to  175,000  cfs,  or  about 

2 

three  times  the  mean  annual  flow  of  48,000  cfs." 

Another  approach  for  determining  dominant  discharge  was  presented  by 
Johnson  In  a paper  entitled  "Current  Dutch  Practice  for  Evaluating  River 
Sediment  Transport  Processes,"  reference  13.  It  Involves  the  development  of 

2 

Reference  2,  p,  11-20. 
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a diagram  depicting  the  total  amount  of  sediment  being  transported  during 
each  Interval  of  stage  experienced  at  a gage  and  the  calculation  of  the  first 
moment  about  the  abscissa  to  locate  the  dominant  water  level.  The  dominant 
discharge  Is  associated  with  this  water  level  by  using  a stage>d1scharge 
rating  curve.  A sample  K-dlagram  Is  shown  In  the  following  figure. 


Fiq.  3.01.  K-Diagrani 


The  abscissa,  K,  Is  calculated  as  follows: 


K 

where  G 
h 
K 
m 

Ah 


(3-07) 

h • Ah 


3 

bed  load  transport,  m /day 
gage  height 

sediment  load  weighted  by  flow  depth  and  class  Interval 

time  Interval  In  days  during  which  stage  was  within  the 
class  Interval,  Ah 

class  Interval  assigned  to  depth  scale 
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The  dimensionless  coefficient,  p,  Is  defined  as  follows: 


(3-08) 


B = width  of  river  channel  at  the  water  surface,  meters 

d = effective  grain  size  diameter  In  the  bed  load  (dgg),  meters 

2 

g » acceleration  of  gravity,  meters/sec 

G ■ bed  load  transport,  units  of  m /second 

p « a constant  related  to  the  least  energy  consumption 
hypothesis  (reference  15) 

a = (Pj  - p)/p  (3-09) 

p = density  of  water 

p^  » density  of  sediment  particles 

The  bed  load  may  be  determined  In  a variety  of  ways.  The  analytical 
procedure  proposed  by  Johnson  Is  the  Meyer-Peter  and  Muller  equation. 
Reference  14  discusses  this  equation,  and  a brief  summary  Is  presented  here. 

The  original  form  of  the  Meyer-Peter  and  Muller  formula.  In  metric 
units,  for  a rectangular  channel  Is 

/vlV3  2/3 

Y • qT  h • s - 0.047  -Yj  • d + 0.25  [J|  (3-10) 

where 

d ■ effective  grain  size  of  bed  material 

2 

g ■ acceleration  of  gravity,  m/sec 

6"  ■ bed  load  transport  In  metric  tons/sec/meter  of  width 
(submerged  weight) 

h - water  depth,  m 


I 

1 
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Q « total,  water  discharge,  liters/sec 

Qp  ■ that  portion  of  the  total  discharge  which  Is  responsible  for 
the  bed  load  transport,  1 1ters/sec 

1/3 

Kg  ■ the  Strickler  roughness  coefficient  for  the  bed  In  m ' /sec 

1/3 

Kg  » the  grain  roughness  In  m ' /sec 
S > slope  (dimensionless) 

3 

Y ■ specific  weight  of  water,  metric  tons  per  m 

y"  ■ submerged  specific  weight  of  sediment  particles  (y-  - y). 
metric  tons/m^ 


The  relationship  between  specific  weight  and  density  Is 
p * Y/g 

Rearranging  terms  In  equation  3-10  leads  to  the  following 


8 • Yj  • d 


3/2 


• VT 


p • 


Q,y  • h • s 


1 3/2 


a.  d 


0.047 


where 


Qk 


Qb 

T 


(3-11) 


(3-12) 


(3-13) 


As  expressed  In  equation  3-12,  the  units  of  bed  load  transport  are 
metric  tons/sec/meter  of  stream  channel  width,  submerged  weight.  Converting 
this  load  to  metric  tons/sec,  submerged  weight  Is  accomplished  by  multiplying 
by  the  stream  channel  width,  B. 


8'y"s 


B • d^^^  V 


g • P 


Q|^  • h • S 


o • d 


- 0.047 
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(3-14) 


In  English  units  bed  load  transport  Is  expressed  as  Ibs/sec,  dry  weight. 
Changing  from  submerged  to  dry  weight  Is  accomplished  by  Including  the  ratio 
y/Yj  In  equation  3-14,  and  redefining  the  units  of  variables  as  follows  . 
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B > water  surface  width  1n  channel,  ft 

d ■ effective  grain  size  (*150) ♦ ft 

g ■ acceleration  of  gravity,  fps 

6 > bed  load  transport,  tons/ day  dry  weight 

h • water  depth,  ft 

Q|^  > dimensionless  coefficient 

S “ slope,  ft/ ft 

3 

Y ■ specific  weight  of  water,  lb/ ft 

3 

Y^  * specific  weight  of  sediment  particles,  1b/ ft 
Yj  • Yj  - Y 

o ■ yJ/y  or  (p  -p)/p  ■ dimensionless  coefficient 


,3/2 


. Vg 


tt 


- 0.047 
« 
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An  even  more  common  set  of  English  units  for  bed  load  transport  Is 
tons/day,  dry  weight 


G - 43.2  • G' 


(3-16) 


In  evaluating  Q|^  It  is  necessary  to  determine  how  much  of  the  total 
water  discharge  Is  used  to  transport  the  bed  load.  This  Is  the  water  moving 
In  the  channel 

Qb  • Kg  • Ag  . Rg^^^. 


and 


J. 

S " "b 


(3-18) 
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Ag  = cross  sectional  area  for  channel  flow,  B • h 
Rg  = hydraulic  radius  of  flow  on  channel  bed,  A/B 
ng  = channel  n- value  includes  grain  and  form  roughness 


The  grain  roughness  may  be  calculated  with  the  Strickler  equation  as 
26 


.1/6 


'90 


(3-19) 


where 


d«.^  = grain  size  for  which  90  percent  of  material  is  finer 


'90 


and 


ng  = 1/Kg  (metric)  = 1.486/Kg  (English) 


Taking  the  first  moment  of  the  K-diagram  about  the  abscissa  involves 


I 


DM 


_ £ K ♦ Ah  • h 


T 


(3-20) 


£ K • Ah 


where 


calculated  with  equation  3-07  for  each  class  interval  of 
depth 


Ah  - the  class  interval  assigned  to  depth 
h > distance  from  the  midpoint  of  Ah  to  the  abscissa 
I the  total  number  of  class  intervals 


f(DW) 


(3-21) 


The  dominant  discharge,  Qg,  is  read  from  a stage-discharge  rating  curve. 
Note  that  this  value  can  change  from  year  to  year.  Its  only  advantage  is  to 
aid  in  simplifying  the  early  phases  of  studies  and  it  should  not  replace  the 
type  of  analysis,  discussed  subsequently,  which  uses  a hydrograph  of  flow. 
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Section  3.07.  The  Impact  of  Tributaries  on  the  Stream  Bed  Profile 


Tributaries  fall  into  two  classes:  (1)  those  transporting  sediment 
that  is  finer  than  the  bed  load  of  the  main  stem,  and  (2)  those  having  bed 
load  material  equal  to  or  coarser  than  that  of  the  main  stem.  The  first 
type  will  assist  the  main  stem  in  transporting  bed  material,  resulting  in 
channel  degradation  and  a decrease  in  slope  downstream  from  the  confluence. 
The  second  type  will  exhibit  the  opposite  trend  with  the  confluence  area 
serving  as  a sink  for  deposition  of  bed  load  until  a flood  on  the  main 
stem  removes  it. 

It  is  Important  to  note  that  the  river  one  observes  reflects  a 
balance  among  the  terms  in  expression  3-01,  especially  downstream  from 
a tributary.  When  man  regulates  the  main  stem  flow  by  a dam,  the  flow 
duration  curve  changes,  which  changes  the  dominant  water  discharge  on 
the  main  stem.  However,  the  tributary  continues  to  bring  in  its  same 
bed  material  load.  Equilibrium  is  upset  and  the  main  stem  often  exhibits 
an  aggradation  trend.  This  trend  results  in  increasing  water  surface 
profile  elevations  and  must  be  taken  into  account  in  the  design  of 
levees  or  the  acquisition  of  lands  to  permit  project  operation. 

Section  3.08.  The  Importance  of  Natural  Levees 

The  typical  ground  profile  across  many  natural  streams  shows  the 
ground  surface  to  slope  away  from  the  channel  in  both  directions. 
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Fia.  .I.n?.  Typical  Cross  Section  Shape 


The  high  ground  adjacent  to  the  channel  1s  referred  to  as  a “natural  levee" 
since  it  is  the  result  of  sediment  depositing  during  flood  events  as  water' 
flows  out  of  the  channel.  The  formation  of  natural  levees  plays  a very  signi- 
ficant role  in  the  overall  development  of  a river  which  is  undergoing  aggrada- 
tion. Unfortunately*  flow  hydraulics  associated  with  the  problem  are  highly 
three-dimensional,  and  analytical  procedures  are  not  available  to  study  the 
problem. 

Qualitatively,  one  impact  of  maintaining  a channel  in  a fixed  location 
appears  to  be  the  continuing  Increase  in  height  of  the  natural  levees.  This 
results  in  increasing  the  channel  capacity  to  the  point  where  the  equilibrium 
between  flow  depth,  width  and  velocity  is  upset.  A consequence  is  the  develop- 
ment of  deep  scour  where  the  flow  is  confined  and  the  subsequent  deposition 
forming  a center  bar  where  the  first  expansion  occurs  downstream. 


3-16 


CHAPTER  4.  COLLECTION  AND  ANALYSIS  OF  DATA 


1 


Section  4.01.  Introduction 

Data  collection  in  alluvial  streams  involves  sampling  the  boundary 
as  well  as  properties  of  the  water-sediment  mixture  that  is  flowing. 

Most  data  collection  programs  are  designed  to  determine  the  annual 
weight  of  suspended  sediment  which  requires  measuring  water  discharge 
and  sediment  concentration.  This  level  of  information  is  essential  for 
calculating  sediment  yield  of  the  watershed.  In  addition  to  the  sus- 
pended sediment  load,  sampling  programs  should  periodically  provide 
data  on  the  amount  of  bed  load  moving.  Particle  size  distribution  in 
the  suspended  load,  in  the  bed  load,  and  in  the  bed  surface  as  well  as 
profiles  along  permanent  sediment  ranges  should  also  be  provided.  If 
degradation  studies  are  contemplated,  sampling  measurements  should 
include  the  particle  size  distribution  of  material  beneath  the  bed 
surface.  All  measurements  of  sediment  loads  should  include  water 
temperature  and  a water  discharge  measurement.  Occasional  measurements 
of  water  surface  slopes  are  desirable.  Each  of  these  types  of  data  is 
discussed  in  the  following  paragraphs. 

Section  4.02.  Suspended  Load  Measurements 

a.  Vertical  and  Lateral  Distribution  of  Sediment 

To  calculate  sediment  discharge  it  is  necessary  to  define  the 
vertical  and  lateral  distribution  of  concentrations  in  the  cross  section 
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and  the  variation  of  the  mean  concentration  with  time.  If  this  distribution 
can  be  defined,  sediment  samples  may  be  collected  routinely  at  a single 
vertical  in  the  cross  section  and  the  mean  concentration  calculated  by 
the  use  of  an  index.  Periodically,  the  distribution  index  should  be 
revised  by  collecting  a comprehensive  set  of  data  at  many  verticals. 
Oftentimes  the  distribution  index  will  vary  with  water  discharge,  there- 
fore it  is  essential  to  analyze  a range  of  water  discharges  before  assigning 
a distribution  index.  The  season  of  the  year  is  another  variable  which 
frequently  correlates  with  the  concentration  distribution  index. 

The  discharge-weighted  mean  concentration  in  a vertical  generally 
is  obtained  from  depth  integrated  samples  collected  with  standard  velocity- 
weighting samplers.  The  horizontal  distribution  of  concentration  may  be 
obtained  from  these  data.  However,  the  determination  of  any  vertical 
distribution  of  concentration  requires  point  sample  data  in  each  vertical. 

The  mean  concentration  in  a cross  section  or  a vertical  is  then  computed 
by  weighting  the  concentration  of  each  individual  point  sample  by  the 
increment  of  discharge  which  it  represents. 

b.  Mean  Concentration  for  Section 

Two  techniques  are  available  for  calculating  the  discharge-weighted 
mean  concentration  in  the  cross  section  from  the  mean  concentrations  of 
the  several  sample  verticals.  If  the  sample  vertical  represents  centroids 
of  equal  discharge  (EDI  method),  the  mean  concentration  is  the  average  of 
the  several  verticals  or  is  the  mean  of  the  composited  samples,  provided 
all  samples  are  of  the  same  volume.  If,  on  the  other  hand,  the  sampled 
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verticals  are  uniformly  spaced  and  the  same  vertical  transient  rate  is 
used  for  all  samples  (ETR  method),  the  mean  concentration  is  the  ratio 
of  the  total  weight  of  sediment  to  the  total  weight  of  water-sediment 
mixture  in  all  the  samples. 

c.  Variation  of  Concentration  with  Time 

Having  obtained  the  discharge  weighted  concentrations  for  the  cross 
sections,  the  next  step  in  computing  sediment  discharge  is  to  translate 
Individual  values  of  concentration  into  a continuous  temporal  concentration 
curve.  This  step  may  be  reasonably  simple  if  values  for  water  discharge 
or  sediment  concentration  do  not  vary  greatly.  However,  at  a new 
station,  lack  of  knowledge  of  these  trends  together  with  the  large 
number  of  variable  conditions  affecting  sediment  erosion  and  transport 
requires  an  intensive  sampling  program,  and  successful  station  operation 
requires  continual  modification  of  the  sampling  program  to  obtain  the 
best  results  commensurate  with  a reasonable  expenditure  of  time  and 
effort.  In  any  case,  concentration  data  should  be  interpreted  and  the 
temporal  graph  prepared  by  personnel  who  have  an  indepth  knowledge  of 
the  sampling  program,  and  of  the  physical  and  cultural  environments 
affecting  the  stream  regimen  and  sediment  sources.  A good  understanding 
of  the  fundamentals  of  sediment  transport  is  essential  (9). 

d.  Calculation  of  Sediment  Load  from  Concentration 
The  sediment  concentration  in  a sample  may  be  determined  as  the 

ratio  of  the  weight  of  the  sediment  to  the  weight  of  the  water-sediment 
mixture.  Because  of  convenience  in  the  laboratory,  it  is  usually 
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expressed  In  parts  per  million  and  defined  as  the  dry  weight  of  sediment 
divided  by  the  weight  of  the  water- sediment  mixture  multiplied  by 
one  million.  It  Is  customary  to  publish  concentrations  as  mg/1,  however, 
so  the  values  determined  In  the  laboratory  must  be  converted  prior  to 
computation  of  the  sediment  discharge.  The  equation  for  sediment 
discharge  Is  as  follows: 

Qs  • Q • Cs  • k (4-01) 

where 

■ the  sediment  discharge  In  tons  per  day 

Q ■ the  water  discharge  In  cubic  feet  per  second 

■ concentration  of  suspended  sediment  In  mg/1 

k ■ coefficient  which  converts  volume  per  second  to 
weight  per  day  as  follows: 

k « .0027  tons  per  day  when  Q Is  expressed 
In  cubic  feet  per  second  or  k « .0864  metric 
tons  per  day  when  Q Is  expressed  1n  cubic 
meters  per  second, 
e.  Conversion  of  Units 

Values  of  the  conversion  factor  C for  converting  parts  per  million 
to  milligrams  per  liter  are  given  In  the  following  table.  When  the 
concentration  exceeds  16,000  ppm  the  following  equation  should  be  used 
to  obtain  equivalent  mg/1. 

Cj  - C • ppm  (4-02) 

where 

ppm  ■ parts  per  million  and 

C ■ the  concentration  of  suspended  sediment  In  mg/1. 
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Table  4.01.  Conversion  Factors*  C,  for  Sediment  i 
Concentration;  Parts  Per  Million  to  Milligrams  Per  Liter 


(The  factors  are  based  on  the  assumption  that  the  density  of  water  is 

1.000  (plus  or  minus  0.005),  the  range  of  temperature  is  0®  - 29®  C,  the 
specific  gravity  of  sediment  is  2.65,  and  the  dissolved-solids  concen- 
tration is  less  than  10,000  ppm.  This  table  supersedes  table  1 in 


Guy  [1969]). 


Concentration 

Range  (ppm)  C 

0 - 15,900  1.00 

16,000  - 46,800  1.02 

46,900  - 76,500  1.04 

76,600  - 105,000  1.06 

106.000  - 133,000  1.08 

134.000  - 159,000  1.10 

160.000  - 185,000  1.12 

186.000  - 210,000  1.14 

211.000  - 233,000  1.16 

234.000  - 256,000  1.18 

257.000  - 279,000  1.20 

280.000  - 300,000  1.22 

301.000  - 321,000  1.24 


Concentration 

Range  (ppm)  C 

322.000  - 341,000  1.26 

342.000  - 361,000  1.28 

362.000  - 380,000  1.30 

381.000  - 399,000  1.32 

400.000  - 416,000  1.34 

417.000  - 434,000  1.36 

435.000  - 451,000  1.38 

452.000  - 467,000  1.40 

468.000  - 483,000  1.42 

484.000  - 498,000  1.44 

499.000  - 514,000  1.46 

515.000  - 528,000  1.48 

529.000  - 542,000  1.50 


This  table  shows  that  concentrations  expressed  in  mg/1  are  the  same  as 
ppm  for  values  up  to  16,000  ppm. 

^United  States  Geological  Survey,  Techniques  of  Water  Resources  Investigations, 
"Computation  of  Fluvial -Sediment  Discharge"  by  George  Poterfleld,  U.  S. 
Government  Printing  Office,  1972,  p.  43. 
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Particle  Size  Analysis  of  the  Suspended  Sediment 
A sufficient  number  of  samples  should  be  analyzed  each  year  to 
determine  representative  values  for  particle  size  distribution  in  the 
water-sediment  mixture.  The  accuracy  of  such  analyses  is  dependent  on, 
among  other  things,  the  quantity  and  physical  characteristics  being 
analyzed.  The  following  table  presents  the  quantity  of  sediment 
required  for  various  methods  of  determining  particle  size  distribution. 

Table  4.02.  Quantity  of  Sediment  for  Particle  Size  Analysis,  in  Grams  (9) 


Method 

Minimum 

Optimum 

Dry  sieve 

SO 

100 

Wet  sieve 

0.05 

1.0 

VA  tube  (1) 

0.05 

1.0  - 7.0 

Pi  pet 

0.8 

3.0  - 5.0 

BW  tube  (2) 

0.5 

0.7  - 1.3 

Table  4*03  shows  the  type  of  data  needed  from  suspended  measurements. 

Of  course,  this  amount  of  detail  is  not  needed  from  every  sample. 

g.  Correlation  Between  Concentration  or  Sediment  Load  and  Water  Discharge 
Attempts  to  correlate  sediment  concentration  with  water  discharge, 
even  when  using  multiple  linear  regression  in  which  other  parameters 
such  as  water  temperature  or  season  of  the  year  were  Included,  have  been 
unsuccessful.  However,  the  mass  of  material  moving,  expressed  in  units 
of  tons  per  d«^.  is  often  correlated  with  water  discharge  to  form  a 
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Table  4.03.  Sunrary  Data  of  Suspended-Sedirient  Particle-Size 
Distribution,  Clearwater  River  at  Spalding,  Idaho 


Sieve 

Size 

(mm) 

Date;  Discharqe  (ft'^/sec);  Concentration  (ir.q/1) 

01-19-74 

28,000 

1441/ 

06-05-74 

63,700 

113 

06-10-74 

46,000 

31 

06-18-74 

117,000 

131 

06-26-74 

67,000 

30 

1.0 

• 

100.00 

100.00 

100.00 

100.00 

.71 

- 

99.65 

93.96 

98.74 

99.07 

.50 

100.00 

- 

- 

- 

- 

.35 

- 

92.01 

77.68 

86.73 

89.44 

.25 

99.6 

79.31 

61.65 

75.63 

83.84 

.10 

- 

67.99 

52.24 

66.79 

30.46 

.12 

9C.3 

62.43 

43.38 

56.34 

76.09 

.09 

- 

56.53 

35.07 

47.33 

70.16 

.06 

97.0 

54.43 

29.65 

42.84 

68.45 

.03 

39 

- 

- 

- 

- 

.016 

70 

- 

- 

- 

- 

.003 

54 

- 

- 

- 

- 

.004 

41 

- 

- 

- 

- 

.002 

36 

- 

- 

- 

- 

—^Conducted  in  Sacramento  Sediment  Laboratory, 

U.S.  Geological  Survey 
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relationship  knovm  as  a sediment  load  rating  curve  (fig.  1 of  Appendix  III). 
Often,  the  scatter  of  data  Is  on  the  order  of  one  log  cycle.  Nevertheless, 
many  sediment  studies  require  a relationship  between  the  average  sediment 
load  and  the  water  discharge,  and  It  Is  possible  to  fit  a curve  through 
such  data  so  the  correct  volume  of  sediment  Is  produced  given  a specified 
water  discharge  hydrograph. 

Sediment  load  rating  curves  are  not  always  suited  for  use  In 
calculating  sediment  yield  because  of  the  wide  scatter  of  data.  On  the 
other  hand,  they  are  convenient  and  are  oftentimes  the  only  data  available. 
Therefore,  they  are  frequently  used  to  estimate  sediment  yield  for 
engineering  projects. 

Section  4.03.  Bed  Load  Measurements 

River  engineers  are  skeptical  about  the  adequacy  of  present  techniques 
for  producing  satisfactory  measurements  of  the  bed  load.  In  terms  of  total 
sediment  yield,  that  portion  of  material  contributed  by  bed  load  Is  usually 
small.  Therefore,  It  has  not  been  of  a great  deal  of  Interest  In  the 
solution  of  practical  problems  In  the  past.  Recent  advances  In  analytical 
capability  afforded  by  simulation  techniques  and  the  electronic  computer 
make  It  possible  to  utilize  bed  load  data  In  doing  channel  aggradation 
and  degradation  studies.  In  these  studies  It  Is  the  bed  material  load 
that  causes  channel  problems.  The  USGS  recently  completed  a three-year 
program  of  Intensive  data  collection  Involving  measuring  the  bed  load  for 
use  In  sediment  distribution  studies  on  the  Clearwater  and  Snake  River 
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arms  of  Lower  Granite  Reservoir  (1972-1974).  The  following  Insight 
was  obtained  from  that  study. 

The  water  sediment  mixture  contained  material  ranging  from  very  fine 
sand  to  256  tnm  cobbles.  (For  the  purpose  of  this  exercise,  bed  load  Is 
defined  as  that  material  moving  within  half  a foot  of  the  channel  bed 
since  that  Is  the  size  of  orifice  In  the  sampler  used  In  the  data 
collection.)  Two  samples  were  collected  at  each  of  10  verticals  going 
across  the  river  and  this  was  repeated  on  the  return  trip  making  a total 
of  40  samples  per  measurement.  Results  showed  that  both  the  temporal 
and  spatial  fluctuation  In  bed  load  Is  greater  than  It  Is  for  suspended 
load,  therefore,  the  concept  of  an  ii.dex  station  and  a reduced  number 
of  samples  does  not  seem  appropriate  for  measuring  bed  load.  Table  4.04 
shows  typical  Information  from  measurements  In  this  program.  The 
Helley-Snilth  sampler  utilized  In  this  program  performed  very  well 
when  compared  with  a continuous  belt  bed  load  sampler  In  another  stream. 
It  Is  presently  undergoing  more  Intensive  calibration  at  the  Federal 
Inter-agency  Sedimentation  Project,  St.  Anthai^' Falls  Hydraulic 
Laboratory,  Minneapolis,  Minnesota. 

Bed  load  measurements  without  a corresponding  particle  size 
analysis  are  not  useful.  Techniques  which  demand  Information  on  the 
bed  load  movement  also  require  a detailed  gradation  of  the  bed  load 
Itself. 
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Table  4.04.  Summary  Data  of  Bedload  Particle-Size  Distribution 
Clearwater  River  near  Spalding,  Idaho 
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— Computed  from  Incomplete  number  of  samples. 
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- Total  sample  weight  » 189  pounds. 


Section  4.04.  Sediment  Yie 1 d 


a.  Calculation  of  Sediment  Yield 

The  best  method  for  determining  annual  sediment  yield  is  to 
accumulate  the  mean  daily  loads.  Usually  the  published  daily  loads 
are  suspended  sediment  only  and  have  to  be  increased  by  5 or  10  percent 
to  include  the  unmeasured  portion  of  sediment  passing  the  gage. 

As  illustrated  in  fig.  2.03, -Sediment  Yield,  the  quantity  of  annual 
sediment  yield  correlates  rather  well  with  water  yield.  Therefore, 
calculating  sediment  yield  from  observed  sediment  loads  will  not 
necessarily  produce  the  long  term  average  or  even  a design  value 
for  sizing  a reservoir.  Therefore,  a sediment  load  curve,  such  as 
presented  in  fig.  1 of  Appendix  III,  should  be  integrated  with  the 
long  term  water  discharge-duration  curve  to  produce  a design  yield 
for  a reservoir  study. 

Land  use  is  a very  important  parameter  in  projecting  sediment 
yield  or  transferring  yield  estimates  from  one  basin  to  another. 

Appendix  III,  entitled  Corps  of  Engineers  Methods  for  Predicting 
Sediment  Yields,  presents  details  of  several  approaches  for  estimating 
sediment  yield. 

Universal  soil  loss  equation.  Some  attempts  have  been  made  to 
calculate  sediment  yield  from  rainfall,  erodability  of  the  soil,  ground 
surface  slope,  lengtii  of  overland  flow,  ground  cover  and  erosion  control 
practices.  Reference  10  is  one  such  example.  This  method  utilizes  the 
following  equation: 

Y^  - R^  • K • S • L • C.  • (4-03) 

e e e f e 
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where 


Yg  * potential  sediment  yield  In  tons/acre 

Rg  = the  rainfall  energy  factor 

Kg  = erodablllty  of  the  soil 

S » slope  of  ground  surface 

L = slope  length 

® cropping  factor 

= erosion  control  factor 
e 

These  terms  appear  to  convey  the  Important  factors  In  the  erosion  process, 
and  the  resulting  quantity  of  sediment  depicts  the  soil  loss  potential  of  a 
specific  small  land  area.  However,  sediment  yield  at  a stream  location 
encompasses  more  than  just  this  soil  loss  potential.  The  eroded  sediment 
material  has  to  be  transported  through  channels  to  an  outflow  point.  Perhaps 
clays  and  silts  In  the  outflow  will  correlate  with  soil  loss  potential 
calculated  In  the  above  manner,  but  the  total  yield  of  sediments  from  a 
drainage  basin  will  depend  primarily  upon  flow  hydraulics  and  the  availability 
of  material  In  and  near  the  stream  channel  or  conveyance  system, 
b.  Transferability  of  Sediment  Yield  Data 

As  one  might  suspect,  the  scatter  In  data  becomes  much  less  when  relating 
annual  sediment  yield  to  annual  water  yield  than  It  Is  at  the  same  gage  for 
dally  values  of  sediment  load  vs.  water  discharge.  Consequently,  annual 
sediment  yield  curves  provide  a very  satisfactory  mechanism  for  transferring 
sediment  data  from  one  basin  to  another  or  from  one  point  to  another  In  the 
same  basin.  This  procedure  Is  not  satisfactory  where  there  Is  a substantial 
difference  In  land  use  between  the  two  gage  sites.  Land  use  In  this  sense 
Is  meant  to  be  the  percentage  of  land  In  the  basin  that  Is  occupied  by 
forest,  pasture,  row  crops  or  by  other  uses.  The  following  figure  Illustrates 
transferring  sediment  yield  data  between  two  points  In  the  same  basin  by  using 


4-12 


tuoi.  40  •uoiniiniai  •Ojommiq  iu*iii|p»8  iiMiuuy 

Fig.  4.01.  Transferring  Sediment  Yield  Data  in  a Basin 


tiie  water  yield  relationship.  Also  a few  observed  data  points  are  shown. 

Ttie  drainage  area  is  17,18b  square  miles  at  Romayor  and  8,146  square  miles 
at  Rosser.  Fig.  4.0?  shov;s  the  unit  water  yield  versus  distance  along  the 
main  stem  river.  Table  4.05  shows  the  actual  calculations  for  distributing 
water  and  sedinient  yields  among  the  various  tributaries  in  this  basin.  In 
general,  only  that  drainage  area  near  streams  or  gullies  contributes 
substantially  to  sediment  yield. 

c . tffect  of  Land  Use  Change  on  Sediment  Yield 

The  impact  of  land  use  on  sediment  yield  can  be  substantial.  For  example. 
When  a strip  mining  operation  enters  a previously  all-forested  watershed, 
sediment  yield  may  increase  one  to  two  orders  of  magnitude.  This  will  impact 
directly  on  reservoirs  in  the  basin  and  perhaps  even  on  the  channel  regime 
i tsel f . 

Section  4.05.  Frequency  of  Measurement 

Fig.  4.03  sliows  tlie  accumulated  weight  of  suspended  sediment  at  Spalding, 
Idaho  as  a function  of  percent  of  time  required.  That  is,  for  water  year 
1974  the  total  period  sampled  amounted  to  about  4 months  (30  percent).  The 
total  suspended  sediment  load  during  that  period  was  1,200,000  tons.  Fifty 
percent  of  this  material  was  transported  during  5 percent  of  the  year 
(18  days).  Horeover,  the  river  transported  more  sediment  during  the  two 
peak  days  of  1974  than  was  transported  during  the  entire  year  in  1973. 

The  Shape  of  the  accumulated  load  curves  in  fig.  4.03  shows  that  essentially 
all  of  the  annual  sediment  load  is  transported  during  flows  which  are 
equaled  or  exceeded  30  percent  of  the  time. 

A sediment  sampling  program  has  to  be  very  responsive  to  flow  con- 
ditions in  the  river.  It  should  include  some  sampling  throughout 
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fig.  4.03.  Accuolulative  Suspended-Sedlaent  Transport  as 
Function  of  Tiae,  Cleaniater  River  «t  Spaulding,  Idaho, 

taken  from  basic  data  prepared  by  U.S.  Geological  Survey 
Boise,  Idaho. 
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the  entire  year.  However,  an  Intensive  effort  should  be  made  during  the 
flood  season.  Initially  a great  variety  of  Information  Is  needed  as 
frequently  as  the  water  discharge  or  flow  characteristics  change. 

Gradually  the  level  of  effort  can  be  decreased  and  a firm  sampling 
schedule  can  be  established. 

Section  4.06.  Sampling  to  Determine  Gradation  of  the  Bed  Surface 

All  sediment  transport  formulas  require  an  effective  grain  size 

and  some,  such  as  the  Einstein  relationship,  require  detailed 
Information  on  the  gradation  of  particle  sizes  for  material  on  the  surface 
of  the  stream  bed.  This  gradation  data,  usually  expressed  as  percent 
finer  vs.  grain  size,  must  be  representative  of  the  entire  stream  bed 
surface  at  that  cross  section.  A further  requirement  Is  that  a 

r 

I sufficient  amount  of  such  data  must  be  available  to  show  how  the 

gradation  of  the  bed  surface  changes  with  changing  water  discharge. 

One  technique  for  developing  a representative  bed  gradation  curve 

* Is  to  collect  a grab  sample  of  material  from  the  stream  bed  at  each 

( 

' vertical  during  a sediment  discharge  measurement  and  analyze  each  sample 

to  determine  Its  gradation  curve.  These  curves  will  usually  vary  a 
great  deal  from  one  vertical  to  the  next,  but  a representative  one  mey 
be  calculated  by  weighting  each  sample  by  the  sediment  load  In  that 
verti cal . 

The  bed  gradation  observed  during  low  flow  Is  not  representative  of 
the  bed  during  a high  flow  event.  Bed  gradation  changes  In  response 
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to  the  sizes  of  sediment  which  are  moving  as  bed  material  load,  and  the 
size  of  the  bed  material  load  varies  with  water  discharge  and  the  supply 
of  material  that  Is  available.  The  supply  of  material  Is  controlled  by 
upstream  sources  and  sinks.  Flow  duration  Is  an  Important  factor  for 
determining  whether  a source  of  material  Is  exhausted  or  not  and, 
consequently,  this  affects  the  bed  gradation.  Therefore,  bed  gradation 
samples  should  be  taken  for  a wide  range  of  water  discharges  and  flow 
durations.  Data  from  each  event  may  be  analyzed  as  discussed  above. 

It  Is  often  Impossible  to  collect  bed  samples  during  a flood 
event.  Two  techniques  are  presented  for  circumventing  this  problem. 

The  first,  and  least  expensive.  Is  to  sample  at  selected  locations  during 
low  flow  periods;  and  the  alternative  Is  to  collect  samples  of  bed  load, 
suspended  load,  water  discharge,  water  velocity,  water  depth,  slope  and 
width,  and  to  calculate  the  gradation  of  bed  material  that  Is  required 
to  transport  the  measured  loads.  These  techniques  are  presented  In  the 
following  paragraphs. 

Recent  experience  with  the  first  technique  on  coarse-bed  streams 
Indicates  that  samples  of  material  near  the  water's  edge  during  periods 
of  low  flow  are  representative  of  the  gradation  In  the  stream  bed  only 
during  periods  of  low  flow.  The  gradation  of  samples  taken  at  about 
midway  between  the  bed  and  the  top  bank  agreed  better  with  calculated 
bed  gradations  In  the  simulation  model  during  flood  flows  when  the 
Inflowing  sediment  load  was  just  transported  without  scour  or  deposition 
of  the  bed. 
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Figure  4.04  entitled.  Gradation  of  Bed  and  Bank  Material,  shows  both 
measured  and  calculated  bed  gradation  on  the  Scioto  River,  Ohio. 

In  this  case,  the  sediment  load  transported  during  flood  discharges 
required  a smaller  effective  diameter  In  the  bed  material  than  that 
transported  during  low  flow  events.  In  both  cases  the  Inflowing 
sediment  load  was  Just  transported.  In  the  simulation  model,  with  neither 
scour  nor  deposition  occurring.  This  comparison  suggests  that  samples 
from  the  stream  bed  should  be  used  when  calculating  sediment  load  for 
low  flow  events,  but  when  calculating  the  load  for  high  flow  events, 
one  should  also  consider  samples  from  the  stream  bank  If  bed  samples 
>kere  not  taken  during  high  discharges. 

The  second  technique  for  developing  bed  gradation  data  Is  to  measure 
the  total  bed  material  load  moving,  to  measure  hydraulic  parameters  and 
to  calculate  the  gradation  of  the  bed  surface  that  Is  required  to 
transport  that  measured  load. 

Section  4.07. 

Sensitivity  of  Transport  Calculations  to  Bed  Surface  Gradation 

It  Is  desirable  to  have  two  samples  at  each  cross  section  when 
utilizing  the  existing  simulation  models  to  calculate  bed  load  transport 
from  bed  gradation  and  hydraulics  of  flow.  However,  In  studies  which 
extend  over  long  distances  such  a data  collection  program  would  require 
too  much  Investment.  Efforts  should  be  concentrated  In  one  short 
reach  with  subsequent  samples  spaced  sufficiently  far  apart  to  provide 
guidance  In  the  calibration  of  the  simulation  model. 
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Fig,  4,04.  Gradation  of  Bed  and  Bank  Material  4.21 


In  general,  bed  samples  should  not  be  combined  Into  one  composite 
value  over  a long  distance.  Such  a technique  filters  out  the  Impact  of 
the  hydraulic  sorting  which  goes  on  continually.  In  some  studies 
hydraulic  sorting  Is  an  Important  consideration. 

Coarse  material  which  Is  present  In  the  bed  samples  should  be 
Included  In  the  analysis.  Rather  than  eliminating  the  large  particles 
to  control  bias.  It  Is  better  to  collect  a large  sample  and  minimize 
the  Impact  of  one  large  particle.  Recent  computations  on  the  Snake 
and  Clearwater  Rivers,  Idaho,  USA,  Indicate  the  representative  bed 
gradation  curve,  defined  as  that  gradation  curve  which  causes  the 
Inflowing  sediment  load  to  be  transported  In  the  simulation  model. 
Appendix  VII,  lies  between  a gradation  curve  developed  as  percent  finer 
by  weight  and  one  which  Is  developed  by  counting  the  number  of  particles 
and  expressing  the  result  as  percent  finer  by  that  number. 

To  a large  extent,  the  variability  between  observed  and  calculated 
bed  material  loads  can  be  attributed  to  bed  surface  particle  sizes  which 
are  not  representative  of  the  prototype  condition.  Sensitivity  of 
calculated  results  to  percentage  of  material  In  the  bed  Is  shown  In 
the  following  table. 

The  first  four  columns  are  typical  of  the  calculations  required  to 
determine  the  sediment  transport  capacity  of  the  bed  material  load  for  a 
given  hydraulic  condition.  The  total  observed  bed  material  load  was 
6000  tons/day.  Column  5 Is  Included  to  Illustrate  the  sensitivity 
of  the  calculated  bed  material  load  to  the  bed  gradation  curve. 
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For  example.  If  the  bed  samples  collected  In  the  field  misrepresent  j 

I 

the  amount  of  Very  Fine  Sand  on  the  stream  bed  by  only  1 percent,  the 
resulting  calculated  bed  material  load  would  have  been  doubled  for 
that  size  fraction.  The  larger  grain  sizes  are  not  affected  this  much. 

This  sensitivity  alone  causes  one  to  question  values  of  bed  material 
load  calculated  by  transport  functions.  However,  particle  sizes  up 
through  medium  sand  move  predominately  as  suspended  load  and  can  be 
sampled  quite  easily.  As  seen  In  table  4.06,  this  Is  the  major  portion 
of  the  total  bed  material  load. 


Section  4.08. 


Beneath  the  Bed  Surface 


Ordinary  data  collection  programs  do  not  address  sampling  of  the 
sediment  material  beneath  the  bed  surface.  Those  situations  Involving 
channel  degradation  and  armoring  require  such  Information,  however. 

Not  only  Is  the  particle  size  analysis  Important,  oftentimes  the  depth 
to  bedrock  Is  an  Important  consideration. 

The  calculations  for  potential  degradation  and  armoring  are  not 
as  sensitive  to  the  amount  of  material  In  the  stream  bed  as  sediment 
transport  calculations  are  because  fine  material  Is  not  Involved. 
However,  It  Is  probably  the  coarsest  2*3  percent  of  material  In  the 
natural  stream  bed  which  accumulates  to  form  an  armor  layer  and  arrest 
degradation  trends. 


I 


Table  4.06.  Sensitivity  of  the  Toffaleti  Sediment 
Transport  Formula  to  Gradation  of  the  hed  Surface 


ir. 


Grain 

Size 

Potential 
Transport 
Capacityd ) 

Gradation 
of  Oed 
Surface 

Transport 

Capacity 

Frror  in 
Oed 

Gradation 

(1) 

tons/ciay 

(2) 

(3) 

tons/day 

(0 

(S) 

VFS 

630000 

.0046 

2030 

6300 

FS 

200000 

.0063 

1053 

2000 

”S 

soooo 

.0370 

1052 

500 

cs 

7700 

.0374 

200 

77 

vcs 

1150 

.0643 

74 

11 

VFG 

310 

.0307 

12 

3 

FG 

130 

.0304 

5 

1 

riG 

60 

.0033 

5 

1 

CG 

26 

.00 

2 

0 

VCG 

10 

.30 

3 

2 

(?) 

.6000 

6224 

8003 

Notes:  (1)  This  is  multiplied  by  the  fraction  of  material  on  the  bed 
surface  to  determine  sediment  load  movinn  in  tons/dav. 


(?)  ^ summation  of  this  column  has  no  meaninq. 
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The  number  and  location  of  samples  required  for  degradation  studies 
are  difficult  and  depend  on  the  stream  and  the  problem  being  analyzed. 

If  the  problem  Involves  local  scour,  the  bed  samples  should  be  taken 
within  the  expected  scour  zone.  If  the  problem  Involves  general  channel 
degradation,  the  samples  should  be  spaced  out  along  the  channel.  The 
significance  of  even  a few  coarse  particles  should  not  be  underestimated. 

Some  locations  show  sediment  gradation  curves  which  do  not  change 
with  depth.  Other  locations  show  a distinct  trend  In  the  gradation  of 
sediment  material  with  respect  to  depth  — particularly  at  several  feet 
below  the  bed  surface.  Where  channel  degradation  Is  anticipated.  It 
Is  essential  to  determine  the  gradation  of  the  material  as  a function 
of  depth  so  that  this  Information  may  be  used  In  calculations  for  the 
armor  layer. 

Section  4.09.  Standardization  of  Sampling  Equipment, 

Techniques  and  Methods  of  Analysis 

The  Committee  on  Sedimentation  of  the  Water  Resources  Council  has 
responsibility  for  standardizing  sediment  sampling  equipment  and 
techniques  In  the  United  States.  This  need  grew  out  of  the  great 
Impact  that  sampling  equipment  and  techniques  have  on  the  scatter  of 
data.  Even  with  standardization.  It  Is  necessary  to  record  the  type 
of  equipment  and  analysis  technique  used  for  each  sample. 
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Section  4.10.  Permanent  Sediment  Ranges 


Because  of  the  uncertainty  In  present  analytical  techniques 
for  handling  sediment  problems  and  to  give  Insight  Into  the  prototype 
behavior  of  sediment  problems,  permanent  sediment  ranges  are  usually 
established  along  reservoirs  and  In  channels  where  aggradation  or 
degradation  trends  are  expected.  These  ranges  are  permanently  marked 
and  resurveyed  periodically  to  determine  the  amount  of  change  In  the 
cross  section  due  to  scour  or  deposition. 

It  Is  essential  that  sediment  ranges  be  close  together  In  places 
where  the  greatest  amount  of  scour  or  deposition  will  occur.  In  the 
case  of  deposition  In  reservoirs,  this  would  be  In  the  upstream  extremity 
of  the  reservoir  where  the  Inflowing  water  establishes  the  backwater 
curve.  In  the  case  of  narrow  reservoirs,  ranges  should  be  located  at 
points  of  expansion.  Ranges  may  be  located  further  apart  through  that 
portion  of  the  reservoir  where  deposits  are  not  expected. 

Imnedlately  downstream  from  the  dam.  It  Is  again  desirable  to  locate 
ranges  close  together  with  the  distance  between  them  Increasing  with 
distance  away  from  the  dam.  Control  locations  should  be  Included  If 
they  are  free  to  degrade.  Sections  should  be  resurveyed  annually 
during  the  early  life  of  a project.  Later  the  time  between  resurveys 
can  be  extended  perhaps  up  to  five  years,  depending  on  how  rapidly 
deposition  or  scour  Is  occurring. 


Sediment  ranges  alone  do  not  provide  useful  data  for  calibrating 


analytical  techniques  such  as  the  simulation  model  In  Appendix  VII. 

The  gradation  and  amount  of  total  Inflowing  sediment  load,  the  Inflowing 
water  discharge,  the  water  temperature,  and  the  operating  policy  of  the 
reservoir  must  also  be  known.  Otherwise,  the  sediment  range  data  will 
only  be  useful  for  monitoring  the  rate  of  deposition  In  the  reservoir. 
This  will  not  be  a linear  function  with  time;  consequently,  projections 
Into  the  future  will  not  be  possible  without  an  analytical  technique. 


4-27 


CHAPTER  5.  RESERVOIR  SEDIMENTATION 


Section  5.01.  Nature  of  Reservoir  Sedimentation  Problems 


The  suninary  statement  In  Chapter  1 addressed  deposition  In  reservoirs 
In  terms  of  the  Conversion  from  kinetic  and  Inertial  energy  to  potential 
energy.  This  chapter  on  reservoir  sedimentation  Identifies  specific 
problem  areas  and  methods  of  solution  conmensurate  with  the  level  of 
detail  required  at  various  stages  between  the  early  planning  and  final 
design  of  a reservoir  project.  Problem  areas  are  associated  with  project 
purposes,  as  well  as  the  physics  of  sediment  movement,  and  the  impact  of 
changes  In  the  bed  elevation  on  the  water  surface  profile.  The  following 
table  correlates  problems  with  project  purposes  in  the  order  of  Increasing 
coirplexity  of  analysis. 

The  first  six  study  objectives,  “Project  Feasibility"  through 
“Aggradation  of  Tributary  Channels"  have  been  successfully  accomplished 
using  analytical  techniques  presented  In  this  volume.  The  "Channelization 
for  Navigation"  requirement  refers  to  constricting  a river  channel  with 
dikes  or  revetments  so  that  navigation  depth  Is  maintained  by  action  of 
the  flowing  water  rather  than  by  dredging.  Good  first-approximations 
to  the  amount  of  constriction  required  have  been  made  using  techniques 
presented  In  this  volume.  However,  final  alignment  of  control 
structures  (dikes  and  revetments)  and  Identification  of  locations 
requiring  over-contraction  should  be  based  on  movable  bed  hydraulic 
model  studies. 
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Shore  erosion  and  beach  Both  a vertical  and  horizontal  distribution  of  sediment  movement 

stability.  including  driving  forces  of  the  3-dimensional  reservoir  current  pattern 

overland  or  local  inflow,  and  action  of  the  ground  water  table. 


Two  and  three  dimensional  mathematical  formulations  of  the  hydro- 
dynamic  equations  are  presently  being  developed  by  several  investigators, 
reference  4.  In  the  future  perhaps  these  will  increase  analytical 
capability  for  density  current  studies.  Presently,  these  studies 
require  hydraulic  model  testing. 

The  study  objectives  Involved  with  siting  facilities  and  locating 
navigation  channels  require  hydraulic  model  analysis.  Shore 
erosion  and  beach  stability  virtually  require  a prototype  scale  of 
testing. 


Section  5.02. 

Data  Requirements  for  Analyzing  Deposition  in  Deep  Reservoirs 


The  amount  of  data  required  to  perform  the  various  types  of  studies 
listed  in  table  5.01  increases  with  each  item  on  the  list.  The  first 
step  is  to  estimate  the  total  amount  of  sediment  that  will  be  available 
for  deposition  during  the  design  life  of  the  project.  This  requires  the 
following  data: 

Table  5.02.  Data  for  Total  Deposition  in  a Reservoir 


Type  of  Data 

Source 

Units 

Design  life  of  reservoir 

Pol  i cy 

Years 

Reservoir  capacity 

Topographic  maps 

m^  or  ac.  ft 

Water  yield  of  watershed 

Streamgage  records,  com- 
putations, or  judgment 

m^/year  or 
ac.  ft/year 

Sediment  yield  of  watershed 

Direct  measurements,  com- 
putations, or  judgment 

m^/year  or 
ac.  ft/year 

CempetUlon  of  sediment  ma- 
tmrial  (n  terms  of  sand,  silt 

ame  c ley 

Direct  measurements 

Percent 

»’•  9^  smdiment 

Direct  measurements  or 

Ib./ft^  or 

• •• 

computations 

tons/m^ 

with  this  level  of  Information  one  can  predict  the  volume  of 
sediment  deposits  using  trap  efficiency  techniques.  This  Is  an  essential 
step  In  the  early  planning  for  a potential  reservoir  project. 

Assuming  flood  control  Is  a project  purpose,  the  next  level  of 
detail  In  reservoir  sedimentation  studies  Is  to  divide  the  total  volume 
of  predicted  deposits  Into  that  volume  settling  In  the  flood  control 
pool  and  that  volume  settling  In  the  remainder  of  the  reservoir. 

Fig.  l.UI  shows,  conceptually,  the  deposition  pattern  to  expect. 

This  level  of  calculation  requires  the  following.  In  addition  to  data 
listed  above: 

Table  5.03.  Data  for  Calculating  Deposition  In  Flood  Control  Pool 


Type  of  Data 

Source 

Units 

Depth  of  flood  control  pool 

- 

ft.  or  m 

Depth  of  reservoir 

- 

ft.  or  m 

Percent  of  time  the  water  level  In  the 
reservoir  Is  at  or  above  the  bottom  of 
the  flood  control  pool 

Calculations 

Percent 

The  next  level  of  detail  In  reservoir  sedimentation  studies  Is  to 
determine  the  location  of  sediment  deposits  In  the  reservoir  so  the  new 
bed  profile  can  be  determined.  This  calculation  addresses  the  land 
acquisitlon/levee  height  and  channelization  problems.  In  addition  to 
data  shown  above,  this  calculation  requires  the  following: 
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Table  5.04.  Data  for  Distributing  Sediment 
Deposits  Along  the  Reservoir 


Type  of  Data 
Cross  sections 

Operating  rule  for  the  reservoir 
elevation 


Source 

Topographic  maps 
Calculations 


Units 
ft.  or  m 


Hydrograph  for  Inflowing  water 
discharge 


Direct  measurements 
or  calculations  cfs  or  cms 


Size  of  sediment  particles  entering  Direct  measurements  nm 
the  reservoir 


Concentration  of  sediment  In  Direct  measurements  mg/1 

Inflowing  water 

Water  temperature  Direct  measurement  "F  or  °C 

or  judgment 


Data  requirements  for  hydraulic  model  studies  are  not  presented  In  this 
volume. 

The  first  three  study  objectives  do  not  require  a great  deal  of 
calculations  and  can  easily  be  accomplished  by  manual  methods.  Dn  the 
other  hand,  determining  the  location  of  sediment  deposits  requires 
voluminous  amounts  of  calculations  and  those  studies  are  best 
conducted  with  automatic  data  processing  equipment. 


Section  S.D3.  Calculation  of  Reservoir  Deposition  Using  Trap  Efficiency 

The  term  "trap  efficiency"  Is  defined  as 
% ■ (fsl  - 1so>«s1 


s-os 


r 


where 

Tg  ■ trap  efficiency  expressed  as  a decimal 
Qsi  ■ inflowing  sediment  load 
Q^o  ” outflowing  sediment  load 

Methods  which  utilize  trap  efficiency  to  calculate  the  volume  of 
sediment  deposits  are  sometimes  referred  to  as  detention-time  methods. 
This  name  stems  from  correlations  which  demonstrate  a relationship 
between  trap  efficiency  and  the  time  required  for  water  to  flow  through 
the  reservoir.  One  such  relationship*  developed  by  Gunnar  Brune,  (3). 
is  shown  in  fig.  5.01. 


MOTC^  Ra»«cduced  from  Gunner  Brwne't 

Artieto,  "Tree  CffiCteneir et  Rtearveire* 
^wblithed  in  Trent  A.G.U.,  «kne  IBGB. 


In  utilizing  Brune's  curve,  flow  through  time  Is  calculated  as 


follows: 

Tj  - S/Y^  (5-02) 

where 

3 

S ■ Capacity  of  reservoir  In  ac.  ft.  or  m 
Tj  ■ flow  through  time  (or  detention  time)  In  years 
« average  annual  water  yield  In  units  of  volume/year 

The  following  example.  Lowsed  Reservoir.!/  Illustrates  the  calculation 
of  total  volume  of  sediment  deposits  using  Brune's  curve. 

Given: 

Design  life  of  the  reservoir  (TIME)  ■ 100  years 

Reservoir  capacity,  top  of  flood 

control  pool  (S)  ■ 3060  acre  feet 

Average  annual  water  yield  (Y)  ■ 36200  acre  feet/year 

Average  annual  sediment  yield  (Q^^)  ” 16000  tons/year 

Composition  of  sediment  material:  clay  ■ 25% 

silt  - 35% 

sand  ■ 40% 

total  ■ 100% 

The  unit  dry  weight  for  sediment  deposits  Is  not  available  for  a 
proposed  reservoir.  Therefore,  the  values  suggested  In  this  volume, 
table  2.03.  Unit  Dry  Weight  for  Sediment  Deposits,  are  utilized  along 
with  equation  2-05  to  calculate  an  average  unit  weight  at  the  end 
of  100  years. 

1/  This  example  was  developed  by  Mr.  Ernest  L.  Pemberton.  United  States 
Bureau  of  Reclamation.  Denver.  Colorado. 
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Table  5.05. 

Calculated  Unit  Dry  Weight  of  Sediment  Deposits 

Sediment 

Material 

^1 

fiCf 

K 

log^o(T) 

Y,+K-log,Q(T) 

% 

TOO 

E£l 

(1) 

(2) 

(3) 

(4) 

(5) 

(6) 

(7) 

Clay 

30 

16 

2 

62 

0.25 

15.5 

Silt 

65 

5.7 

2 

76 

0.35 

26.6 

Sand 

93 

0 

- 

93 

0.40 

37.2 

79.3 

It  Is  convenient  to  convert  sediment  Inflow  from  the  average  annual 
amount  In  tons/year  to  the  total  volume  delivered  during  the  100-year  period 
of  analysis  In  acre  feet. 

V. 


Q tons  . 1 ft^ 
si  year  T^qo 


2000  lb  1 ac.  ft  . 

to?  ?35OT7t3  y«*" 


• ^s1  200000 


• 16000  200000 

"70  * 43660 


Vj,  - 926  ac.  ft. 


(5-03) 


The  next  step  Is  the  calculation  of  trap  efficiency.  Since  this 
coefficient  depends  on  the  capacity  of  the  reservoir  and  the  capacity 
changes  with  time*  It  Is  customary  to  calculate  a value  for  the  beginning 
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and  one  for  the  end  of  the  design  life  of  the  project.  These  calculations 
are  shorn  below,  starting  with  Initial  conditions. 


"d  • 

- 3060/36200 

Tj  • 0.084  years  (5.04) 

Entering  fig.  5.01  with  a detention  tlwe  of  0.084  years,  the  trap 
efficiency  Is  read  directly  as  85  percent.  The  resulting  amount  of 
sediment  deposition,  may  be  calculated  with  the  volume  equivalent  of 
equation  5.01. 


Q - Q ■ T . V 

si  ^so  e d 

Total  Volume  Deposited  ■ % * V^  * 0.85  • 926  ac.  ft. 

■ 787  ac.  ft. 


(5-05) 


The  final  step  Is  to  refine  these  calculations  by  assunlng  the 
reservoir  storage  capacity  Is  decreased  by  787  acre  feet  and  recalculating 
the  trap  efficiency  following  the  same  computation  procedure  as  before. 


The  result  Is  81  percent.  Utilizing  the  beginning  and  ending  trap 
efficiencies  It  follows  that  the  average  trap  efficiency  for  the  100 
year  project  life  Is  83  percent  which  results  In  769  acre  feet  of  sediment 
deposits  In  the  reservoir. 


Grain  size  does  not  enter  directly  Into  these  calculations  for 
trap  efficiency,  although  It  Is  apparent.  Intuitively,  that  the  larger 
grain  sizes  would  deposit  first  and  only  the  smallest  grain  sizes  would 
be  transported  through  the  reservoir.  The  Impact  of  grain  size,  then.  Is 
Inherent  In  the  empirical  data  used  to  develop  the  trap  efficiency 
curve.  On  the  other  hai.d,  approximating  detention  time  by  utilizing 
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total  reservoir  capacity  and  annual  water  yield  does  not  merit  refinement 
of  data  to  Include  particle  size.  Even  so,  the  Brune's  curves  may  be 
used  to  give  a reasonable  estimate  of  the  total  volume  of  sediment 
deposits  to  expect  In  a reservoir. 

In  making  sedimentation  studies  for  Uardanelle  Reservoir  on  the 
Arkansas  River,  the  U.S.  Ani\y  Corps  of  Engineers  District,  Little  Rock, 
Arkansas,  USA,  refined  the  trap  efficiency  calculations  by  establishing 
two  trap  efficiency  curves.  One  of  these  curves  was  representative  of  the 
sand  load  and  the  other  curve  of  the  silt  and  clay  load.  The  reservoir 
was  subdivided  Into  reaches  and  the  storage  In  each  reach  determined  so 
that  for  a given  water  discharge  the  detention  time  could  be  calculated 
for  each  reach  and  the  sediment  deposition  could  be  calculated  separately 
for  sand,  and  silt  and  clay.  The  Little  Rock  District  trap  efficiency 
curve  Is  presented  in  Appendix  IV. 

Section  S.04.  Development  of  Trap  Efficiency  from  Observed  Data 

The  procedure  which  Little  Rock  District  used  In  the  development 
of  their  ^r«p.>if fid ency  curve  Is  contained  In  Appendix  IV,  entitled 
"Procedure  for  Development  of  LRD  Trap  Efficiency  Curve." 

Section  5. 05.  Calculation  of  Deposition  In  the  Flood  Control  Pool 

Fig.  5.02  shows  an  empirical  relationship  which  can  be  used  to 
< 


determine  how  much  of  the  total  sediment  deposition  can  be  expected  In 
the  flood  control  pool  (fig.  1.01).  The  "Flood  Pool  Index"  Is  a measure 


ri5.  5.0?.  Deposition  In  tOe  F1oo<<  Control  Pool 


of  Uie  amount  of  time  the  reservoir  remains  In  the  flood  pool  and  the 
relative  depth  of  that  pool.  The  index  Is  defined  by  equation  5.06. 

Necessary  paraweters  for  using  this  curve  are  listed  In  table  5.03, 
"uata  for  Calculating  Deposition  In  Floou  Control  Pool."  The  calculations 
whlai  require  successive  approximations,  utilize  the  reservoir  depth  versus 
storage  capacity  curve  as  follows. 


The  required  flood  control  storage  at  LoMsed  Is  650  ac.  ft.  The 
Initial  reservoir  capacity  at  the  top  of  flood  control  pool  Is  3060. 
Assuming  no  sediment  deposits  for  the  first  appT*ox1mat1ont  the  reservoir 
volume  at  the  bottom  of  the  flood  control  pool  Mould  be  2410  ac.  ft. 

The  corresponding  pool  elevation  Is  97  feet  msl.  Only  3 feet  of  depth 
Is  required  In  the  flood  control  pool  (I.e.,  100.0  > 97.0}  and  the 
depth  beloM  flood  control  pool  Is  33.7  feet.  The  elevation-duration 
curve,  fig.  5.03a,  shOMS  that  the  reservoir  Is  In  the  flood  control  pool 
only  4.1  percent  of  the  time,  based  upon  average  yearly  conditions. 

A flood  pool  Index  Is  calculated  as 


Flood  Pool  Index 


Depth  of  Flood  Pool 
Depth  BeloM  f^lood  Pool 


Percent  of  Time  In  Flood  Pool 


- (3/33.7)  • 4.1 

■ 0.36  (5,06) 


Entering  fig.  5.02  Nith  this  Index,  the  resulting  amount  of  deposition 
In  the  flood  control  pool  Is  read  as  3 percent.  The  resulting  100 
year  volume  Is  3 percent  of  787  ac.  ft.  which,  rounded  off.  Is  24  ac.  ft. 

This  additional  storage  for  sediment  does  not  change  the  flood  pool 
depth  on  the  second  approximation  and,  therefore.  It  Is  adopted  for 
planning  purposes. 

As  a rough  check,  the  available  conservation  storage  at  the  end  of 
100  years  may  he  determined  as  follows: 

n 

i 


f 


-1  - 


I 
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ELEVATION  (IN  FEET) 


0.  STAGE  DURATION  CURVE 


b.  AREA  AND  CAPACITY  CURVES 


Fig.  5.03.  Duration  and  Capacity  Curves  for  Lowsed  Reservoir 
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storage  at  top  of  conservation  pool 

Sediment  deposition  In  conservation 

pool 

Storage  at  bottom  of  conservation  pool 

Section  5.06.  Distribution  of  Sediment  Deposits  Along  the  Reservoir 

The  planning  or  design  of  a reservoir  requires  an  analysis  to  determine 
how  sediment  deposits  w111  be  distributed  In  the  reservoir.  This  has 
been  the  most  difficult  aspect  of  reservoir  sedimentation  to  deal  with 
because  of  the  complex  interaction  between  hydraulics  of  flow,  reservoir 
operating  policy.  Inflowing  sediment  load,  and  changes  In  the  reservoir 
bed  elevation.  The  traditional  approach  to  analyzing  distribution  of 
deposits  has  relied  on  empirical  methods,  all  of  which  required  a great 
deal  of  simplification  from  the  actual  physical  problem. 

Conceptually,  deposition  starts  In  the  main  channel.  As  flow 
enters  a reservoir,  the  main  channel  fills  at  the  upstream  end  until  the 
elevation  Is  at  or  above  the  former  overbank  elevations  on  either  side. 

Flow  then  shifts  laterally  to  one  side  or  the  other,  but  present  theory 
does  not  predict  the  exact  location.  During  periods  of  high  water 
elevation,  deposition  will  move  upstream.  As  the  reservoir  Is  drawn 
down,  a channel  Is  cut  Into  the  delta  deposits  and  subsequent  deposition 
moves  material  farther  Into  ttie  reservoir.  The  lateral  location  of  the 
channel  may  shift  from  year  to  year,  but  the  hydraulic  characteristics 


« 3060  - (650  +24) 

■ 2386  ac.  ft. 

- 787-24 
* 763  ac.  ft. 

■ 1623  ac.  ft. 
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of  it  will  be  similar  to  those  of  the  natural  channel  existing  prior 
to  Impounding  the  reservoir.  Vegetation  will  cover  the  exposed  delta 
deposits  and  thus  attract  additional  deposition  until  the  delta  takes  on 
characteristics  of  a flood  plain. 

The  diameter  (size)  of  sediment  particles  commonly  transported  by 
streams  ranges  over  five  log  cycles.  As  flow  enters  a reservoir,  coarser 
material  deposits  first  and  the  finest  material  Is  carried  well  Into  the 
reservoir,  even  to  the  point  of  passing  the  dam.  In  order  to  calculate 
the  volume  of  material  which  will  deposit  as  a function  of  distance, 
grain  size  must  be  Included  as  well  as  the  magnitude  of  the  water  discharge 
and  the  operating  policy  of  the  reservoir. 

Reservoir  shape  Is  an  Important  factor  In  calculating  the  deposition 
profile.  For  example,  flow  entering  a wide  reservoir  spreads  out,  thus 
reducing  transport  capacity,  but  the  path  of  expanding  flow  does  not 
necessarily  follow  the  reservoir  boundaries.  It  becomes  a 2-d1mens1onal 
problem  to  calculate  the  flow  distribution  across  the  reservoir  In  order 
to  approximate  transport  capacity  and  therefore  the  resulting  deposition 
pattern.  On  the  other  hand,  flow  entering  a reservoir  that  Is  narrow  has 
a more  uniform  distribution  across  the  section  resulting  In  hydraulic 
conditions  that  are  better  approximated  by  one  dimensional  hydraulic 
theory. 

Flood  waves  attenuate  upon  entering  a reservoir.  Therefore,  their 
sediment  transport  capacity  decreases  from  two  considerations:  (1)  a 
decrease  In  velocity  due  to  the  Increase  In  flow  area  and  (2)  a decrease 
In  velocity  due  to  a decrease  In  water  discharge  (attenuation)  resulting 
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from  reservoir  storage.  As  reservoir  storage  Is  depleted  by  the  sedi- 
ment deposits  In  the  delta,  the  Impact  of  attenuation  on  transport  capacity 
diminishes.  The  resulting  configuration,  therefore,  is  assumed  to 
depend  upon  the  first  consideration  whereas  the  time  for  delta 
development  Is  Influenced  somewhat  by  the  second  consideration. 

A very  useful  technique  for  estimating  the  deposition  of  suspended 
sediment  In  deep  reservoirs  Is  presented  In  Appendix  VI.  This 
computer  program  does  not  actually  change  bed  elevations;  consequently 
the  user  should  run  two  conditions  before  accepting  the  results: 

(1)  Initial  Led  elevation  condition,  and  (2)  calculated  bed  condition 
for  end  of  reservoir  life,  determined  from  (1). 

The  Engineering  Technical  Letter,  "Distribution  of  Reservoir  Sediment 
Deposits,"  ETL  1110-2-64,  dated  7 July  1969,  Department  of  the  Anny, 

Office  of  the  Chief  of  Engineers,  Washington,  DC  20314,  by  Mr.  Brice 
L.  Hobbs  Is  Included  In  Appendix  V to  give  Insight  Into  reservoir  sedimenta- 
tion problems.  In  this  ETL  Mr.  Hobbs  points  out  several  weaknesses  In 
"present  empirical  methods"  (1969)  for  distributing  sediment  deposits. 

These  Involve  "total  sediment  loads,"  "average  trap  efficiencies,"  and 
"gross  volumes  of  sediment  trapped."  The  weaknesses  have  been  overcome 
to  a large  extent  by  using  analytical  methods  that  can  be  applied  with 
the  aid  of  the  electronic  cornputer. 

The  final  approach  presented  In  this  chapter  Is  based  on  digital 
simulation,  using  the  computer  program  "Scour  and  Deposition  In  Rivers 
and  Reservoirs,"  described  In  Appendix  VII.  This  approach  eliminates 
the  requirement  of  working  with  many  o^  the  averages  and  approximations 
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of  the  reservoir  behavior  and  Increases  the  engineer's  capability  to 
approach  sediment  distribution  problems  from  the  standpoint  of 
simulation,  that  Is,  from  the  standpoint  of  Including  a great  deal  of 
detail  about  reservoir  operation,  the  Inflowing  sediment  load,  grain 
size  of  the  material,  the  Inflowing  water  discharge,  hydrograph 
and  reservoir  shape,  and  changes  In  bed  elevation  with  respect  to  time. 

However,  It  Is  a one  dimensional  steady  flow  technique  and  does  not 
attempt  to  locate  the  position  of  the  channel  through  delta  deposits. 

It  does  recognize,  however,  that  ultimately  the  channel  will  exist  and 
will  be  similar  (In  geometric  size)  to  the  natural  channel  that  existed 
prior  to  Impounding  the  reservoir.  The  results  obtained  enable  prediction 
of  an  ultimate  deposition  profile  through  the  delta  deposits. 

Thus,  the  Impact  of  deposition  on  the  water  surface  profile  and  on  the 
upstream  limits  of  deposition  can  be  determined. 

Probably  the  most  difficult  problems  Involving  the  distribution 
of  sediment  deposits  are  associated  with  the  flow  through  shallow 
reservoirs.  Therefore,  the  example  chosen  to  Illustrate  the  aforementioned 
simulation  techniques  Is  a reservoir  project  which  Is  sufficiently 
shallow  so  that  during  some  Inflow  periods  the  deposits  are  actually 
entrained  from  the  bed  and  moved  further  down  Into  the  reservoir. 

The  project  Is  shown  on  fig*  5.04.  This  Is  a multipurpose 
project  Including  hydroelectric  power  and  navigation  among  other 
project  purposes.  The  problem  Is  complicated  by  the  fact  that  levees 
are  required  for  flood  protection  In  the  vicinity  of  the  upstream 


limit  of  the  reservoir  project.  Therefore,  it  is  necessary  to  estimate 
the  impact  of  future  sediment  deposits  on  levee  heights. 

Basically,  the  data  required  are  cross  sections  and  Manning 
n-values  for  use  in  backwater  computations;  the  inflowing  sediment  load 
in  tons  per  day  as  a function  of  the  water  discharge,  the  gradation  of 
material  in  the  stream  bed;  the  operating  rules  at  the  dam,  the  water 
inflow  hydrograph  to  the  reservoir  and  the  water  temperature.  Briefly, 
the  sequence  of  events  in  the  simulation  study  are  as  follows:  The 
hydrograph  of  flows  is  converted  to  a histogram  for  which  each  discharge 
has  an  associated  duration  in  days.  The  first  discharge  event  is  entered 
and  a water  surface  profile  calculated  to  determine  the  velocity,  slope, 
depth  and  width  at  each  cross  section.  The  inflowing  sediment  load 
associated  with  that  water  discharge  is  entered  at  the  upstream  end 
and  routed  through  the  reservoir  using  the  hydraulic  parameters 
mentioned  above.  Deposition  or  scour  is  then  calculated  for  each  reach 
based  on  sediment  inflow  and  outflow  for  that  reach.  Finally,  the 
amount  of  change  in  the  bed  elevation  is  calculated  for  each  cross 
section  and  cross  section  coordinates  in  the  movable  bed  are  changed. 

This  completes  the  first  pass  through  the  calculations.  A sample  of  the 
data  Involved  is  illustrated  in  the  following  figures  and  tables. 

Fig.  &.05  shows  a typical  cross  section  with  the  limits  of 
movable  bed  assigned.  The  user  specifies  the  portion  of  the  cross  section 
over  which  deposition  will  occur. 

Hydraulic  roughness  is  assumed  to  be  the  same  throughout  the  life 
of  the  project.  This  is  a rather  significant  assumption  since  hydraulic 
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roughness  depends  on  both  the  grain  size  and  bed  form.  The  assumption 
Is  that,  as  the  new  equilibrium  bed  Is  approached,  both  grain  size  and 
bed  form  will  return  toward  those  existing  In  the  natural  state.  One 
should  evaluate  the  sensitivity  of  such  an  assumption  by  varying  the 
n-value.  Better  hydraulic  roughness  theories,  than  are  presently 
available  are  needed. 

Table  5.06  shows  the  type  and  level  of  detail  available  for 
describing  the  Inflowing  sediment  load.  Both  bed  load  and  suspended 
load  are  considered.  This  table  corresponds  to  a water  discharge  of 
100,000  cfs.  In  the  actual  study  similar  Information  was  developed  for 
water  discharges  varying  from  5000  to  200,000  cfs.  Neither  of  the 
rivers  entering  the  reservoir  are  high  contributors  of  sediment. 

The  Clearwater,  for  example,  conveys  only  about  100  acre  feet  a year 
and  the  Snake  produces  only  about  320  acre  feet  per  year.  Nevertheless, 
over  the  project  life,  both  streams  have  the  potential  of  depositing 
a sufficient  depth  of  material  to  make  a substantial  Impact  on  the  water 
surface  profile  elevation  for  the  levee  design  flood.  Another  point  about 
the  Inflowing  sediment  load  relates  to  the  relative  amounts  of  bed  load 
and  suspended  load.  Since  the  bed  load  comprises  only  about  2 percent 
of  the  total  Inflowing  load,  one  Is  tempted  to  rely  only  on  suspended 
load  measurements  and  neglect  the  bed  load  altogether.  However,  for 
a water  discharge  of  100,000  cfs,  calculations  show  that  90  percent 
of  the  gravel  (bed  load)  deposits  upstream  from  mile  2.6  and  the  sand  does 
not  deposit  until  the  flow  reaches  about  mile  2,  fig.  5.06.  This 


Table  5.06.  Distribution  of  Sediment  Load  by  Grain  Size  Class 

Clearwater  River 

Water  Discharge,  cfs  .,35.000  Total  Bed  Load,  tons/day 

Total  Suspended  Load,  tons/da 
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sorting  of  grain  sizes  Is  particularly  significant  over  the  life  of 
the  project  where  the  volume  of  bed  load  will  be  sufficient  to  make 
a substantial  Impact  on  the  water  surface  profile  and  must  therefore 
be  Identified  separately. 

The  gradation  of  the  stream  bed  Is  not  particularly  significant 
In  deposition  studies  since  that  gradation  Is  being  continually 
changed  by  the  material  depositing.  In  fact*  gradation  of  the  stream 
bed  and  the  detailed  Information  on  Inflowing  sediment  load,  shown  In 
table  5.06,  are  redundant  Information.  One  Is  not  free  to  specify 
both  the  gradation  In  the  stream  bed  and  the  detailed  Information  about 
the  sediment  Inflow  since  the  two  data  sources  are  linked  together  by 
sediment  transport  theory.  Therefore,  rely  on  measured  bed  material 
loads  and  not  on  values  calculated  from  bed  gradation  curves. 

In  most  studies,  bed  load  measurements  are  not  available  and 
the  bed  load  Is  calculated  by  measuring  gradation  of  the  stream  bed. 

This  technique,  when  attempted  for  this  study,  produced  an  amount  of  bed 
material  load  equal  to  twice  the  measured  amount.  However,  because  of 
the  great  variability  In  grain  sizes  deposited  on  the  bed,  this  factor 
could  have  Just  as  easily  been  10  or  20.  This  problem  Is  acute  in 
armored  streams  such  as  the  Clearwater  and  the  Snake,  because  their 
capacity  to  transport  sediment  far  exceeds  the  supply  of  sediment  avail- 
able. The  computations  become  quite  sensitive  to  the  percentage  of 
material  In  the  bed,  and  It  Is  extremely  difficult  to  sample  a bed 
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for  Mhich  the  grain  size  ranges  from  very  fine  sand  to  cobbles  and 
larger  and  have  any  assurance  that  the  sample  Is  representative  of  the 
entire  bed.  Therefore,  In  Initial  problem  areas  such  as  this  Involving 
levees  In  an  urban  area,  measure  all  the  data  presented  In  table  5.06. 

The  Inflowing  water  discharge  hydrograph  was  analyzed  In 
periods  varying  from  one  day  during  the  peak  flood  season  up  to  monthly 
values  during  the  period  of  low  flow.  A sample  period  Is  shown  In 
fig*  5.07.  Water  temperature  was  obtained  from  records  at  a nearby 
stream  gaging  station  and  the  operating  policy  of  the  reservoir  was 
specified.  In  this  case,  to  be  a constant  pool  elevation  738  for  the 
entire  50  years  of  water  discharge  hydrograph. 

The  results  of  the  analysis  are  shown  on  fig.  5.08,  Deposition 
In  Lower  Granite  Reservoir.  In  addition  to  calculating  the  average 
bed  profile  and  the  resulting  Impact  on  the  levee  design  flood  profile, 
the  simulation  technique  calculated  the  trap  efficiency  over  the 
50  year  life  of  the  project  to  be  30  percent.  Also,  the  amount  of 
maintenance  dredging  required  to  support  navigation  up  to  mile  2 was 
calculated.  Other  simulation  runs  were  made  to  test  the  Impact  of 
the  maintenance  dredging  on  both  the  sediment  deposits  and  the  water 
surface  profile. 

Section  5,07.  Impact  of  Natural  Aggradation  Trends 

Using  a simulation  technique,  such  as  the  Scour  and  Deposition 
program  In  Appendix  VII,  to  calculate  the  distribution  of  sediment  deposits 
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Fig,  5.0y.  Deposition  in  Lower  firanite  Reservoir 
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DISTANCE  IN  MILES 


ir  *.  reservoir  requires  a three  step  procedure.  The  first  step  Is  to 
calibrate  the  model  to  the  particular  problem  area.  The  second  step  Is 
to  perform  a base  test.  This  Is  a test  which  demonstrates  that  the 
simulation  technique  and  the  calibrated  digital  model  can  duplicate 
the  behavior  of  the  study  reach  during  a historical  period.  The  third 
step  Is  to  run  the  desired  studies. 

In  performing  the  base  test  It  often  becomes  apparent  that  the 
study  area  Is  a natural  point  of  aggradation  even  without  a reservoir 
project.  This  does  not  Imply  that  the  reservoir,  when  Impounded,  will 
trap  an  extraordinarily  large  amount  of  sediment.  Two  points  should 
be  made.  In  the  first  place,  sediment  yield  from  the  watershed  has 
the  greatest  Impact  on  the  amount  and  distribution  of  deposits  In  a 
reservoir.  In  the  second  place,  aggradation  trends  In  natural  rivers 
should  not  be  associated  with  geologic  time.  Reaches  will  tend  to 
aggrade  for  a time  and  often  the  trend  will  reverse  and  they  will 
degrade.  Therefore,  the  analysis  of  sediment  deposits  for  a potential 
reservoir  site  which  demonstrates  natural  trends  of  aggradation  may  be 
made  using  the  simulation  techniques  presented  In  the  previous  section 
and  In  the  same  manner.  However,  extra  care  should  be  taken  to  establish 
why  such  a trend  exists.  Perhaps  It  reflects  a recent  Increase  In 
sediment  yield  from  the  basin  In  which  case  historic  yield  curves  should 
be  modified  accordingly. 
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Section  5.08.  Impact  of  Upstream  Reservoirs  In  the  Basin 


1 


It  Is  Important  to  Identify  and  locate  any  reservoirs  In  a basin 
where  a sediment  study  Is  to  be  made.  The  projects  upstream  from  the 
point  of  analysis  potentially  modify  both  the  sediment  yield  and  the 
water  discharge  duration  curve.  Both  are  very  significant. 

The  date  of  Impoundment  Is  essential  Information.  With  this*  one 
may  coordinate  observed  Inflowing  sediment  loads  with  whatever  conditions 
existed  In  the  basin  during  the  periods  selected  for  calibration  and 
verification. 

Useful  Information  on  the  density  of  sediment  deposits  and  the 
gradation  of  sediment  deposits  are  often  available  from  other  reservoirs 
In  the  basin.  Finally*  other  reservoirs  offer  a check  on  sediment  yield. 

Section  5.09.  Erosion  Control  Measures 

Erosion  control  measures  refer  to  land  use  and  farming  practices 
In  the  basin.  One  Is  always  a bit  uncertain  about  how  much  credit 
to  give  erosion  control  measures  when  projecting  sediment  yield  Into 
the  future.  It  Is  very  likely  that  no  credit  should  be  given  for  the 
sand  sizes  and  larger  material.  On  the  other  hand*  the  silt  and  the 
clay  sizes  of  material  m^y  undergo  reduction  by  erosion  control  measures. 
Over  the  long  run*  It  will  be  the  management  of  those  measures  that 
Insures  they  will  continue  to  be  effective. 
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Fanning  practices  conducive  to  erosion  control  were  Implemented 
many  years  ago.  Construction  and  Industrial  practices,  on  the  other 
hand,  are  just  now  beginning  to  Implenent  erosion  control  measures. 
Possible  techniques  are  In  the  form  of  mulch,  of  seeding  disturbed  land 
surfaces,  of  stabilizing  the  soil  with  chemicals  and  of  utilizing 
terraces  and  sediment  traps  to  prevent  material  from  entering  the 
natural  stream  system.  Any  one  construction  project  Is  going  to  be 
short  lived  compared  to  the  life  of  a water  resource  project.  However, 
many  areas  undergo  continual  construction  activity  and  one  cannot 
Ignore  the  Impact  this  has  on  sediment  washing  Into  the  natural  stream 
system.  Potential  Increases  of  50  to  100  times  natural  erosion  rates 
can  be  associated  with  construction  activities. 

Section  5.10.  Land  Use  Changes  In  the  Basin 

Land  use  can  Impact  directly  on  sediment  yield  from  the  basin. 

For  example,  land  being  converted  from  pasture  or  forested  area 
to  row  crop  area  can  generate  10  to  20  times  more  sediment  runoff. 

In  most  sediment  studies  It  Is  customary  to  project  existing  land 
uses  through  the  life  of  the  project.  That  Is,  unless  some  explicit 
alternative  future  Is  being  considered,  the  sediment  yield  from  the 
watershed  Is  considered  to  remain  constant  over  the  life  of  the  project. 
A recent  example  where  land  use  Impacted  substantially  on  sediment  yield 
occurred  when  forested  land  was  converted  to  strip  mining.  Because  of 


the  close  proximity  to  the  main  channel,  the  impact  of  this  land  use 
change  was  manifest  immediately  by  an  increase  in  the  sediment  being 
deposited  in  a local  reservoir. 
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Section  5.11.  Reservoir  Sedimentation  Ranges  for  Post  Project  Surveys 

Cross  sections  should  be  located  as  for  water  surface  profile 
calculations.  These  define  points  of  contraction  and  expansion  and 
points  where  bed  change  Is  expected  or  Is  existing.  Many,  but  not 
necessarily  all,  of  the  same  points  would  be  valuable  for  use  In 
monitoring  the  rate  of  sediment  deposition  In  the  reservoir.  Permanent 
ranges,  called  sediment  ranges,  are  usually  established  for  this  purpose. 
It  Is  essential  that  these  ranges  be  periodically  surveyed  at  the  same 
location.  Therefore,  permanent  range  markers  should  be  Installed 
to  insure  the  range  can  always  be  located.  It  is  always  feasible  to 
discontinue  ranges  that  prove  to  be  non-descriptive,  but  beware  of 
the  tendency  to  establish  too  few.  Simulation  studies  provide  useful 
guidance  In  establishing  sediment  ranges. 

The  frequency  of  which  reservoir  ranges  are  resurveyed  depends 
on  how  much  sediment  deposition  Is  occurring.  In  the  Initial  operation 
of  a project,  annual  resurveys  are  essential.  Later  on,  as  the 
reservoir  delta  develops,  the  change  In  delta  elevation  from  one  year 
to  the  next  Is  not  substantial  and  five  or  perhaps  even  more  years  may 
lapse  between  resurveys. 
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The  volume  of  sediment  deposits  is  readily  determined  by  analyzing 
successive  surveys  of  the  sediment  ranges.  However,  it  is  essential  to 
know  the  inflowing  water  discharge,  temperature,  the  inflowing  sediment 
load,  gradation  of  material  in  the  inflowing  sediment  load  and  in  the 
deposits,  and  the  operating  levels  of  the  reservoir  in  order  to  utilize 
data  collected  from  sediment  range  surveys  to  verify  or  improve  analytical 
techniques.  This  is  an  important  step  in  a project  and  one  which  is 
usually  omitted  because  study  funds  are  not  available  once  the  project 
goes  into  operation.  Those  responsible  for  the  operation  and  maintenance 
j of  projects  are  urged  to  establish  a procedure  by  which  the  testing  of 

\ i 

I analytical  techniques  can  be  accomplished  as  data  becomes  available  I 

j 

during  the  life  of  the  project.  ] 

Section  5.12.  Channel  Behavior  Downstream  from  the  Dam  | 

i 

I 

Degradation  lowers  the  stream  bed  and,  consequently,  the  tail water  I 

rating  curve  at  the  dam.  This  is  important  from  the  standpoint  of  the 
dam  foundation  in  the  case  of  structures  built  on  piling,  and  also  with 
regard  to  impact  on  stilling  basin  performance.  An  important  aspect  of 
degradation  studies  is  to  locate  the  extent  of  degradation  downstream 
from  the  dam.  Rock  outcrops  or  other  hard  bottom  controls  may  reduce  the 
vertical  amount  of  degradation,  but  not  the  distance  over  which  degradation 
problems  exist.  This  distance  is  controlled  by  the  modified  flow  duration 
curve  and  availability  of  bed  material  to  replenish  material  removed  by 
the  reservoir. 
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The  computer  program,  "Scour  and  Deposition  In  Rivers  and 
Reservoirs,"  develops  the  armor  layer  by  the  selective  removal  of  the 
finer  grain  sizes  from  the  bed.  Coarser  sizes  are  removed  at  a much 
slower  rate  according  to  transport  equations,  and  as  a result  the 
coarser  sizes  collect  on  the  stream  bed  as  degradation  continues.  The 
coarse  particles  shield  the  finer  particles  underneath  and  reduce  the 
rate  of  degradation.  The  program  may  be  used  to  calculate  both  the 
distribution  of  deposits  In  the  reservoir  and  degradation  of  the 
channel  downstream  from  the  dam  In  one  continuous  analysis.  This  Is 
particularly  useful  when  some  sediment  passes  the  dam. 

Degradation  studies  are  more  complicated  than  the  calculation  of 
deposition  In  reservoirs  because  degradation  studies  encompass  both 
degradation  and  aggradation  as  flow  passes  from  one  reach  to  another 
In  the  natural  river  downstream  from  the  dam.  Nevertheless,  the  scour 
and  deposition  program  provides  a tool  which  Is  proving  to  be  very 
useful  for  analyzing  degradation  problems  as  Illustrated  In  fig. 

5.09. 

The  degradation  anticipated  downstream  from  Ft.  Randall  Dam  was 
arrested  by  armoring  of  the  bed  surface.  The  size  of  sediment  In  the 
armor  layer  composed  less  than  five  percent  of  the  coarsest  material 
sampled  In  the  natural  stream  bed.  Once  established,  flow  regulation 
prevented  the  armor  from  being  disturbed. 

A rock  outcrop  or  another  reservoir  Is  often  suggested  as  a location 
for  establishing  the  downstream  end  of  a degradation  study  area. 


5-33 


TAILWATER  ELEVATION, IN  FEET  MSL 


DEPOSITION  IN  RIVERS  AND  RESERVOIRS? 

note:  SASIC  INFORMATION  FROM  U.  S.  DEPARTMENT  OF 
ASRICULTURE,  PROCEEOINSS  OF  THE  FEDERAL 
INTER-ASENCY  SEDIMENTATION  CONFERENCE, 
MISCELLANEOUS  PUBLICATION  NO.  STO , tSSS. 


Fig.  5.09. 


Degradation  Downstream  From  Ft.  Randall  Dam 
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Certainly  a reservoir  is  appropriate.  However,  a rock  outcrop  only 
limits  the  magnitude  of  degradation  upstream  from  it  and  not  the 
downstream  migration  of  the  degradation  process.  Two  points  are 
worthy  of  consideration  relative  to  degradation  studies. 

The  first  is  concerned  with  magnitude.  A rock  control  does  not 
stabilize  the  upstream  bed  at  the  elevation  of  the  control.  Bed 
material  will  be  scoured  upstream  from  the  control  until  flow  energy 
is  in  balance  with  resistance  of  the  remaining  bed  material. 

The  second  point  concerns  migration  of  channel  degradation  past 
the  rock  control.  If  other  structures,  such  as  bridge  piers,  water 
intakes,  or  diffusers,  or  if  tributaries  are  present  along  the 
stream,  the  downstream  extent  of  the  degradation  study  area  should 
include  these.  A sufficiently  long  reach  should  be  established  to 
include  the  transition  from  degradation  to  aggradation  commonly 
observed  downstream  from  existing  projects. 

It  is  possible  that  no  degradation  will  occur.  In  fact,  a 
reduction  in  channel  capacity  has  been  observed  in  numerous  cases. 
This  implies  aggradation  is  occurring,  and  probably  to  narrow  the 
channel  in  response  to  the  modified  flow  duration  curve.  In  other 
cases,  the  modified  peak  discharges  lose  transport  capacity  at  major 
tributary  junctions  causing  aggradation. 


CHAPTER  6.  AGGRADATION  AND  DEGRADATION  IN  FREE  FLOWING  STREAMS 


Section  6.01.  Introduction 


The  terns  aggradation  and  degradation  generally  refer  to  trends 
In  behavior  of  the  stream  bed  profile.  An  aggrading  stream  Is  one 
on  which  the  bed  profile  Is  tending  to  become  steeper*  whereas  a 
degrading  stream  is  one  for  which  the  bed  profile  is  tending  to 
become  flatter. 

Streams  aggrade  or  degrade  In  discrete  reaches.  The  process 
alternates  between  aggradation  and  degradation  as  the  flow  passes  from 
one  reach  to  the  next.  The  fact  that  a reach  of  a river  scours  during 
the  passing  of  a flood  event  does  not  make  that  a degrading  reach.  A 
degrading  reach  is  one  which  projects  a net  lowering  of  the  bed 
elevations  with  respect  to  time  or,  perhaps,  a net  widening  of  the  river 
channel  with  respect  to  time.  The  rate  at  which  aggradation  proceeds  depends 
on  sediment  yield,  water  yield,  and  grain  size  of  sediment  particles. 

These  discrete  reaches  are  not  fixed  In  space.  They  tend  to  shift 
back  and  forth  with  respect  to  time,  with  respect  to  the  magnitude  of 
water  discharge  and  with  respect  to  land  use  changes.  Land  use  changes 
are  Important  because  they  generally  Impact  on  sediment  production 
from  the  basin.  Other  changes  which  are  equally  Important  are  the 
construction  of  reservoir  projects  In  the  basin  and  modification  of 
the  water  discharge  duration  curve  (by  Interbasin  transfer  of  flow  or 
reservoirs). 
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Interbasin  transfer  of  water  Is  usually  associated  with  large  scale 
projects.  However,  this  same  concept  Is  applicable,  and  the  consequences 
more  obvious,  on  small  scale  projects  such  as  the  construction  of 
subdivisions,  or  other  urban  developments.  Natural  runoff  patterns  are 
altered,  flow  Is  collected  Into  drains  and  other  conveyances  and 
degradation  problems  occur  when  these  flows  return  to  minor,  natural 
conveyance  channels. 

In  terms  of  Chapters  1 and  3 of  this  volume,  the  flow  duration  curve 
has  been  altered  from  that  naturally  forming  the  conveyance  channel. 

Section  6.02.  Degrading  Reaches 

The  obvious  explanation  as  to  why  a reach  degrades  Is  that  the 
Inflowing  sediment  load  Is  too  small  for  the  transport  potential  of  the 
stream  through  this  reach.  The  reach  Is  providing  a sediment  supply 
for  downstream  points.  Intensive  meander  Is  generally  not  associated 
with  a degrading  reach  but  rather  the  tendency  Is  for  the  stream  to 
become  Incised  and  perhaps  to  develop  cutoffs. 

The  magnitude  of  degradation  Is  arrested  by  several  mechanisms.  A 
rock  outcrop  or  sill  across  the  stream  will  tend  to  stabilize  the  bed. 

The  development  of  an  armor  layer  has  the  same  effect.  Oftentimes, 
the  hydraulic  conditions  downstream  from  the  degrading  reach  will  change 
due  to  a downstream  aggradation  condition  and  this  Influence  extends 
upstream.  The  classical  examples  of  degrading  reaches  occur  downstream 
from  dams. 
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The  stability  of  sediment  material  Is  often  associated  with 
tractive  forces  In  accordance  with  the  following  figure.  Tractive 
force  Is  defined  as 

T » yds  (6-01) 

where 

D > water  depth 

S ■ energy  slope 

Y “ unit  weight  of  water 

It  should  be  pointed  out  that  degradation  Is  a function  not  only  of  the 
stability  of  material  on  the  stream  bed,  but  also  the  amount  of  the 
Inflowing  sediment  load.  A stream  bed  can  be  In  equilibrium  because 
material  is  being  deposited  as  rapidly  as  It  Is  being  removed.  The 
tractive  force  concept  does  not  Identify  this  type  of  equilibrium 
condition. 

Section  6.03.  Aggrading  Reaches 

The  problem  causing  aggradation  Is  just  the  opposite  of  that  causing 
degradation.  Reaches  aggrade  because  the  Inflowing  sediment  load  Is 
greater  than  can  be  transported  out  of  the  reach  by  the  hydraulics 
of  flow.  A tendency  to  meander  Is  oftentimes  associated  with  an 
aggrading  reach,  especially  In  streams  having  fine  sediment  with  a non- 
coheslve  bed  and  banks.  Flow  passing  through  an  aggrading  reach 
causes  hydraulic  sorting  of  grain  sizes.  Therefore,  there  Is  always 
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a deficiency  of  the  coarser  sizes  as  flow  leaves  the  downstream  end 
of  the  reach.  The  aggrading  reach  provides  a temporary  stopping  place 
for  sediment  material  moving  downstream. 

Aggradation  is  arrested  by  reducing  the  inflowing  sediment  load, 
by  sorting  to  produce  smaller  particle  sizes  in  the  water-sediment 
mixture  and  by  shifting. in  the  downstream  control  due  to  a degrading 
condition  at  a downstream  point. 

Section  6.04.  Shifts  In  the  Stage-Discharge  Rating  Curve 

Trends  in  the  stage-discharge  curve  are  good  indicators  of  reach 
behavior.  In  alluvial  streams,  these  curves  typically  exhibit  a wide 
scatter  of  data  which  should  not  be  construed  as  sample  error.  The 
range  in  scatter,  between  high  and  low  data,  is  more  appropriately 
attributed  to  aggradation,  degradation,  changing  bed  forms  in  the 
stream  channel , or  sediment  transport  rates  produced  by  water  temperature 
variances.  The  growth  of  natural  levees  will  increase  the  discharge 
at  which  overbanks  become  effective.  These  factors  should  be 
considered  in  utilizing  the  stage  discharge  rating  curves  in  the 
design  of  engineering  projects,  such  as  levees.  It  is  not  appropriate 
to  develop  an  average  line  of  best  fit  through  the  data,  but  rather 
an  upper  envelope  curve  is  more  appropriate  for  use  in  calculating 
the  height  of  levees.  Navigation  depth  and  rip  rap  design,  on  the 
other  hand,  should  utilize  a lower  envelope  curve  through  the  scatter 
of  points.  If  the  analysis  of  curves  developed  from  successive  years 
of  stage -discharge  data  demonstrates  a trend,  then  subsequent  design 
studies  must  project  the  impact  of  this  trend  on  project  design. 
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Periodic  measuring  of  profiles  across  the  river  will  Identify 
trends  In  lowering  of  the  bed  or  enlargement  of  the  channel.  A meander 
does  not  represent  an  enlargement  of  the  channel  since  both  banks  are 
shifting  In  the  same  direction.  An  enlargement  of  the  channel  would 
be  represented  when  the  two  banks  shift  farther  apart. 

Section  6.05.  Analysis  of  Aggradation  and  Degradation 

The  method  for  analyzing  aggradation  and  degradation  presented  In 
Appendix  VII  of  this  volume  Is  a movable  bed,  digital  model.  It  can  be 
compared  with  a fixed  bed  digital  model  as  shown  In  figs.  6.02  and 
6.03.  Relative  to  fig.  6.02,  If  the  topography  of  the  river  valley 
Is  known  and  the  hydraulic  roughness  Is  known,  the  engineer  may 
specify  a water  discharge  and  starting  elevation,  and  he  may  then 
calculate  the  water  surface  profile.  It  Is  a direct  calculation; 
there  Is  no  feedback  such  as  that  shown  In  fig.  6.03  by  the  double 
arrows.  In  fig.  6.03.  sediment  load  has  been  added,  and  additional 
Information  on  hydraulic  roughness  results.  There  Is  a great  deal  of 
uncertainty  about  how  to  predict  the  change  In  bed  form  and 
resulting  Impact  on  n-values,  but  It  Is  at  this  point  that  the  computer 
program  presented  In  Appendix  VII  breaks  Into  the  loop  for  movable  bed 
calculations.  Basically,  In  utilizing  that  program,  the  assumption  Is 
made  that  n-values  are  either  constant  or  they  vary  In  the  vertical  with 
water  discharge.  The  basic  assumption,  then.  Is  that  bed  form  can  be 
related  to  water  discharge. 
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Fig.  6.02.  fixed  bed  MODEL 
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In  proceeding  with  the  analysis  of  aggradation  or  degradation.  It  Is 
essential  to  follow  the  same  guidelines  given  In  Chapter  5. 

1.  Calibrate  to  an  observed  field  condition 

2.  Verify  that  the  calibrated  digital  model  will  reproduce 
the  behavior  of  the  river  during  a period  of  time,  and 

3.  Analyze  the  degradation  problem  at  hand. 

Analysis  of  mobile  boundai*y  problems  In  free  flowing  streams  adds 

a great  deal  of  complexity  to  sedimentation  studies.  Input  data  Is 

essentially  the  same  as  required  for  reservoir  deposition  studies; 

however,  having  the  capability  for  aggrading  and/or  degrading  complicates 

the  calibration  and  base  testing  phase  of  studies.  In  summary,  the 

required  Input  data  are  as  follows: 

Initial  channel  geometry  Including  a detailed  description  of 
natural  levees  along  the  channel 

Water  discharge  hydrograph 

Inflowing  sediment  load  curve 

Gradation  of  material  In  the  stream  bed 

Water  temperature 

Sediment  transport  formula 

a.  Calibration  of  the  Simulation  Model 

It  Is  essential  to  calibrate  the  water  surface  profile  calculations 
Just  as  for  a fixed  bed  study.  The  techniques  are  the  same.  A range 
of  discharges  varying  from  a low.  In  channel,  discharge  to  the  peak 
flood  of  Interest  should  be  analyzed  to  calibrate  the  n-values  as  the 
first  step. 
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In  fixed  bed  studies*  one  does  not  have  to  concern  himself  with 
the  distribution  of  flow  between  the  channel  and  the  overbanks  as  long 
as  the  distribution  Is  reasonable.  In  movable  bed  studies.  It  Is  a 
portion  of  the  flow  which  travels  In  the  channel  that  transports  the 
bed  material.  For  this  reason.  It  Is  essential  that  the  correct 
distribution  of  flow  be  known.  Two  parameters  that  may  be  evaluated 
are  (1)  cross  section  shape  Itself,  especially  the  elevation  at  which 
water  spills  onto  the  flood  plains  and  (2)  the  distribution  of  n-values 
between  channel  and  overbank.  This  n- value  Information  Is  seldom  ever 
known;  therefore.  It  becomes  a matter  of  judgment  to  select  the  proper 
n- value  for  channel  and  for  overbanks. 

The  third  major  Item  In  the  calibration  Is  to  determine  and  calibrate 
the  Inflowing  sediment  load.  The  amount  and  also  the  distribution  with 
respect  to  grain  size  must  be  determined.  The  final  type  of  data 
requiring  calibration  Is  the  gradation  of  material  in  the  stream  bed. 

The  magnitude  of  water  discharge  selected  for  use  during 
calibration  Is  very  Important.  If  the  flow  Is  too  low,  the  simulation 
model  simulates  a low  water  bed  which  often  will  have  aggradation  and 
degradation  trends  that  do  not  appear  to  match  the  behavior  measured  In 
the  prototype  cross  sections.  If  the  water  discharge  Is  too  large, 
aggradation  and  degradation  trends  may  shift  from  those  observed  In 
the  low  water  bed,  but  again  they  may  persist  In  a manner  that  does 
not  follow  the  overall  behavior  of  the  prototype.  In  order  to  simulate 
the  behavior  of  the  prototype  over  a period  of  time,  one  might 
consider  using  the  dominant  discharge  as  presented  In  Section  3.06. 
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It  will  be  necessary  to  evaluate  the  water  surface  profile  computations 
for  discharges  both  higher  and  lower  than  this  doarinant  discharge. 

However,  primary  emphasis  may  be  placed  on  the  dominant  discharge  to 
calibrate  the  Inflowing  sediment  load,  distribution  by  grain  size  In 
the  Inflowing  load  and  the  representative  gradation  of  the  bed  surface 
Itself. 

The  dominant  discharge  transports  considerably  more  sediment  material 
than  the  average  dally  volume,  and  therefore  It  Is  necessary  to  ratio  the 
time  scale  to  simulate  the  correct  volume  of  sediment  material  moving 
during  a year.  This  Is  quite  straightforward  once  the  total  sediment 
yield  and  the  sediment  load  associated  with  the  dominant  water  discharge 
are  known.  It  should  be  pointed  out  that  sediment  volume  calculations 
are  related  to  the  total  sediment  load  carried  by  the  dominant  discharge 
and  not  just  the  bed  load  when  calculating  a time  scale. 

By  definition  of  dominant  discharge  one  expects  the  calculated 
behavior  of  the  river  channel  to  follow  the  same  trends  as  observed  In 
the  prototype.  If  this  Is  not  the  case,  one  must  re-evaluate  analytical 
model  calibration  for  n-values,  the  amount  of  material  In  the  Inflowing 
sediment  load,  the  distribution  by  grain  size  of  material  In  the 
Inflowing  sediment  load,  and  the  gradation  of  material  In  the  stream 
bed. 

These  calculations  assiane  the  flow  Is  distributed  properly  between 
the  channel  and  overbanks.  The  dominant  discharge  Itself  Is  an  Inbanks 
flow.  It  will  be  a rather  high  flow,  approximating  that  of  bank  full, 
oftentimes.  Specifying  cross  section  geometry  such  that  natural  levees 
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are  simulated  Is  critical  In  this  calibration  study.  If  at  one  point 
or  another  the  water  discharge  does  flow  into  the  overbanks • the 
correct  distribution  of  n-values  between  the  channel  and  the  overbanks 
Is  essential.  A representative  geometry  and  the  proper  distribution  of 
n-values  must  be  developed  before  problems  with  sediment  transport  may 
be  resolved. 

The  following  performance  criteria  Is  suggested  for  use  In  determining 
when  a model  Is  properly  calibrated.  For  a discharge  equal  to  the 
dominant  discharge,  trends  In  stream  bed  elevation  should  approximate 
the  behavior  of  the  prototype.  That  Is.  degradation  trends  In  the  analytical 
model  should  correspond  In  magnitude  and  location  to  those  In  the 
prototype  and,  likewise,  aggradation  trends  In  the  model  should  correspond  In 
magnitude  and  location  to  those  In  the  prototype.  The  sediment 
load  moving  from  point  to  point  In  the  model  may  vary  a great  deal. 

However,  at  points  corresponding  to  gage  locations,  the  sediment 
load  passing  should  match  that  observed  In  the  prototype.  Once  the 
model  Is  calibrated  for  the  dominant  discharge.  It  Is  necessary  to 
test  performance  at  a high  discharge  and  at  a low  discharge.  One  can 
only  check  n-values  and  distribution  of  flow  since,  for  any  single 
dlschar^  other  than  the  dominant  discharge,  the  model  will  not 
necessarily  respond  to  the  trends  observed  In  the  prototype.  Calibra- 
tion Is  completed  with  this  phase  of  testing. 

b.  Base  Test 

The  next  step  In  problem  analysis  Is  to  verify  that  the  calibrated 
model  will  reproduce  behavioral  trends  observed  In  the  prototype.  This 
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Implies  knowing  a starting  condition,  calculating  over  a period  of  time, 
and  comparing  the  results  of  the  calculations  with  an  ending  condition 
In  the  prototype.  Often,  there  are  not  enough  data  to  Identify  a 
starting  and  an  ending  condition.  In  these  cases,  complete  model 
verification  Is  not  possible.  There  are  oftentimes  trends  that  one 
can  observe.  For  example,  the  magnitude  of  any  scour  or  deposition 
should  approximate  that  observed  In  the  prototype  and  the  calculated 
gradation  of  the  bed  surface  shoulti  match  that  observed  In  the  prototype. 

The  process  of  making  small  adjustments  In  the  model  data  In  order 
to  achieve  verification  Is  often  called  fine  tuning.  It  Is  In  this 
phase  of  problem  analysis  that  the  engineer's  understanding  of  river 
morphology  Is  most  essential.  Any  of  the  variables  In  the  Input  data 
are  susceptible  to  manipulation  to  some  degree  or  another.  However, 
adopted  values  must  be  reasonable  and  within  the  range  of  normal 
engineering  experience.  After  the  engineer  has  verified  that  the  model 
can  reproduce  observed  events,  the  study  Is  ready  to  proceed  to  the  third 
and  final  phase. 

c.  Production  Run 

It  Is  most  Important  that  results  from  problem  analysis  be  compared 
with  the  base  test  results  established  In  the  verification  runs.  This 
comparison  will  usually  give  a better  response  than  absolute  magnitudes 
of  values  calculated  In  the  production  run.  The  value  of  simulation 
Is  that  any  of  the  Input  data  can  be  modified  and  the  Impact 
measured  directly  on  either  stream  bed  aggradation  and  degradation 
or  the  water  surface  profile  elevations. 
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Section  6,06,  Sediment  Movement  In  Urban  Conveyance  Systems 

Typically,  urban  conveyance  systems  for  storm  water  runoff  consist 
of  storm  drains  In  pipes  Interconnected  with  open  channels.  On  the 
fringes  of  urban  areas  It  Is  not  unusual  for  flow  to  be  collected  into  > 

a well  designed  conveyance  system  In  a subdivision  only  to  pass  out  Into 
a natural  channel  of  limited  capacity  which  becomes  taxed  with  having 

i 

to  transport  considerably  more  water  than  before  the  subdivision  was 
constructed.  Erosion  of  the  natural  channel  Is  accelerated.  These 
small  conveyance  channels  obey  the  same  regime  concepts  as  those  discussed 
earlier  for  larger  rivers. 

Deposition  In  the  urban  conveyance  system  Itself  Is  a comnon 
occurrence.  Typical  criteria  for  designing  conveyance  systems  specify 
that  velocities  and/or  bed  shear  stresses  will  be  sufficiently  high 
to  prevent  material  from  depositing  In  the  section.  Tractive  force  Is 
a typical  design  criteria.  As  discussed  In  paragraph  6.02,  tractive 
force  alone  Is  not  sufficient  to  determine  If  sediment  will  deposit. 

Tractive  force  Is  a parameter  for  determining  transport  capacity;  whereas 

deposition  depends  on  both  the  transport  capacity  and  the  sediment  I 

1 

yield  from  the  basin. 
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Because  of  the  large  amount  of  construction  activity  in  urban  areas, 
there  is  always  an  abundance  of  sediment  material  available  for  washoff. 
When  the  rate  of  washoff  exceeds  the  transport  capacity  of  the  conveyance 
system,  deposition  results,  and  the  conveyance  system  takes  on  the 
characteristics  of  an  alluvial  stream  in  which  the  bed  slope  is  in 
balance  with  sediment  material  in  motion.  Appropriate  construction 
practices,  including  mulching,  seeding,  terracing,  chemical  applications 
and  detention  storage  ponds,  will  help  to  reduce  the  deposition  problem. 

The  analysis  of  sediment  problems  in  urban  drainage  conveyance 
systems  follows  the  same  procedures  discussed  earlier  for  aggradation 
and  degradation  of  free-flowing  streams.  Again,  the  two  primary 
components  of  the  problem  are  the  determination  of  sediment  yield 
and  the  calculation  of  sediment  transport  in  the  conveyance  system. 

The  sediment  yield  calculations  are  more  complex  than  those  for  a 
major  river  because  the  conveyance  system  is  much  closer  to  the  source 
of  erosion  and  therefore  the  spatial  distribution  of  sediment  becomes 
an  important  consideration.  Transport,  on  the  other  hand,  follows 
the  same  physical  laws  as  that  in  large  rivers.  Transport  theory 
cannot  be  applied  directly  since  there  is  usually  no  bed  gradation. 
However,  it  can  be  calculated  using  procedures  such  as  the  Scour  and 
Deposition  program  in  Appendix  VII. 

One  technique  for  calculating  land  surface  erosion  is  the  universal 
soils  loss  equation  presented  in  Chapter  2 of  this  volume.  The  major 


disadvantage  of  this  technique  Is  that  transport  of  the  sediment  materiel 
Is  not  Included  In  the  formulation.  Yet*  one  application  of  the 
equation  quickly  demonstrates  that  erosion  of  the  land  surface  Is  not 
synonymous  with  sediment  yield  from  the  basin.  Therefore,  an  additional 
factor  called  the  sediment  delivery  ratio  Is  necessary  In  order  to 
convert  from  land  surface  erosion  to  sediment  yield  at  the  outflow 
point.  All  of  the  considerations  of  sediment  transport  are  Included 
In  this  sediment  delivery  ratio  factor.  A better  procedure  would  be 
to  couple  this  soil  erosion  equation  with  the  sediment  transport  models 
and  eliminate  the  need  for  the  sediment  delivery  ratio  factor. 

The  use  of  sediment  detention  basins  requires  analytical  techniques 
that  permit  the  design  of  the  size  of  the  basin  to  effect  a prescribed 
trap  efficiency.  The  detention-time  trap  efficiency  concept  Is 
appropriate  for  deposition  In  the  sediment  basin*  however*  grain  size 
of  the  sediment  material  Is  an  Important  consideration*  (reference  9). 

Section  6.07.  Local  Scour  at  Structures 


The  difference  between  scour  and  degradation  lies  In  the  extent* 
and  the  trend  with  respect  to  time*  of  changes  In  the  stream  bed 
elevation.  Laurson  and  Tock  (reference  11)  demonstrated  that  grain 
size  was  not  a major  factor  In  the  magnitude  of  scour  around  bridge 
piers  and  abutments*  but  rather  the  dimensions  of  the  pier*  Its  alignment 
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with  respect  to  the  flow  and  its  shape  were  the  important  parameters. 
Their  basic  relationship  is  shown  in  the  following  figure,  and  their 
concept  is  presented  in  the  next  few  paragraphs. 


O l.O  2.0  3.0  4.0  5.0  6.0 

Dapth  of  flow 

Width  of  Pior  ot  Stroom  Bod 


Peg.  6.04.  Scour  Around  A Pier 

In  using  this  design  curve,  the  following  information  on  angle 
of  attack  is  pertinent: 
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Local  scour  at  piers.  All  of  the  available  data  were 


adjusted  to  scour  conditions  at  a zero  angle  of  attack 
for  a rectangular  pier  and  plotted  non-dimensional ly 
as  equilibrium  depth  of  scour  against  depth  of  flow 
with  the  width  of  pier  at  stream  bed  used  as  the 
repeating  variable.  The  design  curve  was  drawn  somewhat 
conservatively  and  represents  the  predicted  local  depth 
of  scour  for  a rectangular  pier,  the  sides  of  which  are 
aligned  with  the  current,  as  a function  of  the  depth  of 
flow  and  the  width  of  the  pier  at  stream  bed. 

In  the  laboratory  Investigation  It  was  found  that  the 
angle  between  the  axis  of  the  pier  and  the  direction 
of  the  current  Is  more  important  than  any  other 
geometric  detail  of  solid  or  webbed  piers.  If  the 

« 

pier  Is  skewed  to  the  current,  the  scour  depth  predicted 
must  be  multiplied  by  a factor  greater  than  unity. 

Values  of  this  factor  are  presented  as  a function 
of  the  angle  of  attack  and  the  length-width  ratio  of 
the  pier  at  stream  bed."  (Reference  11,  p.  42) 

Angle  of  attack  and  pier  shape  Information  are  shown  In  the 
following  figure  and  table. 
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PROFILE 

Fig.  6.06,  Scour  at  Bridge  Piers  and  Other  Obstructions 
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Consider,  as  an  example  application,  fig.  6.06  which  shows  a round 

nose  pier  having  a width  of  four  feet  and  located  at  an  angle  of  24^ 
to  the  approaching  flow.  The  pier  Is  36  feet  long.  The  water  depth, 

corresponding  to  a flood  frequency  of  25  years  recurrence  Interval,  Is 

12  feet. 

Entering  fig.  6.04  with  y/b  equal  to  3.0  results  In  a ratio  of 

d^/b  equal  to  2.1.  Adjusting  for  pier  angle  requires  entering  fig.  6.05 

at  24®  with  a L/b  ratio  of  9.  The  resulting  correction  factor  Is  2.6. 
The  depth  of  scour  can  now  be  calculated  as 

ds  - K^i  • (dj/b)  . b 

* 2.6  .2.1-4 

- 22  ft.  (6-02) 

Note  that  the  pier  shape  coefficient,  K^,  Is  only  utilized  when 
the  angle  of  attack  Is  zero. 

This  water  depth  does  not  necessarily  cause  the  most  severe 
condition  of  scour.  In  fact,  a range  of  depths  should  be  tested. 

The  following  table  presents  the  probability  of  at  least  1 occurrence 
of  the  design  flood  during  the  design  life  of  a structure. 

(Reference  11,  page  46). 

Table  6.02.  Chance  of  1 -Occurrence  of  the  Design  Flood 


Design 

Life, 

Years 

25 

Exceedance  Frequency,  Years 
50 

100 

25 

.64 

.40 

.22 

50 

.87 

.64 

.39 

100 

« 

00 

00 

.63 
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These  data  are  useful  for  Interpreting  the  relative  severity  of  the 


scour  problem  when  a range  of  floods  is  utilized. 

Turbulence  In  the  flow  has  a major  Impact  on  scour.  However*  just 
the  effect  of  the  pier  on  concentration  of  flow  Is  substantial.  'When 
the  unit  discharge  Increases  due  to  deflection  of  flow  around  the 
pier*  point  velocities  also  Increase.  This  redistribution  of  flow 
contributes  substantially  to  the  local  scour  problem.  The  design  curves 
presented  earlier  contain  the  Impact  of  both  turbulence  and  the 
redistribution  of  flow  on  scour  depth. 

The  Scour  and  Deposition  model  In  Appendix  VII  Is  not  appropriate 
for  use  In  determining  the  scour  around  bridge  piers.  It  Is  useful* 
however*  for  establishing  the  stream  bed  elevation  between  the  bridge 
abutments.  The  magnitude  of  local  scour  due  to  the  presence  of  the 
pier  Is  In  addition  to  and  measured  relative  to  this  general  bed 
elevation  between  the  bridge  abutments. 
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Shortest  axis  of  a sediment  particle. 

Cross  sectional  area  of  channel  floM. 

Intermediate  length  axis  of  a sediment  particle  or 
pier  width  or  diameter  for  local  scour  calculations. 

Top  width. 

Longest  axis  of  a sediment  particle. 

Factor  for  converting  ppm  to  mg/1. 

Drag  coefficient. 

Cropping  factor  In  the  Universal  Soil  Loss  equation. 
Concentration  of  suspended  sediment,  mg/1. 

Sieve  diameter  of  a sediment  particle. 

Water  depth. 

The  particle  diameter  for  which  90  percent  of  particles, 
by  weight.  In  the  sample  are  finer. 

Depth  of  local  scour  at  a bridge  pier. 

Dominant  water  depth  calculated  for  obtaining  the  dominant 
water  discharge. 

The  grain  size  at  which  50  percent,  by  weight,  of  the  sediment 
grains  are  smaller. 

Acceleration  of  gravity. 

Bed  load  transport  In  tons/day.  dry  weight. 

Total  bed  load  transport  In  metric  tons/sec.  or  Ibs/sec. 
submerged  weight. 

Bed  load  transport  In  metric  tons/sec/meter  of  width, 
submerged  weight. 

Sage  height. 

Coefficient  fbr  converting  volume  per  second  to  weight  per 
day. 

Sediment  load  weighted  by  flow  depth  and  class  Interval  In 
dominant  discharge  calculations. 


r' 


j 


K' 


•s 


L 


m 


Consolidation  coefficient  for  silt  deposits. 

Consolidation  coefficient  for  clay  deposits. 

Coefficient  for  angle  of  attack  between  flow  and  bridge 
pier  for.  analyzing  local  scour. 

1/3 

Strickler  roughness  coefficient*  m /sec. 

Soils  erodablllty  In  the  Universal  Soil  Loss  equation. 

1/3 

The  grain  roughness*  m ' /sec. 

Coefficient  for  pier  shape  for  analysis  of  local  scour. 

Length  term  specifically  defined  where  used. 

Time  Interval  In  days  during  which  the  stage  was  within 
class  Interval  Ah  In  the  dominant  discharge  calculations. 

Manning  n-value  for  grain  and  form  roughness  In  channel. 

Grain  roughness. 

A dimensionless  coefficient  defined  by  equation  3-08. 
Erosion  control  factor  In  the  Universal  Soil  Loss  equation. 
Parts  per  million. 

Water  discharge. 

That  portion  of  the  total  water  discharge  which  Is 
responsible  for  bed  load  transport. 

Dominant  discharge. 

A dimensionless  coefficient  In  the  Meyer-Peter  and  Muller 
equation*  defined  by  equation  3-13. 

Sediment  load. 

Inflowing  sediment  load. 

Outflowing  sediment  load. 

Reynolds  number. 
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r 


Rg  Hydraulic  radius  of  flour  on  channel  bed. 

Rainfall  energy  In  the  Universal  Soil  Loss  equation. 

S Energy  slope  or  storage  capacity. 

Detention  time  or  flow  through  time. 

Trap  efficiency. 

V Average  velocity  of  flow. 

V.  Total  sediment  delivered  to  a reservoir  during  the  project 

life  (50  or  100  years). 

w Fall  velocity  of  a sediment  particle. 

Y Potential  sediment  yield  from  land  surface  erosion  of  a 

® test  plot. 

Average  annual  water  yield. 

Y Specific  weight  of  fluid, 

T*  Initial,  submerged  unit  weight  of  a silt  deposit. 

y"  Initial,  submerged  unit  weight  of  a clay  deposit. 

Y.  Specific  weight  of  a sediment  particle. 

Yj"  Submerged  specific  weight  of  sediment  particles  (y^-y), 

Ysand  Submerged  unit  weight  of  a sand  deposit. 

Yj  Representative  unit  weight  of  a sand/si It/cl ay  deposit  T 

. years  after  deposition  occurred. 

Ah  The  class  Interval  assigned  to  water  stage  In  the  dominant 

discharge  calculations. 

V Kinematic  viscosity  of  fluid, 

p Density  of  the  fluid. 

Density  of  a sediment  particle. 

0 (p,-p)/p* 

T Tractive  fbrce  on  stream  bed. 

: 
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CORPS  OF  ErlGINEER  METHODS  FOR  PREDICTIHG  SEDIM.EHT  YIELDS 


By  Robert  H.  Livesey  _]_/ 

IHTRODUCTIO.T 

The  methods  used  by  the  Corps  of  Engineers  for  predicting  sediment 
yields  are,  in  general,  based  upon  empirical  relationships  but  vary  in 
scope  and  procedure  depending  upon  the  complexity  of  the  individual 
water  resource  project  plan  or  design.  Due  to  the  diverse  nature  of 
these  projects,  both  in  design  magnituce  and  geographic  location,  a 
standard  method  approach  for  design  application  is  not  er.ployed  ttirough- 
out  the  Corps.  Instead,  the  individual  District  Offices  make  a sensi- 
tivity appraisal  to  evaluate  the  impact  of  all  sedimentation  influences 
on  a specific  project  plan.  From  this  first  approxiPiation  analysis,  the 
scope  of  the  sedimentation  problem  is  defineo.  This  cefinicion  then 
becomes  the  basis  for  selection  of  feasible  metiiods  to  be  used  in 
establishing  the  true  ragnitude  of  the  problem  coriponents  and  design 
solution  criteria.  Where  it  is  apparent  that  modification  of  a nrfjtnod  mig.it 
be  practical  to  produce  an  improvement  in  design  evaluation,  sucit  modifica- 
tion is  encouraged.  For  this  reason,  a variety  of  proceuures  are  developed 
and  employed  throughout  the  Corps  but  they  all  relate  closely  to  one 
of  the  three  basic  empirical  approaches  for  predicting  sediment 
yield;  namely  (1)  pleasuring  the  yield  rate  directly  by  sediment  sampling 
or  reservoir  surveys,  (2)  extrapolation  of  such  measured  uata  to  un- 
measured drainages  by  various  correlation  and  probability  techniques, 
or  (3)  establishment  of  identifiable  physiographic  watershed  or  stream 
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flow  characteristics  that  permit  development  of  predictive  equations. 
Theoretical  approaches  to  the  prediction  of  sediment  yields  have  been 
occasionally  employed  for  special  circumstances  where  empirical  relation- 
ships were  weak  or  confidence  lacking  but  such  procedures  are  not 
common. 


DEVELOPMEUT  OF  PREDICTION  METHODS 

The  earliest  record  of  sediment  sampling  in  the  United  States  dates 
back  to  1838  when  the  Corps  of  Engineers  was  engaged  in  navigation  channel 
work  on  the  lower  Mississippi  River.  During  the  next  100  years,  the 
need  for  sediment  predictions  related  almost  entirely  to  river  navigation 
or  estuary  maintenance  work.  It  was  not  until  after  passage  of  the 
Flood  Control  Acts  of  1928  and  1936,  when  the  Corps  started  the  planning, 
design  and  construction  of  multiple  purpose  reservoirs,  that  the  need 
for  sediment  yield  predictions  developed.  Typical  of  this  initial 
phase  of  sediment  yield  investigations  was  STRAUB'S  work  that  is  well 
documented  in  the  1933  Missouri  River  basin  308  Report.  His  development 
of  the  sediment  rating  curve  method  was  later  amplified  by  CAMPBELL  and 
BAUAER  in  the  1940's,  and  MILLER  in  the  I950's,  into  the  popular  rating 
curve-flow  duration  method.  After  STRAUB'S  308  work  the  emphasis  on 
documenting  sediment  yield  rates  shifted  in  the  1940's  to  reservoir 
survey  measurements  and  the  relation  of  sediment  yield  to  contributing 
drainage  areas,  reservoir  capacities,  stream  density  or  slope,  and 
runoff.  The  early  work  of  BROWN  and  GOTTSCHALK  is  typical  of  this 
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period.  But  this  work,  like  STRAUB'S,  was  considered  professionally 
weak  because  It  related  sediment  yield  to  only  a few  of  the  n-.auy 
contributing  factors.  Next,  during  the  early  ISGO's,  efforts  were 
concentrated  on  the  expansion  of  MUSGKAVt'S  definition  of  quantitative 
factors  for  small  land  units  to  the  drainage  increments  of  large 
river  control  projects.  These  evaluations  atteii^pted,  without  much 
success,  to  relate  many  of  MUSGRAVE'S  factors  on  a regional  or  annual 
basis  In  lieu  of  local  or  seasonal  definitions.  During  this  same 
period,  sediment  sampling  and  reservoir  survey  measurement 
techniques  were  enhanced.  Long  term  basin  runoff  characteristics  were 
also  Identified  to  Improve  confidence  In  the  sediment  rating  curve  - 
flow  duration  method.  However,  by  the  late  1950's,  a shift  In  project 
planning  to  smaller  drainage  areas  started.  The  definition  of  local 
drainage  controls  and  urban  runoff  assumed  greater  importance;  the 
"big  dam"  criteria  for  yield  predictions  was  no  longer  vogue.  Tnis 
change  required  a downward  extrapolation  toward  t.-ie  upper  limlis  of  SCS 
criteria.  To  meet  this  need,  the  number  of  Corps  sedii'-ent  load  stations 
were  doubled,  plans  were  Implemented  to  document  urban  runoff 
characteristics  and  correlation  techniques  concentrated  on  qualifying  the 
adequacy  of  short  term  records.  As  the  environmental  Issues  of  the 
late  1960's  developed  their  Impetus,  design  criteria  and  needs  mush- 
roomed Into  the  broad  fields  of  v/ater  quality  control,  biological  repro- 
duction, eutrophication  acceleration  and  nost  recently  v/astewater  management. 
The  proper  prediction  methods  for  most  of  these  latter  aspects  renain 
unqualified  at  the  present  time. 
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PREDICTION  CRITERION 


Corps'  study  investigation  or  design  criterion  for  individual 
projects  dictates  that  alternative  methods  for  sediment  yield 
prediction  be  considered  and  evaluated.  These  empirical  methods  generally 
fall  into  two  basic  categories:  (1)  the  extrapolation  of  measured 
records  and  (2)  the  use  of  predictive  equations.  Most  work  related  to 
the  planning,  design  and  operation  of  reservoirs  is  based  upon  the 
former  method  while  channel  alignment  and  stabilization  aspects  relate  to 
the  latter  niettiod.  Although  the  Corps  has  no  standard  method,  the  use  of 
the  sediment  rating  curve  - flow  duration  technique,  or  some  ramification, 
has  had  the  widest  application.  Prior  to  discussing  the  characteristics 
of  various  methods,  some  comment  must  be  focused  on  the  appraisal  tech- 
niques which  precede  the  selecting  of  the  alternative  methods  to  be 
considered.  The  following  are  common  project  evaluation  criteria: 

1)  Sensitivity  Evaluation.  This  appraisal  attempts  to  bring 
the  scope  of  sedimentation  influences  into  focus  with  the  over-all 
project  purpose  and  plan.  Later  detailed  analyses  define  the  real 
magnitude  and  occmsI orally  dictate  a shift  in  study  emphasis;  but, 
early  in  both  the  planning  and  design  phases,  all  potential  problems  are 
identified  and  priorities  established  regarding  their  importance  to  the 
various  design  features.  When  considering  current  environmental  aspects, 
it  is  not  unconnton  to  schedule  a later,  second  sensitivity  evaluation 

in  case  a re-orientation  of  project  priorities  or  efforts  is  necessary. 

2)  Identifying  Basic  Data.  This  operation  consists  of  a records 

and  literature  search  to  Identify  available  streamflow  or  sediment  measure- 
ment records,  previous  related  study  data  or  pertinent  research  activity 
work. 
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3)  Hydrology.  This  analysis  usually  constitutes  the  initial  design 
effort.  Early  in  this  work  preliminary  sediment  yield  predictions  are 
made  to  satisfy  first  approximation  storage  and  flow  routing 
requirements.  Later,  as  alternate  project  operation  schemes  are  developed, 
the  design  yield  rates  are  incorporated  into  the  final  hydrologic  analysis. 

The  development  of  long-term  flow-duration  data  is  another  important 
contribution. 

4)  Geotnorphic  Characteristics.  The  factors  considered  in  this 
analysis  are  quite  varied  and  broad;  the  degree  of  investigation  depends 
upon  need.  General  needs  include  such  items  as  geologic  variations, 
soil  classifications  and  characteristics,  channel  dimensions,  composi- 
tion of  stream  bank  and  bed  materials  and  stream  slopes. 

5)  Basin  Reconnaissance.  At  least  one  field  reconnaissance  of 
the  contributing  drainage  is  made  early  in  the  preliminary  design  phase. 

Usually  it  is  a "wind  shield"  type  survey  of  sufficient  scope  to  identify 
major  or  contrasting  features  of  the  drainage  and  permit  judgment 
comparison  of  variations  in  yield  rates  and  channel  dimensions.  J 

PREDICTION  METHODS  | 

The  following  discussion  of  the  various  sediment  yield  prediction 
methods  used  in  Corps  of  Engineer  studies  and  designs  will  be  divided 
into  the  two  basic  categories  previously  mentioned.  Descriptive  conments 
will  be  limited  to  a summary  nature  but  a study  or  design  nemorandum 
reference  is  given  for  each  method,  fiost  references  can  provide  corplete 
details  on  one  or  more  methods.  An  attempt  v/as  made  to  include  with 
each  discussion  a reference  plate  which  Identified  the  salient  features  of 
the  method  or  the  output  results. 

I 
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The  first  category  involves  the  extrapolation  of  measured  records 
and  is  divided  into  the  three  major  measurement  classifications: 
sediment  loads,  reservoir  surveys  and  reconnaissance  inspections. 
Relevant  methods  for  each  include: 

1 ) Sediment  Load  Measurements 

a)  Sediment  Rating  Curve  Method.  This  basic,  older  method  is 
usually  associated  with  a flow-duration  analysis  but  occasionally  special 
circumstances  still  require  its  use.  An  example  would  involve  instan- 
taneous units  of  flow  and  concentration  rather  than  mean  daily  values. 
These  applications  usually  relate  to  a near  constant  or  limited  range  of 
flows,  such  as  for  seasonal  or  monthly  variations  between  run-of-river 
reservoirs  within  a large  system.  In  such  instances  the  minor 
incremental  flow  and  sediment  contributions,  including  their  duration 
and  frequency,  are  usually  obscured  by  the  large  base  flow.  The  method 
involves  the  plotting  of  measured  suspended  sediment  load  values  versus 
equivalent  units  of  discharge  for  desired  time  periods  and  defining 

the  mean  curve.  An  example  is  shown  in  Figure  1.  The  original  use  of 
this  method  was  developed  for  the  1933  Missouri  Dasin  30S  report  with 
further  enhancement  by  CAMPBELL  and  BAUDER  in  their  paper,  A Rating- 
Curve  Method  for  Determining  Silt-Oischarge  of  Streams  as  published 
in  Transactions  American  Geophysical  Union,  Volume  30,  August  1949. 

b)  Sediment  Rating  Curve-Flow  Duration  Method.  This  popular 
method  combines  the  above  rating  curve  principle  with  the  measured 
stream-flow  record  to  develop  a probability  correlation  between  the 
sediment  and  water  discharge  of  a stream.  The  method  consists  of  a 
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determination  of  suspended  sediment  load  values  from  the  rating  curve  for 
corresponding  increments  of  discharge  from  a flow  duration  curve. 
Multiplication  of  the  suspended  sediment  load  and  discharge  increments 
by  the  time  percentage  interval  gives  a daily  occurrence  value. 

Totaling  these  daily  average  values  produces  the  mean  daily  discharge 
and  suspended  sediment  load  for  the  year.  Further  multiplication  of 
these  mean  daily  values  by  the  number  of  days  in  the  year  gives  average 
annual  rates.  Addition  of  unmeasured  suspend'  1 or  bed  load  estimates 
results  in  a total  average  annual  sediment  load  value.  Variations  in 
this  method  permit  development  of  long-term  rate  estimates  based  upon 
seasonal  or  short  periods  of  record.  Figures  1,  2 and  3 record  the 
principal  details  of  this  method.  For  more  complete  details,  including 
an  evaluation  of  the  techniques  of  this  iiiethod,  see  An  Analysis  of  the 
Flow-Duration,  Sediment-Rating  Curve  Method  of  Computing  Sediment  Yield 
by  Carl  R.  Miller  as  published  by  the  Bureau  of  Reclamation  in  April  1951. 

c)  Sediment  Discharge-Soil  Type  Relationships.  This  method 
relies  upon  a runoff-sediment  load  record  to  obtain  a correlation  of 
sediment  yield  according  to  soil  classification  and  cultivated  areas. 

River  basins  are  divided  by  soil  types  and  annual  surface  runoff- 
sediment  discharge  curves  developed  for  each  classification  according  to 
the  degree  of  cultivated  acreages.  An  example  of  the  relationship 
developed  for  13  drainages  of  mixed  loess  and  glacial  soils  is  shown 
in  Figure  4.  A comparable  correlation  was  possible  for  residual 
limestone,  sandstone  and  shale  soils  but  in  loessial  terrain  the  results 
were  indeterminate.  For  further  information  refer  to  an  utipublisned  report 


Hill. 


d)  Dominant  Basin  Characteristics.  The  similarity  of  the 
dominant  physical  characteristics  of  a drainage  basin  versus  the  measured 


sediment  production  is  the  basis  for  this  method.  The  dominant 
characteristics  included  land  use,  relief  and  topography,  climate, 
water  and  soil  types.  Land  resource  areas  are  used  to  group  the 
defined  individual  sediment  yield  rates  into  comparable  area  categories. 
Both  suspended  sediment  load  and  reservoir  sedimentation  survey  records 
are  used  to  establish  yield  rates  by  drainage  area  or  time  Increments 
for  a given  base  period.  The  flow-duration  principle  is  applied  to 
adjust  short-term  records  to  the  base  period.  Such  adjustments 
require  establishment  of  sediment  discharge  to  streamflow  relationships 
for  the  period  of  measurements  and  then  correlating  this  data  to  the 
long-term  flow  regimen  of  the  stream.  The  method  has  produced 
indications  of  sediment  yield  trends  with  time  in  several  instances. 
Figures  5,  b and  7 depict  the  general  features  of  this  method.  For 
details  see  Sediment  Yields  in  the  Upper  Mississippi  River  Basin  by 
Frank  0.  Mack  as  published  in  Proceedings  of  a Seminar  on  Sediment 
Transport  in  Rivers  and  Reservoirs  in  April  1970  by  the  Hydrologic 
Engineering  Center,  Corps  of  Engineers  at  Davis,  California. 

e)  Sediment  Yield  by  Isogram  Intervals.  Except  for  the 
degree  of  individual  basin  analysis,  a similarity  exists  betvfeen  this 
and  the  preceding  method.  This  method  recognizes  the  dominant  physical 
characteristics  and  measured  sediment  production  records  of  the  basin, 
but,  in  addition,  relies  upon  personal  knowledge  and  engineering 
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judgment  to  evaluate  tlie  sediment  yield  characteristics  of  a basin. 

The  method  was  developed  for  use  as  a task  force  expedient  by  a group 
of  inter-agency  sedimentation  specialists  to  document  seditrent  yield 
rates  for  large  river  basins.  Yield  rates  for  standara  periods  of 
time  are  derived  by  extrapolation  of  shorter  period  records  by  one 
of  three  procedures;  comparing  sediment  load-water  discharge  relations 
between  periods  of  record  and  the  standard  period,  (derivation  of  sediment- 
water  regression  curves  for  increments  of  drainage  area  or  evaluating 
relations  between  intermittent  sediment  measurements  made  over  short- 
time  periods.  The  final  delineation  of  isogram  lines  are  based  upon 
group  experience  and  judgment.  A typical  end-proouct  of  this  method 
is  shown  on  Figure  1),  Examples  of  this  method  can  be  found  in  any 
one  of  the  seven  sub-basin  sedimentation  reports  prepared  by  the  Task 
Force  on  Sedimentation  for  the  fiissouri  River  basin  Comprehensive 
Framework  Study,  as  submitted  to  the  Missouri  Basin  Inter-Agency 
Comnittee  in  1968  and  1969. 

f)  Sediment  Yield  Durino  Urban  Expansion.  The  techniques  of 
this  metliod  are  still  in  the  developmental  stage.  The  basic  premise 
concerns  transition  of  rural  lands  to  urban  usage  over  given  time 
periods.  It  assumes  that  sediment  yield  rates  accelerate  from  agricul- 
tural values  to  a high  peak  during  landscaping  or  construction,  then 
decline  to  a lower  plateau  as  the  land  "heals"  and  finally  level  off 
at  some  low  stable  rate  representative  of  business  or  residential 
lands.  A projection  of  urban  expansion  limits,  as  provided  by  the 
local  metro  planning  authority,  serves  as  a base  for  converting 
contributing  drainages  from  rural  use  to  single  family,  multi -family 
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or  cormercial  usage.  Integration  of  varying  yield  rates  for  area 
increments  under  various  stages  of  development  permits  a continuous 
assessment  over  the ‘design  life  of  the  project.  Judgment  extrapolation 
of  limited  urban  runoff  and  sediment  yield  measurements  is  currently 
necessary  but  data  collection  programs  that  concentrate  on  stonn 
runoff  measurements  can  quickly  improve  this  limitation.  A generalized 
schematic  outline  of  this  method,  as  being  developed  by  the  Omaha 
District,  is  shown  in  Figures  9 and  10. 

2)  Reservoir  Sedimentation  Surveys 

a)  Sediment  Yield  Per  Unit  of  Drainage.  The  application  of  this 
method  is  common  because  of  its  simplicity  in  relating  measured  rates  of 
sediment  yield  to  the  contributing  drainage  area  increrijent.  dumerous 
correlations  are  possible  witiiin  certain  ranges  of  drainage  area  by 

soil  types,  runoff  volumes,  watershed-capacity  ratios,  dominant  discharge, 
land  use,  physiographic  areas  and  many  other  parameters.  Most  Corps 
applications  of  these  yield  rates  pertain  to  contributing  drainages 
greater  than  100  square  miles  so  correlation  with  the  conventional 
soil  loss  parameters  is  not  conwon.  The  prijiciple  source  of  reference 
data  is  ARS  Mi  sc  Publication  No.  1143,  Summary  of  Reservoir  Sediment 
Deposition  Surveys  Made  in  tiie  United  States  through  196b.  or  related 
reporting  Form  17S7.  A typical  example  of  this  method  can  be  noted  in 
Figure  11. 

b)  Yield  Production  for  Debris  basins.  This  is  a special 
application  used  to  determine  the  sediment  yield  into  flood  control 
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debris  basins  in  mountainous  type  terrain.  The  method  was  developed 
from  observed  debris  volumes  that  reflect  ground  conditions  influenced 
by  prior  rain  runoff  and  areas  subjected  to  partial  or  complete 
''burns.''  Influencing  factors  include  size  and  shape  of  drainage  area; 
steepness  of  canyons  and  side  slopes;  geological  characteristics;  type 
and  density  of  plant  cover;  recency  of  burns;  and  frequency,  duration 
and  intensity  of  storms.  Measured  debris  volumes  are  adjusted  to  a 
conrion  base  and  curves  developed  for  separate  corrections  of  the  major 
factors  affecting  debris  production.  Figure  12  summarizes  the  details 
of  this  niethod.  Further  information  is  available  in  A New  rifethod  of 
Estimating  Debris-Storage  Requirements  for  Debris  3asins  by  Fred  E.  Tatum 
as  published  in  the  Proceedings  of  the  1963  Federal  Inter-Agency 
Sedimentation  Conference,  ARS  Hisc.  Pub.  No.  970. 

3)  Reconnaissance  Inspections.  The  follov/ing  methods  are 
directed  toward  establishing  preliminary  estimates  of  sediment  yield  for 
large  drainage  areas.  On  occasion,  the  investigation  details  have  been 
expanded  to  cover  studies  of  design  scope  for  small  to  moderate  drainages. 
Their  basic  premise  consists  of  a quick  but  detailed  reconnaissance 
type  inspection  of  the  contributing  drainage  area  by  two  or  ntore  Sed- 
imentation Specialists  who,  by  experience,  are  capable  of  making  judgment 
estimates  of  sediment  yield  rates.  During  the  field  reconnaissance 
they  collectively  establish  representative  point  rates  for  incremental 
portions  of  major  drainages  within  the  over-all  study  basin.  This 
technique  is  particularly  applicable  for  a degree  assessment  of 
contributing  versus  non-contributing  drainage  as  influenced  by  soil 
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management  practices,  smaller  reservoirs  or  ponds  or  Irrigation 
diversion  projects.  If  the  basin  Is  relatively  small,  perhaps  less 
than  1000  square  mile,  the  estimates  for  even  third  or  fourth  order 
streams  can  become  quite  detailed.  For  larger  basins,  selected  streams 
might  be  covered  In  more  detail  and  the  remainder  left  to  a 
random  choice  of  Inspection.  The  end  product  Is  usually  similar  to 
that  shown  In  Figure  6. 

a)  Interpolation  of  Rates  Within  a Basin.  This  method 
requires  several  points  of  measured  sediment  yield,  by  either  sediment 
sampling  or  reservoir  surveys,  within  the  basin  drainage.  One  of  these 
points  should  be  located  near  the  mouth  of  the  basin  to  reflect  the 
total  measured  yield  from  the  drainage.  During  the  field  reconnaissance 
these  measured  rates  are  used  as  a comparative  guide  for  estimating 
yield  rates  for  small  Increments  of  the  unmeasured  drainages.  When 
sufficient  point  estimates  are  established,  a yield  contour  map  Is 
developed.  Using  digitizing  or  planimetering  processes,  drainage  area 
Increments  of  equal  yield  rates  are  totaled  for  the  major  drainages 
within  the  basin.  A summation  of  these  totals  and  division  by  the 
contributing  drainage  area  value  gives  an  average  sediment  yield  rate 
for  the  subject  Increment.  These  Increment  rates  are  checked  against 
the  measured  Increment  rates  for  verification.  If  they  are  not  com- 
parable within  reason,  adjustments  to  selected  point  estimates  are 
Justified  to  bring  the  Integrated  total  Into  balance  with  measured 
data. 

b)  Extrapolation  to  Unmeasured  Watersheds.  The  basic 
procedure  Is  similar  to  that  above  except  that  a comparison  between  the 
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total  estimated  and  measured  rates  for  a basin  Is  not  possible.  Prior 
to  the  field  Inspection  of  the  unmeasured  drainage,  the  reconnaissance 
team  usually  makes  a preliminary  Inspection  of  the  measured  drainages 
being  used  as  the  extrapolation  base.  This  visual  Inspection  requires 
additional  time  and  effort  but  serves  as  an  effective  means  for 
comparative  extrapolation.  The  validity  of  this  method  Is  dependent 
upon  the  degree  of  extrapolation  but  apparent  satisfactory  results  have 
been  produced  within  a restricted  time  period. 

The  second  category  Involves  predictive  equation  methods.  Most 
of  the  Individual  methods  discussed  below  apply  to  the  solution  of 
specific  problems.  They  differ  from  the  preceding  methods  In  that 
the  predicted  sediment  yield  relates  primarily  to  channel  contributions 
rather  than  from  a watershed  drainage.  The  Corps  use  of  predictive 
equations  for  determining  watershed  yields  Is  very  limited. 

1)  Sediment  Transport  Relationships.  There  are  a 
variety  of  methods  In  this  classification  but  the  most  common  Is  the 
EINSTEIN  approach,  with  one  of  Its  many  modifications,  or  the  more 
recent  TOFFALETI  procedure.  Their  connection  with  sediment  yield 
predictions  usually  relate  to  channel  stabilization  projects  Involving 
aggradation  or  degradation  problems.  But  their  application  Is  also 
connon  In  establHhlng  the  magnitude  or  rate  of  unmeasured  suspended 
or  bed  sediment  load  values.  Estimates  of  such  values  are  extremely 
Important  In  certain  Instances  when  establishing  yield  rates  from 
measured  suspended  sediment  load  records,  such  as  Is  required  In  the 
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rating  curve-flow  duration  method.  An  excellent  discussion  of  the 
ElfiSTEIM  and  TOFFALETI  methods,  plus  others,  and  a listing  of  complete 
references  can  be  found  In  the  Journal  of  the  Hydraulics  Division, 

ASCE,  April  1971,  under  Paper  C07f-  titled  Sediment  Transportation 
Mechanics:  Sediment  Discharge  Formulas 

2)  Uetentlon-TIme  Method.  This  method  is  used  to 
predict  the  volume  of  sediment  passing  through  a series  of  run-of-the- 
river  type  reservoir  projects.  It  is  based  upon  empirical  relationships 
between  tlie  detention  time  of  flows  passing  through  the  reservoir  system 
and  the  percentage  of  sediment  deposited.  Detention  time  is  defined  as 
the  ratio  of  reservoir  storage  to  the  Inflow  discharge  at  any  given 
tinie.  Curves  of  detention  time  versus  percent  of  wash  load  and  bed 
material  load  deposited  are  developed  to  determine  deposition-flow 
through  volumes  for  incremental  reaches.  As  the  reservoir  volume  Is 
depleted  by  deposition,  the  detention  time  is  reduced  and  the  yield 
rate  per  unit  of  flow  increases.  Figure  13  shows  a typical  relationship 
between  detention  time  in  liours  and  percent  of  load  deposited  for  a 
series  of  reservoirs.  Reference  details  for  this  method  can  be  found 

in  Dardanelle  Reservoir  Design  Memorandum  Ho.  6,  Part  IV  Sedimeritatlon 
prepared  in  October  1957  by  the  Little  Rock  District. 

3)  Hydraulic  Elements.  This  is  another  method  for 
detemlning  run-of-the-river  reservoir  yield  rates.  Individual  reservoirs 
are  divided  into  reaches  and  hydraulic  elements  are  determined  for  various 
discharges  by  back-water  computations.  The  sediment  load  Is  then  related 
to  the  velocity,  depth  and  slope,  or  a combination  of  such  eleiiients.  The 
Inflow-outflow  sediment  loads  are  computed  on  a step  procedure  from  reach 
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to  reach  with  backwater  conputations  repeated  as  necessary  to  obtain 
nevj  hydraulic  eleivent  values  as  deposits  accumulated  with  tine.  The 
method  is  tine  consuming  but  permits  recognition  of  the  individual 
element  changes  or  trends,  an  example  is  shown  in  Figures  14  and  lb. 

(Note:  These  figures  could  not  be  reduced  in  size  and  consequently 
are  not  included.  The  method  depicted  is  embodied  in  the  corr^'Uter  program 
presented  in  Appendix  VII  of  this  volume.  W.  A.  Thomas)  Tne 
reference  for  this  method  is  the  same  as  given  above  for  tiie  detention- 
time  method. 

4)  Bank  Caving-Uash  Load  Yields.  This  is  another  special 
application  method  used  to  estimate  the  excess  rate  of  bank-caving 
over  bank  building  relative  to  an  increase  in  the  wash  load  being 
transported.  It  applies,  again,  to  run-of-river  reservoirs  where  the 
erosion  of  concave  banks  in  bends  nay  continue  at  the  natural  high  flow 
rates  due  to  the  duration  of  artificial  bank  full  stages  under  operating 
conditions.  The  change  in  bank  caving-bank  building  equilibrium 
produces  an  additional  supply  of  wash  load  due  to  channel  widening. 

The  method  is  admittedly  over-simplified  for  the  complex  phenomena 
of  bank  caving  but  it  offers  a systematic  method  of  recognizing  the 
increase  in  wash  load  due  to  bank  caving.  An  example  of  this  method 
is  given  in  Figures  16  and  17.  Reference  details  can  be  found  in  a 
Supplement  To  Dardanelle  Reservoir  Design  Menwrandum  No.  6-4  prepared 
by  the  Little  Rock  District  in  January  1959. 
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5)  Sediment  Delivery  Ratio.  This  category  covers  both 
the  sediment  delivery  and  sheet  erosion  prediction  methods  developed  by 
the  Department  of  Agriculture,  The  application  of  these  methods  to 
Corps  projects  is  generally  limited  to  small  watersheds  of  less  than 

25  to  50  square  miles.  Adoption  of  the  MUSGRAVE  equation  is  probably 
still  preferred  over  the  Universal  Soil  Loss  equation  for  smaller 
drainages.  However,  more  useful  are  the  various  empirical  equations 
developed  by  such  authors  as  ANDERSON.  BARNES,  BRUNE,  GLYMPH,  GOTTSCHALK, 
HEINEMANN,  KOHLER,  MANER,  PIEST  and  others.  Due  to  the  limited  applica- 
tion of  these  methods,  reference  is  made  to  the  following  sources  for 
procedural  details:  Studies  of  Sediment  Yields  From  Watersheds 
by  Louis  Glymph  or  Sediment  Sources  and  Sediment  Yields  prepared  by 
Robert  F.  Piest  as  Chapter  IV  of  the  ASCE  Manual  on  Sedimentation 
Engineering  and  published  as  paper  7337  in  the  Journal  of  the  Hydraulics 
Division,  June  1970. 

6)  Tailwater  Degradation.  Several  methods  are  Included 
in  this  grouping.  Their  principal  function  is  to  predict  degradation 
trends;  but,  as  part  of  the  computational  procedure,  sediment  yield 
values  for  the  degrading  reach  are  developed.  Factors  considered  in 
their  application  include  composition  of  the  bed  material  and  its 
coarsening  with  time,  the  magnitude  of  future  flows  and  changes  in  flow 
characteristics  such  as  channel  shape,  depths,  velocity  and  slope. 

Three  common  methods  include  use  of  the  EINSTEIN  bed  load  function  for 
sediment  transport;  use  of  the  KALINSKE  formula  to  compute  bed  load  plus 
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the  development  of  relationships  between  natural  and  modified  bed  i 

'( 

i 

material  load  transport  as  degradation  progresses;  and  determination 
of  the  thickness  of  an  armor  layer  and  the  depth  of  scour  at  which 
this  layer  would  form.  This  latter  technique  assumes  degradation  would 
cease  when  a layer  of  "non-moving"  particles  forms  a sufficient  armor 
to  prevent  the  effective  leaching  of  finer  particles  from  the  underlying 
bed  material.  Procedural  details  on  these  methods  can  be  found  In  a 
paper  on  Sedimentation  Studies  for  Robert  S.  Kerr  Lock  and  Dam. 

Arkansas  River  Basin  by  Howard  0.  Reese  as  published  In  the  Hydrologic 
Engineering  Center's  Proceedings  of  a Seminar  on  Sediment  Transport  In 
Rivers  and  Reservoirs,  April  1970. 

FUTURE  HEEDS 

During  the  past  decade  a shift  In  emphasis  has  taken  place  within 
the  Corps  regarding  the  need  for  sediment  yield  predictions.  Prior  to 
the  m1d-1950's,  sediment  was  viewed  primarily  as  a malignant  growth 
that  reduced  the  effectiveness  of  reservoirs,  floodways,  navigation 
channels  and  harbors.  This  was  also  the  period  of  "big  dam"  planning 
and  construction  where  sediment  depletion  rates  played  a relatively 
minor  role  In  design  because  of  the  volimilnous  storage  allotted  for 
multiple  purpose  use.  The  need  for  sediment  yield  predictions  for 
large  drainage  areas  has  essentially  vanished.  As  an  Indication,  about 
15  years  ago  the  Corps  was  operating  135  sediment  load  stations  with 
43%  having  drainage  areas  greater  than  5000  square  miles,  27%  In  the 
500-5000  range  and  30%  less  than  500  square  miles.  Presently  this  number 
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of  stations  has  doubled  with  a comparative  shift  In  percentage  ratio 
of  25:37:38.  Almost  half  of  the  current  stations  are  operated  for 
planning  or  design  purposes.  For  example*  during  1969  the  Corps  had 
under  construction  23  reservoir  projects  for  hydro-electric  power  and 
flood  control,  64  reservoir  projects  for  flood  control  and  multi-purpose 
use  and  84  local  flood  control  protection  projects.  Now  the  Corps 
emphasis  seems  to  be  focused  on  projects  with  sediment  contributing 
drainages  that  generally  vary  within  the  500  to  2500  square  mile  range. 
But  If  our  prediction  approach  Is  to  continue  on  an  empirical  basis, 
long-term  data  records  for  drainage  areas  within  this  bracket  are 
Inadequate,  particularly  for  reservoir  survey  data.  It  Is  estimated 
that  there  are  some  28,000  reservoirs  In  the  United  States 
yet  we  have  sediment  yield  records  on  only  4X.  But  more  significant 
Is  the  fact  that  of  the  1200  Individual  reservoirs  listed  In  the 
1965  summary  of  reservoir  survey  data,  80t  of  the  documented 
record  ranges  below  50  square  miles  and  90X  below  500  square  miles.  It 
Is  apparent  that,  figuratively  speaking,  a scarcity  of  data  exists  In 
the  no-man's  land  that  Is  bracketed  with  voluminous  SCS  records  on  the 
low  side  and  adequate  Corps,  USBR  and  TVA  experience  on  the  high  side. 

In  essence,  the  basis  for  future  Corps  yield  predictions  by  empirical 
methods  will  be  somewhat  handicapped  until  data  records  within  this 
bracket  are  expanded  by  measurements  or  enhanced  by  correlation 
techniques. 

TodiO^,  the  dirty  word  "sediment*  has  dual  connotations;  It  must  now 
be  recognized  from  both  a beneficial  and  detrimental  point  of  view. 

On  one  hand  sediments  rank  as  a major  cause  of  water  pollution,  but  on 
the  other  hand  they  play  a dominating  role  In  water  quality  control 
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due  to  their  assimilation  capabilities.  Apparently  they  also  serve 
similar  dual  roles  as  catalytic  or  transporting  agents  in  physical, 
chemical  or  biological  processes.  With  the  current  focus  of  Corps 
activities  in  areas  of  environmental  control,  urban  development  or 
expansion  and  wastev^ater  nanagenent,  the  recognition  of  such  aspects  is 
receiving  prime  attention  in  planning  and  design,  but  unanswered 
questions  continue  to  outnumber  even  qualified  answers.  There  is  an 
unquestionable  need  for  further  amplification  of  key  sedimentation 
influences  in  certain  environmental  processes  before  proceeding  with 
detailed  planning  or  design  applications. 

The  imediate  needs  of  the  Corps  of  Engineers  in  expansion  of 
sediment  yield  prediction  methods  will  probably  be  focused  along  these 
two  major  channels:  1)  definition  of  empirical  relationsliips  for 
drainage  areas  of  moderate  size,  and  2)  establisniient  of  tiie  role 

I 

sediments  play  in  the  complex  environmental  processes.  The  application  | 

of  computer-oriented  nietliods  for  mathematical  sii.iulations  or  modeling  | 

will  undoubtedly  play  a key  role  in  tiie  solution  of  some  of  these 

problems.  Past  experience  tiowever,  has  clearly  demonstrated  that  one 

or  two  standard  methods  or  universal  equations,  regardless  of  their 

complexity,  will  not  meet  the  diverse  needs  of  modern  day  engineering, 

planning  or  design.  For  this  reason,  our  efforts  will  continue  to 

concentrate  along  the  lines  of  practical  needs  to  resolve  particular 

problems.  But  in  doing  so,  we  intend  to  also  continue  our  policy  of 

improving  available  methods  or  techniques  by  modification,  regardless  of 

their  degree  of  sophistication,  as  the  needs  of  the  problem  warrant. 
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PMTURM.  WASH  LOAD  IN  1,000,000  TONS/YEAR 


NOTE:  This  Information  was  presented  In  Project  Design  Memorandum  No.  6, 
Sedimentation-Part  IV,  Section  II.  Deposition  In  Dardanelle  Reservoir, 
Arkansas  River  and  Tributaries,  Multiple-Purpose  Plan,  Arkansas  and 
Oklahoma,  U.S.  Amy  Engineer  District,  Little  Rock,  Corps  of  Engineers, 
Little  Rock,  Arkansas,  October  1957. 


SEDIMENT  DATA  FOR  KEOKUK  (POOL  HO.  19) 
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and  sediment  range  resurveys  were  made  In  1928,  1938,  and  1946.  By  1946, 
after  33  years  of  operation,  the  reservoir  had  lost  167,600  acre-feet, 
or  35  percent  of  Its  capacity.  The  results  of  these  surveys  are  shown 
In  "Reservoir  Sedimentation  Data  Sunmary"  and  In  Table  12  herein  with 
the  average  Inflow  for  the  periods  between  surveys. 

40.  In  1956  the  Rock  Island  District  spudded  at  eight  locations 
spaced  throughout  the  reservoir  and  obtained  unit  weights  and  grain  size 
distribution  of  the  deposits.  From  those  data  an  average  dry  weight  of 
55  pounds  per  cubic  foot  was  computed.  The  deposits  In  a reservoir 
become  denser  with  age.  Using  the  terminal  unit  weight  of  55  pounds 
per  cubic  foot  and  the  formula  listed  below  by  Lane  and  Koelzer  (2) 

as  a guide,  unit  dry  weights  of  50  to  54  pounds  were  assumed  for  the 
periods  as  shown  in  Table  12.  The  formula  Is: 

Weight  of  sand  * 93  pounds/cu.  ft. 

Weight  of  silt  * 65  + 5.7  log  (No.  of  years  of  deposit) 

Weight  of  clay  *30  + 16  log  (No.  of  years  of  deposit) 

The  size  distribution  of  the  deposits  used  were  10  percent  sand,  53  percent 
silt,  and  37  percent  clay.  The  volumes  of  deposits  were  then  converted  to 
weight  In  tons  as  shown  In  Table  12.  Conversion  of  the  deposits  to 
weight  was  necessary  to  compare  with  the  sediment  load  Inflow  which  was 
of  course  sampled  In  terms  of  percent  of  sediment  by  weight  of  water 
mixture. 

41.  The  next  step  was  to  estimate  the  sediment  Inflow  to  Keokuk 

for  the  different  periods.  The  flow-duration  curves  were  determined  for  the 
periods  as  shown  on  plate  20.  The  Rock  Island  District  has  sampled  the 


suspended  load  at  Burlington,  Iowa,  near  the  head  of  the  Keokuk  pool 
during  the  years  1944  to  1950  and  determined  the  average  annual  suspended 
load.  The  Soil  Conservation  Service  has  also  made  similar  estimates.  An 
average  load  curve  Is  shown  on  Plate  21  which  Is  Intended  to  Include  both 
suspended  and  bed  load.  By  applying  the  f1ow>durat1on  curves  on  Plate  20 
to  the  average  load  curve  on  Plate  21,  the  average  annual  Inflowing  load 
for  each  period  between  range  surveys  was  computed.  As  the  load  curve 
on  Plate  21  Is  based  on  sediment  measurements  after  construction  of  the 
upstream  pools  Nos.  11  to  17,  It  was  necessary  to  correct  the  loads 
for  the  periods  1913-28,  1928-38,  and  1938-46  to  allow  for  the  sediment 
which  was  trapped  In  those  pools  (deposits  In  those  pools  are  described 
In  par.  43).  Finally,  the  silt-clay  loads  for  the  respective  periods 
were  separated  from  the  total  load  as  shown  In  Table  12.  Grain-size 
analysis  of  samples  by  the  Rock  Island  District  Indicate  that  the  load 
Is  composed  of  about  10  percent  sand  and  90  percent  silt-clay.  The 
primary  Interest  In  Keokuk  herein  Is  the  silt-clay  fraction.  The 
percentages  of  Inflowing  silt-clay  load  deposited  In  the  reservoir  for 
the  four  periods  are  shown  In  Table  12. 

42.  In  order  to  develop  a relationship  between  detention  time  and 
percent  of  sediment  deposited  for  the  respective  periods  between  resurveys 
It  Is  necessary  to  know  the  representative  detention  time  to  plot  since 
It  varies  with  each  discharge.  It  was  found  that  the  best  trial-plotting 
position  was  for  a discharge  corresponding  to  the  percent  of  time  when 
about  half  the  load  was  transported,  or  10  percent  of  the  time  for  Keokuk. 
The  detention  times  for  flows  occurring  10  percent  of  the  time  are  shown 
In  Table  12.  The  plotted  points  of  detention  time  versus  percent  of 
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TABU  13 
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the  s1U-c1a^  load  deposited  are  shown  on  Plate  22  with  a curve  drawn 
approximately  through  the  points.  An  example  of  the  trial  computations 
to  check  the  curve  Is  shown  In  Table  13  for  the  1913-28  period.  The 
percent  of  time  and  corresponding  discharges*  sediment  load*  and 
detention  times  are  listed  In  the  table.  The  detention  time  Is  the 
ratio  of  the  average  available  level  storage  for  the  period  divided  by 
the  discharges  as  noted  by  footnote.  The  corresponding  percentages  of 
sediment  deposited  are  obtained  from  the  silt-clay  curve  on  Plate  22  and 
multiplied  by  the  loads  to  obtain  the  deposits.  The  sum  of  these 
deposits  Is  6*601*000  tons  per  year  compared  to  the  estimated  value 
of  7*146*000  tons  from  the  range  survey*  or  a discrepancy  of  about  8 
percent  for  the  1913-28  period.  Similarly*  the  other  three  periods 
were  tested  and  the  computed  deposits  are  shown  In  the  lower  part  of 
Table  12.  They  may  be  compared  with  the  silt-clay  deposits  from  the 
range  surveys  shown  Imnedlately  above.  It  will  be  noted  that  the  silt- 
clay  curve  on  Plate  22  Is  poorly  defined  for  detention  times  greater 
than  about  50  hours.  However*  that  part  of  the  curve  corresponds 
to  the  smaller  discharges  and  therefore  the  smaller  sediment  loads* 
so  that  this  part  of  the  curve  Is  relatively  Insensitive  to  error. 

50.  Garrison  Reservoir.  Garrison  Is  a large  multiple-purpose 
project  on  the  middle  Missouri  River.  On  7-9  July  1954  the  suspended 
load  was  sampled  and  the  hydraulic  elements  were  measured  at  a number  of 
ranges  starting  at  the  head  of  backwater  and  terminating  In  the  lake 
where  most  of  the  load  had  been  deposited.  The  discharge  at  the  time 
was  about  30*000  cubic  feet  per  second.  The  results  of  this  survey  are 
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given  In  reference  (29).  Similar  surveys  were  made  at  Fort  Peck  and  Fort 
Randall.  The  differences  in  the  concentration  at  successive  ranges 
indicate  the  sediment-  load  deposited.  The  percent  deposited  upstream 
from  any  range  may  be  computed  from  the  data  as  well  as  the  storages 
and  detention  times. 

51.  A summary  of  this  information  for  Garrison  Reservoir  is  shown 
in  Table  17.  The  percent  deposits  of  sand  and  silt-clay  are  shown 
separately.  The  detention  times  versus  percent  of  deposits  for  these 
two  size  fractions  are  shown  on  Plate  22,  representing  accumulated  values 
above  the  respective  ranges.  The  plot  of  the  sand  fraction  follows  the 
trend  of  the  data  shown  for  Imperial  Dam.  The  silt-clay  fraction  checks 
the  silt-clay  curve  shown  for  short  detention  times,  but  indicates  a 
generally  steeper  curve  as  will  be  noted  from  the  points  for  detention 
times  of  29  to  52  hours.  The  plotted  values  showing  minus  deposits,  or 
scour,  may  be  the  result  of  error  in  measuring  the  concentration,  although 
scour  will  sometimes  occur  with  short  detention  times.  The  data  are  not 
fully  comparable  to  the  other  information  shown.  Storages  under  the 
backwater  curve  were  included  in  computing  the  detention  time  and  the 
cross  sections  used  omitted  dead  water  areas.  Also,  Garrison  is  a large 
storage  reservoir  and  the  velocities  decrease  rapidly  as  the  flow 
reaches  the  deeper  part  of  the  lake.  This  may  account  in  part  for  the 
steep  trend  of  the  plotted  points,  or  the  more  rapid  fall  out  of  the 
suspended  load.  The  data  for  Fort  Peck  and  Fort  Randall  indicate  a 
similar  trend  and  are  not  shown.  While  this  type  of  information  is 
limited,  it  is  of  considerable  interest  and  contributes  to  the  available 
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records.  Further  reference  is  made  to  the  report  on  these  projects  In 
connection  with  another  method  considered  herein  for  estimated  sediment 
deposits  In  Oardanelle  Reservoir. 
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TABLE  17 
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FORECASTING  DISTRIBUTION  OF 
SEDIMENT  DEPOSITS  IN  LARGE  RESERVOIRS 
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DEPARIMENl  OF  THE  ARMY 
Office  of  the  Chief  of  Engineers 
ENGCW-EY  Washington,  D.C.  20315 

Engineer  Technical 
Letter  No.  1110-2-64 

ENGINEERING  AND  DESIGN 
Distribution  of  Reservoir  Sediment  Deposits 

1.  Purpose  and  Scope.  The  purpose  of  this  ETL  is  to  present  a discussion 
of  problems  associated  with  reservoir  sediment  accumulations  and  to  review 
pertinent  planning  and  design  considerations  relating  to  anticipated  distri- 
butions of  reservoir  deposits.  This  letter  is  applicable  to  all  Divisions 
and  Districts  concerned  with  civil  works  activities. 

2.  References ; a.  Paper  entitled  "Forecasting  Distribution  of  Sediment 
Deposits  in  Large  Reservoirs,"  by  Brice  L.  Hobbs  (Appendix  I). 

b.  EM  1110-2-4000,  "Reservoir  Sedimentation  Investigations  Program," 

IS  November  1961. 

3.  Discussion;  a.  The  Importance  of  forecasting  the  probable  distribution 
of  reservoir  sediment  deposits  is  often  overlooked  in  planning  and  design 
investigations.  Many  recent  project  design  reports  present  only  one  or  two 
brief  paragraphs  concerning  the  gross  volume  of  storage  space  depletion 
anticipated  during  the  period  considered  for  economic  investigation  and 
make  no  reference  to  possible  serious  problems  that  might  be  created  by 
adverse  distributions  of  sediment  deposits  in  particular  areas  or  elevation 
zones.  Conditions  at  certain  reservoirs  which  have  been  in  operation  for 
long  periods  point  up  the  need  for  making  sediment  distribution  studies  in 
connection  with  most  reservoir  design  Investigations;  this  is  particularly 
true  of  all  projects  on  streanu'.  which  transport  substantial  quantities  of 
sediment. 

b.  In  most  of  the  check  coiif>utations  made  in  the  Office  of  the  Chief 
of  Engineers  in  reviewing  the  methodology  and  results  of  field  measurements , 
the  results  obtained  by  the  "Pool-Elevation  Duration  Method,"  presented  with 
Appendix  I,  have  been  found  to  give  results  conparing  favorably  with  the 
measured  values.  However,  there  are  a few  Instances  where  the  values  esti- 
mated for  certain  elevation  zones  depart  appreciably  from  the  measured  values 
reported  in  the  "Reservoir  Sediment  Data  Summaries"  submitted  in  accordance 
with  EM  1110-2-4000  (ref  2b).  This  is  also  true  of  estimates  obtained  by 
other  methods.  The  reasons  for  some  of  these  discrepancies  are  obscure  while 
in  others,  examination  of  the  records  yield  logical  explanations.  For 
example,  it  is  doubtful  that  conditions  of  deposition  in  Jemez  Canyon 
Reservoir  in  New  Mexico  could  have  been  accurately  predicted  by  any  of  the 
currently  available  enpirical  methods,  since  substantial  quantities  of 
material  have  accumulated  in  reservoir-elevation  zones  high  above  the 
maximum  experienced  pool  elevation,  and  it  appears  likely  that  some  of 
the  deposition  is  wholly  unrelated  to  reservoir  effects. 
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4.  Action  to  be  Taken.  The  information  presented  In  i^pendlx  I discusses 
various  types  of  upstream  sedimentation  problems  that  may  be  Induced  by  the 
construction  and  operation  of  a dam  and  reservoir.  The  computational  pro- 
cedures outlined  In  the  report,  and  those  described  In  references  therein, 
are  representative  of  methods  available  for  estimating  the  locations  of 
future  sediment  deposits  In  reservoirs.  The  practices  and  techniques 
described  In  Appendix  1 have  been  found  by  long  experience  to  give  generally 
satisfactory  results  In  estimating  future  sediment  conditions.  Accordingly, 
It  Is  suggested  that  Corps  representatives  utilise  the  methods  and  techniques 
outlined  In  the  Inclosed  report,  where  feasible.  In  connection  with  sediment 
Investigations. 

FOR  THE  CHIEF  OF  ENGINEERS: 


I Appendix 

Forecasting  Distribution  of 
Sediment  Deposits  In  Large 
Reservoirs,  10  Feb  69 


NBNmL  E.  J 
Chief,  Engl 
Civil  Works 
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FORECASTING  DISTRIBUTION  OF  SEDIMENT  DEPOSITS 
IN  LARGE  RESERVOIRS 
by  Brice  L.  Hobbs 

Reservoir  impoundments  disrupt  the  natural  order  of  the  processes  of 
sediment  transportation  in  streams  and  the  results  range  from  those 
that  are  relatively  insignificant  to  those  where  undesirable  deposi- 
tion is  expected  to  seriously  affect  the  utility  of  the  project  within 
a relatively  short  period  of  time.  It  is  the  purpose  of  this  discus- 
sion to  characterize  some  of  the  problems,  discuss  certain  practical 
considerations  regarding  the  importance  of  forecast  estimates  and  to 
present  some  approaches  for  approximating  future  distributions  of  sedi- 
ment deposits  in  large  artificial  reservoirs. 

RESERVOIR  CLASSIFICATION  ACCORDING  TO  SIZE 


Usually  there  is  a considerable  degree  of  ambiguity  in  designation  of 
a reservoir  as  "large"  or  "small".  The  capacities  of  the  reservoirs 
considered  herein  range  from  about  60,000to  20,000,000  acre-feet  at 
elevations  of  the  spillway  crests.  Perhaps  the  descriptive  terms 
"large"  or  "small"  should  not  be  emphasized  in  discussions  of  reser- 
voir sedimentation;  comparisons  of  relationships  such  as  capacity- 
inflow  ratios  are  more  meaningful.  However,  the  size  of  a reservoir 
(and  therefore,  the  areal  changes  associated  with  pool  fluctuations) 
is  one  of  the  more  important  Influences  affecting  the  distribution  of 
sediment  deposits.  Accordingly,  the  vague  demarcation  between  large 
and  small  reservoirs  is  likely  to  persist. 

DISTRIBUTION  PR0BLE^S 

General  - Information  on  sedimentation  contained  in  many  reservoir  plan- 
ning and  design  reports  suggest  that  those  responsible  attach  little 
importance  to  the  possible  affects  of  adverse  distribution  of  the  sedi- 
ment. Often  times  the  only  information  presented  consists  of  a brief 
statement  to  the  effect  that  estimated  depletion  for  50  or  100  years 
will  represent  only  a small  fraction  of  the  gross  storage.  In  this 
connection  it  should  be  mentioned  that  there  have  been  some  important 
problems  associated  with  local  deposltswhere  depletion  of  gross  stor- 
age is  not  expected  to  be  serious  for  1000  years  or  longer. 

Depletion  of  Storage  Space  - If  volumetric  reductions  of  reservoir  stor- 
age space  allocated  for  various  purposes  represented  the  only  problems 
associated  with  reservoir  sedimentation,  forecast  information  of  fractional 
distributions  of  total  deposits  would  not  serve  any  particularly  useful 
purpose  even  where  rapid  gross  depletion  is  anticipated.  If  such  were 
the  case  it  would  only  be  necessary  to  make  appropriate  reallocations 
as  would  be  indicated  by  periodic  resurveys.  However,  the  significance 


W Sedimentation  specialist,  Corps  of  Engineers,  Office  of  the  Chief 
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of  storage  depletion  and  other  related  problems  depend  generally  upon 
the  average  sedlisentatlon  rates  and  progressive  distribution  o£  deposits. 
On  alluvial  streams,  It  Is  usually  Important  to  have  forecast  estimates 
of  the  probable  distributions  of  deposits  both  with  respect  to  areal 
location  and  volumetric  accumulations  in  various  elevation  zones.  Such 
information  is  useful  in  connection  with  planning. and  design  considera- 
tions to  assure  that  serious  encroachments  upon  space  allocated  for  pur- 
poses other  than  sediment  retention  will  not  occur  during  the  period  used 
for  economic  analysis  of  the  project.  Some  important  sediment  distribu- 
tion problems  are  discussed  further  In  the  following  subparagraphs. 

Aggradation  of  Tributary  Channels  - A reservoir  on  an  alluvial  stream 
is  one  of  the  more  Important  manmade  influences  which  may  affect  channel 
conditions.  The  aggradation  of  channels  which  sometimes  occurs  above 
reservoirs  is  an  extension  of  the  reservoir  sedimentation  processes 
which  may  adversely  affect  drainage  conditions  and  aggravate  flooding 
problems  on  adjacent  lands.  Relatively  small  fractions  of  the  total 
accumulations  are  usually  involved  in  the  aggradation  of  channels  in 
the  reaches  of  reservoir-backwater  influences  above  established  pool 
levels  and  future  dimensions  of  aggraded  channels  cannot  be  accurately 
forecast  by  known  methods. 

Aesthetic  Effects  - Regardless  of  the  need  for  sediment  distribution 
estimates  for  other  purposes.  It  is  occasionally  important  to  foresee 
future  conditions  which  might  be  unsightly  and  therefore  objectionable 
to  people  residing  nearby. 

Depletion  of  Storage  Space  In  Single-Purpose  Reservoirs  - Normally  a 
reasonable  estimate  of  the  total  volume  of  sediment  anticipated  during 
the  period  used  for  economic  considerations  is  all  that  is  necessary 
for  establishing  storage  requirements  in  single-purpose  reservoirs,  and 
advance  information  regarding  the  locations  of  Che  deposits  is  usually 
not  needed.  Exceptions  may  be  found  in  cases  where  substantial  inactive 
storage  is  required  in  reservoirs  operated  primarily  for  power  production. 

Depletion  of  Storage  Space  in  Multiple-Purpose  Reservoirs  - In  cases 
where  sediment  yields  are  appreciable,  advance  information  of  probable 
future  distributions  of  sediment  deposits  is  important  in  connection 
with  planning  and  design  considerations  of  storage  depletion  regardless 
of  the  project  purposes.  Misjudgments  Involved  in  Che  initial  alloca- 
tions of  storage  space  cannot  always  be  satisfactorily  rectified  by 
reallocations  of  space  remaining  at  some  future  date.  For  example,  head 
limitations  might  preclude  lowering  Che  elevation  of  the  minimum  power 
pool . 

Depletion  of  Space  Where  Water  is  Stored  for  Recreational  Purposes 
Recently,  there  has  been  a rapid  increase  in  demands  for  storage  of 
water  for  recreational  activities  in  artificial  lakes.  The  needs  are 
usually  satisfied  by:  use  of  water  stored  primarily  for  other  purposes; 
provisions  for  perpetual  storage  of  a given  volume  of  water  regardless 
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of  pool  elevation;  specific  allocations  of  storage  below  a given  pool 
elevation;  or  arrangement  for  regulation  so  as  to  provide  for  a minimum 
pool  having  storage  not  exceeding  that  provided  for  conservation  and 
the  undepleted  space  initially  reserved  for  sediment.  There  is  general 
agreement  regarding  the  importance  of  recreational  needs,  therefore, 
the  problems  that  may  be  expected  to  result  from  unfavorable  sediment 
distributions  should  be  recognized.  For  example,  a plan  to  continuously 
provide  a small  pool  of  fixed  volume  in  the  lowest  elevation  zone  of 
remaining  space  may  become  completely  unsatisfactory  for  the  planned 
activity  relatively  early  in  the  life  of  the  project.  Also,  decisions 
are  often  made,  after  completion  of  the  design  stage  and  without  benefit 
of  additional  engineering  study,  to  regulate  a reservoir  so  as  to  uti- 
lize space  reserved  for  sediment  deposits  for  recreational  or  other  con- 
servation purposes.  In  such  cases  there  is  no  opportunity  for  changing 
the  total  storage,  therefore,  the  effects  of  the  change  on  sediment 
distribution  expected  to  result  from  the  change  in  regulation  procedures 
should  be  carefully  examined. 

Shore  Erosion  - Shore  erosion  and  bank  caving  processes  frequently 
create  beach  and  boat  harbor  problems.  Movement  of  material  by  these 
processes  may  cause  an  exchange  of  storage  space  between  elevation-zones 
or  a net  storage  loss  or  both. 

Utilization  of  Delta  and  Backswarop  Areas  - Interests  opposing  the  con- 
struction of  reservoirs  frequently  cite  sedimentation  as  one  of  the 
horrible  consequences  of  these  developments  and  the  general  public  is 
led  to  believe  that  the  results  are  always  entirely  bad.  Actually, 
the  program  for  wildlife  propogation,  by  the  U.  S.  Fish  and  Wildlife 
Service,  in  the  delta  areas  of  Denison  Reservoir  is  reported  to  be 
quite  successful.  This  represents  a type  of  planning  problem  that  has 
not  had  proper  consideration  in  the  past. 

FACTORS  AFFECTING  DEPOSIT  ACCUMULATIONS 

General  - The  factors  involved  in  reservoir  sedimentation  processes  are 
numerous.  The  influences  most  frequently  mentioned  in  qualitative  dis- 
cussions!.^ are;  (1)  Reservoir  size  and  shape;  (2)  Sediment  quantities 
and  characteristics;  (3)  Sediment  sources;  (4)  Progressive  vegetative 
growth  on  frequently  exposed  deposits;  (5)  Consolidation  of  deposits; 

(6)  Magnitudes,  frequency  and  sequences  of  hydrologic  events;  and  (7) 
Reservoir  regulation  practices.  These  factors  and  other  influences 
cofflingle  in  ever  changing  combinations  to  produce  the  distribution  of 
deposits  at  any  given  time. 

Dominant  Factors  - As  indicated  above,  the  oistribution  as  well  as  the 
quantities  o^  sediment  involved,  are  sensitive  to  numerous  factor  combi- 
nations which  include  unpredictable  sequences  of  flood  events  and  pool 
elevations  coincident  with  high  sediment  inflows.  Regulation  is  one  of 
the  dominant  influences  affecting  deposition  in  reservoirs.  In  fact, 
the  first  five  factors  listed  above  are  governed  in  some  degree  by  pool 
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fluctuations.  The  charts  and  diagrams  shown  on  Incls  1 and  2 illustrate 
the  effects  of  regulation  on  sediment  deposition  in  a large  multiple- 
purpose  reservoir. 

Incl,  No,  1 shows  estimat€>d  sand  loads  delivered  by  a design  flood  which 
was  developed  for  a large  reservoir  on  an  alluvial  stream.  The  flood 
was  assumed  to  have  started  at  a time  when  the  pool  level  was  at  the 
bottom  of  the  flood  control  pool.  Attention  is  directed  to  the  cumula- 
tive sand  curve  which  shows  that  97  percent  of  the  sand  transported  by 
the  flood  would  be  delivered  to  the  reservoir  before  the  maximum  pool 
elevation  was  reached.  The  sand  load  graph  is  based  upon  an  average 
rating  curve  and  therefore,  has  the  same  shape  and  peak  time  as  the 
discharge  hydrograph;  actually  maximum  suspended  sediment  concentra- 
tions quite  often  precede  peak  discharge  rates  and  under  such  condi- 
tions, the  cumulative  sand  curve  would  be  in  a position  to  the  left  of 
that  shown  on  Incl  1. 


Coincidental  values  of  inflow  and  pool  elevations  from  graphs  on  Incl  1 
were  used  to  develop  the  general  illustration  shown  on  Incl  2.  The  locus 
of  the  upstream  limits  of  backwater  effects  in  this  hypothetical  situa- 
tion further  demonstrates  why  most  of  the  sediment  delivered  to  a large 
flood-control  reservoir,  by  any  given  flood,  is  transported  to  areas 
below  the  higher  pool  elevations  attained.  Also,  it  affords  some  insight 
to  tendencies  for  material  previously  deposited  to  be  redistributed. 
Reference  Nos.  \J  through  5/  contain  considerable  additional  information 
regarding  the  manner  in  which  sediment  is  transported  into  and  deposited 
in  a reservoir. 

THODS  FOR  ESTIMATING  DISTRIBUTION  OF  SEDIMENT  DEPOSITS 

General  - There  are  several  methods  described  in  the  literature  for 
estimating  future  distributions  of  sediment  deposits.  For  purposes  of 
this  discussion,  it  is  considered  satisfactory  to  classify  these  methods 
as:  (1)  analytical  methods—^*  2.^*  — ^ which  utilize  procedures  based 

upon  theoretical  concepts  of  hydraulics  and  sediment  transport  and  (2) 
purely  empirical  methods2^.  All  methods  are  subject  to  limitations 
imposed  by  necessary  simplifications  involving  uncertain  assumptions 
regarding  sequences  of  significant  hydrologic  events  and  lack  of 
accountability  for  other  dominant  influences  which  help  to  produce 
individual  deposition  patterns.  For  example,  large  tributaries  occasion- 
ally transport  5 to  10  times  the  average  annual  sediment  load  in  a 10- 
day  period  and  the  range  of  extremes  in  "low-order"  tributaries  is  often 
much  greater.  There  is  no  way  of  predicting  when  during  the  project 
life,  extreme  events  may  occur.  The  uncertainties  and  difficulties  not 
withstanding  it  is  Important  to  make  the  best  possible  distribution 
determination  commensurate  with  probable  seriousness  of  anticipated 
problems  and  practical  considerations  of  the  purpose  of  the  estimate, 
available  data  and  time  allowed  for  the  study. 

Choice  of  Methods  - Where  sediment  deposition  is  expected  to  have  a 
major  effect  upon  the  design  and  operation  of  a reservoir  project,  it 
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is  prudent  to  use  more  than  one  approach  so  that  the  results  of  some- 
what independent  determinations  can  be  used  as  guidance  for  judgment 
in  allowances  for  conservatism.  In  many  project  investigations  it  is 
impracticable  to  attempt  detailed  analyses,  to  consider  the  probable 
effects  of  many  individual  influences,  because  of  insufficient  data  and 
a number  of  other  limitations.  With  present  knowledge,  the  writer 
believes  that  in  most  cases  the  results  obtained  by  empirical  methods 
are  sufficiently  reliable  for  engineering  purposes.  However,  there  are 
many  cases  where  thoroughly  rigorous  sedimentation  analyses  are  impera- 
tive. Dardanelle  Reservoir—^  is  one  example.  Also,  there  are  cases 
where  serious  tributary  channel  and  flood-plain  aggradation  may  be 
expected  to  accompany  delta  developments.  In  some  of  these  situations, 
it  may  be  found  reasonable  to  use  empirical  methods  to  establish  the 
initial  condition  assumed  for  the  starting  point  in  an  analytical  study 
of  patterns  of  continued  deposition. 

Empirical  Methods  - In  practice,  empirical  methods  are  being  used  in 
most  instances  where  sediment  distribution  estimates  are  required. 

These  methods  have  the  advantage  of  simplicity  but  sacrifice  considera- 
tion of  some  important  influences  and  when  followed  implicitly,  these 
methods  often  produce  misleading  results.  This  is  true  of  the  method 
presented  herewith  as  well  as  all  of  those  described  in  literature 
known  to  the  writer.  All  depend  upon  the  same  basic  requirements  for 
estimates  of  total  sediment  loads,  average  trap  efficiencies  and  gross 
volumes  of  sediment  trapped  during  the  period  under  consideration. 

None  delineate  developments  at  individual  tributaries.  The  methods 
described  by  Borland  and  Millev^^  are  probably  most  widely  used  in  the 
United  States;  it  is  believed  that  the  principal  weakness  of  these 
methods  is  that  no  particular  consideration  is  given  to  the  sediment 
characteristics  or  the  operational  effects  of  the  reservoirs.  The 
method  presented  in  Incl  3 also  has  obvious  weaknesses  but  considera- 
tion is  given  to  sediment  characteristics  and  regulatory  influences. 

.Pool-Elevation  Duration  Method  - This  empirical  method  attempts  to 
account  for  the  influences  of  reservoir  regulation  inherent  in  a pool- 
elevation  duration  curve,  general  effects  of  the  fraction  of  sand  mate- 
rials involved  and  the  size  and  shape  of  the  reservoir.  The  approach 
is  based  upon  the  idea  that  pool  elevation  and  the  size  characteristics 
of  the  sediment  are  two  of  the  most  important  factors  Influencing  deposi- 
tion in  given  elevation  zones.  It  also  embodies  the  thesis  that:  (1) 
over  a long  period  of  time,  sediment  delivered  by  medium  and  moderate 
floods  will  establish  some  statistical  order  of  coincidence  with  pool 
elevations  between  the  maximum  and  minimum  and  (2)  that  regulation  of 
the  rare  floods  (and  therefore,  the  distribution  of  sediment  deposited 
in  the  higher  elevation  zones)  will  be  similar.  This  suggests  that 
there  may  be  some  reasonably  definable  relationships  between  duration 
of  a given  pool  and  the  amount  of  sediment  that  will  be  deposited  above 
and  below  the  elevation  of  that  pool.  Step  procedures  for  one  simple 
approach  are  presented  in  Incl  3.  This  index  method  is  based  upon 
limited  data  and  its  reliability  for  general  application  will  require 
much  more  testing.  There  are  several  refinements  that  appear  worthy  of 
further  study. 

I 
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CONCLUDING  OBSERVATIONS 

1.  The  Reservoir  Sediment  Data  Summaries  assembled  and  published  by 
the  Committee  on  Sedimentation,  Water  Resources  Council,  represent  the 
most  complete  source  of  available  information  on  sediment  distribution 
but  relatively  few  of  these  data  summaries  contain  information  regard- 
ing the  sediment  size  characteristics.  Also,  many  of  the  formal  sedi- 
ment survey  reports,  which  present  these  summaries  and  other  information, 
contain  little  useful  data  on  particle  sizes. 

2.  Some  thought  might  well  be  given  to  the  utility  of  certain  informa- 
tion presented  in  accordance  with  the  instructions  issued  for  preparing 
certain  items  of  the  Reservoir  Sediment  Data  Summaries,  For  example, 
in  some  large  reservoirs  the  capacity  at  the  elevation  of  the  spillway 
crest  elevation  may  bear  little  or  no  relationship  to  the  trap-efficiency; 
some  other  capacity  value  would  often  be  more  meaningful  in  capacity- 
inflow  ratios. 

3.  A great  mass  of  reservoir  sediment  data  has  been  collected  on  a 
routine  basis  a,9C9rd)ing  to  somewhat  arbitrary  observation  schedules. 
Analytical  studies  of  these  data  to  establish  useful  relationships  for 
general  application  and  to  identify  data  collection  deficiencies  have 
lagged  far  behind  the  data  collection  programs.  It  is  believ&d  that 
some  of  the  money  being  used  for  surveys  should  be  diverted  for  funding 
some  intensive  analytical  studies  of  the  information  now  available. 

4.  Without  the  aid  of  electronic  computers,  it  would  generally  be 
impracticable  to  attempt  rigorous  analytical  investigations  aimed  at 
predicting  future  distributions  of  sediment  deposits  in  large  reservoirs 
having  a large  number  of  tributaries.  Important  advances  have  been  made 
in  computer  applications  in  hydrologic  analyses  relating  to  flood  inflows 
and  regulation  of  reservoirs.  Also,  considerable  progress  has  been  made 
in  the  use  of  electronic  computers  in  sediment  discharge  determinations.^^ 
The  possibilities  for  Integrating  these  techniques  appears  to  offer  the 
greatest  promise  for  development  of  improved  analytical  procedures  for 
investigating  reservoir  sediment  distributions. 


Incls 

1.  Illustration-Coincidental 
Sand  Inflows  & Reservoir 
Pool  Elevations 

2.  Illustration-Locus  of  Up- 
stream Limits  of  Reservoir 
Effects 

3.  Paper-"Forecasting  Sediment 
Distribution  in  Large  Reser- 
voirs" 
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IN  LARGE  RESERVOIRS 


POOL- ELEVATION  DURATION  METHOD 


INCLOSURE  3 ^ 

FORECASTING  SEDIMENT  DISTRIBUTIONS 


1 . Required  Information. 

a.  Pool  elevation  charts  developed  in  connection  with  operation  studies. 

b.  Reservoir  capacity  table.  (Table  No.  1) 

c.  Estimate  of  total  volume  of  sediment  expected  to  accumulate  in 
reservoir  during  period  under  consideration. 

d.  Estimate  of  sand  fraction  of  total  deposit. 

2.  Estimating  procedure  (use  Ft.  Peck  Reservoir  data  for  explanation). 

a.  Develop  the  pool  elevation  duration  curve  from  pool  elevation 
graphs  (Curve  1,  Fig.  1). 

b.  Plot  first  differences  of  capacity  for  depth  increments  (five- 
foot  increments)  on  log-log  paper.  (Fig.  2) 

c.  Draw  estimated  distribution  curve  on  Fig.  3 with  position  based 
on  the  sand-percent  scale  (four  percent  for  Ft.  Peck)  and  judgment  of  the 
plotting  positions  of  points  determined  from  measurements  at  other  reser- 
voirs. The  right  envelope  position  was  selected  because  of  the  low  per- 
centage of  sand  and  the  large  capacities  of  pools  in  the  operating  range 
(from  about  110,000  to  19,000,000  acre-feet).  (Sand  scale  shown  on  Fig.  3 
is  explained  in  paragraph  3 below.) 

d.  Prepare  Table  No.  2 as  follows; 

(1)  Tabulate  time  durations  (10  percent,  20  percent  ...  95  percent 
and  100  percent)  in  column  No.  1. 

(2)  Tabulate  pool  elevations  corresponding  to  the  durations  in 
column  No.  2 (obtain  values  from  Curve  No.  1,  Fig.  1). 

(3)  Tabulate  first  differences  of  capacity  (obtained  from  Fig.  2) 
in  column  No.  4. 


y Incl  3 with  paper  "Forecasting  Distribution  of  Sediment  Deposits  in 
Large  Reservoirs,"  by  Brice  L.  Hobbs,  10  Feb  69. 
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(4)  Conipuit'  ratios  (first  differences  of  capacity  divided  by 
the  first  dif fi-rence  i^f  capacity  corresponding  to  Che  pool  elevation 
that  is  exceeded  only  five  percent  of  the  time)  and  tabulate  in  column 
No.  5. 

(5)  Enter  the  chart  (Fig.  3)  with  ratios  from  column  No.  5 
and  tabulate  the  corr»'sponding  values  of  percent  of  total  sediment 
deposits  in  column  No.  6.  These  values  represent  the  estimated  distri- 
bution of  deposits.  Measured  values  are  tabulated  in  column  No.  7. 

3.  The  sand-perci-n t scale  on  Fig.  3 is  plotted  from  values  taken 
from  Fig.  4 which  is  a correlation  of  percentages  of  sand  with  total 
depos i ts . 

BRICE  L.  HOBBS 

Hydrology  6 Hydraulics  Branch 
Engineering  Division,  CW.  OCE 
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TABLE  NO.  1 


FORT  PECK  RESERVOIR 
CONDENSED  AREA-CAPACITY  TABLE 
(Based  on  1961  Aggradation  Survey) 


ELEV 

DEPTH 

AREA 

CAPACITY 

II.  s . 1 . ) 

(Ft.  ) 

(Acres  ) 

(Ac re-Feet ) 

2033 

0 

0 

0 

2035 

2 

103 

113 

2040 

7 

402 

1,214 

2045 

12 

1,075 

5,002 

2050 

17 

1,652 

11,109 

2055 

22 

2,305 

21,423 

2060 

27 

4, 149 

36,870 

2070 

37 

10,672 

106,662 

2080 

47 

16,714 

245,371 

2090 

57 

22,966 

440,692 

2100 

67 

29,732 

702,113 

2110 

77 

38,458 

1,042,665 

2120 

87 

50,560 

1,484,307 

2130 

97 

61,391 

2,044,261 

2140 

107 

71,243 

2,709,084 

2150 

117 

81,944 

3,474,396 

2160 

127 

92,712 

4,346,056 

2170 

137 

106,393 

5,335,418 

2180 

147 

122,028 

6,485,415 

2190 

157 

936,912 

7,777,395 

2200 

167 

152,792 

9,222,634 

2210 

177 

170,021 

10,839,099 

2220 

187 

187,829 

12,625,547 

2230 

197 

206,874 

14,600,015 

2240 

207 

226,827 

16,771,900 

2250 

217 

246,919 

19,138,48? 

2260* 

227 

270,200 

21,704,68*. 

^Extrapolated  above  elevation  2250 


ite  ratio  is  the  increnental  first  RESERVOIR  S' 
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TABLE  NO.  2 

ESTIMATE  OF  DISTRIBUTION  OF  SEDIMENT  DEPOSITS 
IN  FORT  PECK  RESERVOIR 
MISSOURI  RIVER,  MONTANA 


FIRST  DIFF  SEDIMENT 


POOL  ELEV. 
DURATION 
(Percent 
of  Time)J/ 

ELEV 
(Ft  MSL) 

DEPTH 

(Ft) 

OF  CAPACITY 
(Ac-Ft/5-Ft 
Depth 
Increment) 

RATloi/ 
(Col.  4 + 
1.125.000) 

DISTRIBUTION 

ESTIMATED 
(f  %) 

MEASURED 

(i  %)y 

(1) 

(2) 

(3) 

(4) 

(5) 

(6) 

(7) 

10 

2,117 

84.0 

236,000 

ty 

0.27^  . 

30.0 

32.0 

20 

2,116.5 

143.5 

580,000 

0.52 

65.0 

78.5 

30 

2,195 

162.0 

722,600 

0.64 

79.0 

92.0 

40 

2,208 

175.0 

828,000 

0.74 

88.0 

95.5 

50 

2,212 

179.0 

862,000 

0.77 

90.0 

97.8 

60 

2,218 

185.0 

915,000 

0.81 

94.0 

98.8 

70 

2,225 

192.0 

987 , 000 

0.88 

97.0 

99.5 

80 

2,230 

197.0 

1,030,000 

0.92 

98.5 

99.8 

90 

2,236 

203.0 

1,090,000 

0.97 

99* 

95 

2,240 

207.0 

1,125,000 

1.00 

99.5 

99.95 

100 

2,248 

215.0 

1,200,000 

1.07 

100.0 

100.0 

i^ercent  of 

time  pool 

was  at 

or  below  corresponding  elevation  shown  in 

Co 1 uan  No . 2 . 

Clo  is  I.O  at  the  95  percent  pool. 

i(/a  lues  from  Item  No.  26  of  Reservoir  Sediment  Data  Summary  (See  Curve  3 
Fig  1). 
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